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Infrared Spectra of Pentachlorofluoroethane, 1, 2- 
Dichlorotetrafluoroethane, and 1-Bromo-2-Fluoroethane 


Delia Simpson ' and Earle K. Plyler 


The infrared absorption spectra of pentachlorofluoroethane and 1,2-dichlorotetrafluro- 


ethane have been measured from 2 to 40 yu. 


25 to 40 uw, which had not 


been previously 


Several bands were found in the region from 


investigated. The spectrum of 1-bromo-2- 


fluoroethane has also been measured in the spectral region from 2 to 40 uw, and many bands 


were observed. 
molecules are given. 


much interest has been shown 
of fluorinated hydrocarbons. 
* for example, have given the 


In recent years 
in the spectroscopy 
Glockler and Sage,’ 


Raman spectra of such compounds and Smith, Niel- | 


sen, and coworkers **® have given the infrared and 
Raman spectra of many fluorinated ethanes, along 
with frequency assignments for some of them. In 
this report we present spectra and assignments for 
three fluorinated hydrocarbons, pentachlorofluoro- 
ethane (CCLCCLF), 1,2-dichlorotetrafluoroethane 
(CCIF.CCIF,), and 1-bromo-2-fluoroethane (CH, 
BrCH,F). The Raman spectrum of CCIF,CCIF, has 
been reported by Glockler and Sage (see footnote 2) 
and the infrared spectra of CClCCIL,F and CCIF, 
CCIF, have been reported by Smith and Nielsen (see 
footnotes 4 and 5). We have extended the infrared 
spectra of the last two named molecules to 40 » and 
present the spectrum of the third in the 2 to 40 u 
range. 

The spectra of these compounds were obtained on 
a Perkin-Elmer model 12C_ spectrometer. The 
prisms used in the various spectral ranges were: | 
to 5 uw, lithium fluoride; 2 to 16 uw, sodium chloride; 
14 to 24 w, potassium bromide; and 23 to 40 u, thal- 
lium bromide-iodide. 

Figure 1 shows the spectrum of CCI,CCLF. It 
was obtained by comparing the absorption of the 
compound with the absorption of air. These obser- 
vations were then reduced to the percent transmit- 
tance of the compound. Figures 2 and 3 present the 
spectra of CCIF,CCIF, and CH,BrCH,F obtained in 
a similar manner. These compounds were examined 
in various states and under various experimental con- 
ditions, the conditions being described on the figures 
and in the captions. 

Because these molecules do not form a closely re- 
lated series of substituted ethanes, it would not be 
expected that their spectra would show many simi- 
larities. A general feature of the three spectra is 
that most of the intense bands fall between 8 and 20 u. 

In tables 1 through 3 are given the observed infra- 
red bands, the available Raman data, and these fre- 


quency assignments for the three molecules. These 
Newnham College, Cambridge University, Cambridge, England 
(+. Glockler and C. Sage, J. Chem. Phys. 8, 291 (1940 
* (i. Glockler and C. Sage, J. Chem. Phys. 9, 387 (1941 
* D.C. Smith, J. Reed Nielsen, L. H. Berryman, H. H. Claassen, and R L 
Hudson, NRL Report 3567 
). C. Smith, M. Alpert, R. A. Saunders, G. M. Brown, and N. B. Moran, 


NRL Report 3924 


A list of the observed bands and the frequency assignments for the three 


assignments were made on the basis of analogy with 
the assignments of closely related structures which 
have previously been studied. Table 4 compares 
these assignments for CH,BrCH,F with the wave 
number ranges for the various vibrations of three 
molecules. (CH,Cl)., (CH,Br)., and CH,BrCH,Cl. 
These values are taken from the assignments of J. U. 
Brown and W.S. Sheppard, who have kindly allowed 
the use of the results of their work, which now await 
publication. These four molecules listed in table 4 
are expected to exist in two isomeric forms, trans and 
gauche. In this table are given the differences in 
going from the gauche to the trans form for individual 
vibrations. The assignments for CCI,CCLF and 


CCIF,CCIF, were made in the same way by com- 
paring their vibrations with related molecules. 
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Figure 1. Infrared absorption spectrum of CCIsCChF. 


The spectrum of this substance, a solid at room temperature, was obtained by 
making a saturated solution in CS». Insert A was obtained by mixing 1 part of 
the saturated solution with 4 parts of CS»; insert B was obtained by mixing I 
part of the saturated solution with 20 parts of CS). Dotted lines represent regions 
of intense atmospheric absorption and of absorption of the solvent 


Way 


Figure 2. Infrared absorption spectrum of CCIF,CCIPF». 
his spectrum is of the vapor stat« The pressure was atmospheric, except for 
the inserts which were obtained by removing some of the vapor. The cell thick 
ness employed for the different regions are marked on the figure Dotted lines 
ire used in regions of intense atmospheric absorption 
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The Raman frequencies are those of Glockler and Sage 
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Modified 1913 Reference Tables for 
Iron-Constantan Thermocouples 


Robert J. Corruccini and Henry Shenker 


An investigation was made of the characteristics of iron-constantan thermocouples 
typical of those being supplied to reproduce the standard temperature-electromotive force 


relationship that was established commercially in 1913. 


Using the calibration of the thermo- 


couple that most nearly matched the 1913 reference table over the range to which iron- 


constantan thermocouples usually are limited (32 


to 1,400° F), a slightly different tempera- 


ture-electromotive force relationship was derived which is more nearly realizable by com- 


mercially available materials. 


tables covering the range, —320° to + 1,600 


The new relationship is presented in four mutually consistent 


(—196° to +871° C), in both Fahrenheit and 


Celsius units and with both temperature and electromotive force as arguments. 


1. Introduction 


The application of thermocouples to pyrometry has 
heen accompanied by the development of indicating 
instruments (millivoltmeters, potentiometers) with 
scales calibrated as to indicate temperature 
directly. Such seale calibration requires a knowl- 
edge of the temperature-emf relationship of the 
associated thermocouple. In general, the thermo- 
couple is required to be readily replaceable, whereas 
the associated indicator is relatively long-lived. It is 
not practical to alter the scale of the indicator to 
allow for changes in the characteristics of the thermo- 
couples used with it. Consequently, each manu- 
facturer of a given type of thermoelectric pyrometer 
calibrates the seales of the indicator in accordance 
with a definite temperature-emf relationship (refer- 
ence table) and undertakes to furnish replacement 
thermocouples that satisfy this relationship within 
certain narrow limits. 

Several reference tables have existed in the past for 
each of the common thermocouples, each supplier 
adopting, for a given type of thermocouple, a table 
or sometimes two) that was thought to represent 
closely the thermocouples supplied by that company. 
Thus, a partial survey in 1937 showed the existence 
of four reference tables for iron-constantan thermo- 
couples, differing from each other by as much as 7 
percent. 

In the vears 1933 to 1938, the National Bureau of 
Standards published reference tables for all of the 
common thermocouples [{1, 2, 3].'* These tables 
were based on laboratory calibrations of numerous 
representative samples of these thermocouples. The 
tables in RP530 for platinum-platinum rhodium and 
those in RP1080 for copper-constantan were based 
on older tables and soon were widely accepted. The 
tables for chromel-alumel in RP767 likewise were 
widely accepted. Here, the situation was simplified 
by the fact that all chromel-P and alumel originates 
with one producer who prepares these alloys to a 
single set of specifications. 


so 


Figures in brackets indicate the literature references at the end of this paper. 
? Most of the same tables are now available in expanded and revised form in 
NBS Circular 508, Reference Tables for Thermocouples, by Henry Shenker, 
John I. Lauritzen, Jr., and Robert J. Corruccini, issued May 7, 1951. 
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On the other hand, the tables for iron-constantan 
thermocouples published in RP1080 were the result 
of an investigation that began after another Govern- 
ment agency had noted that the NBS calibrations of 
iron-constantan thermocouples being purchased in the 
early 1930's did not match closely any existing refer- 
ence table. The reference table in RP1080 repre- 
sented an average of the commercially available 
thermocouples tested at that time and, hence, did not 
resemble closely any of the older tables for this type 
of thermocouple. RP1080 showed the manufae- 
turers the necessity of either (1) altering the scales 
of new indicators, and of some existing ones, to con- 
form to a temperature-emf relationship that was 
representative of the thermocouples actually being 
produced (RP1080), or (2) to develop new thermo- 
couple materials that would have, within the desired 
tolerances, the temperature-emf relationship em- 
bodied in the scales of the existing indicators. Most 
of the manufacturers chose the latter course. Users 
with large investments in installed indicators were 
satisfied with this development, and the RP1080 
table has not displaced the previous tables to any 
great extent, although it is the only table used for 
iron-constantan thermocouples in the field of military 
aircraft. 

A further difference between the situations for 
iron-constantan and the other common thermo- 
couples is related to the manner of procurement. 
Thermoelectric platinum and platinum-rhodium are 
prepared to specifications having the end use in view. 
Constantan and the purer grades of copper likewise 
are prepared with careful attention to the electrical 
properties of the product. On the other hand, such 
a small proportion of the total production of iron 
goes to thermocouples, that large producers of iron 
cannot be interested in controlling the properties of 
their product so as to suit it to this use. The manu- 
facturer of pyrometers thus finds it more economical 
to select from the available production of commercial 
iron certain lots that happen to be suitable for 
thermoelectric use, rather than to have a suitable 
material prepared to specifications at a premium 
price by a producer of specialty metals. In classify- 
ing lots of iron, chemical tests are helpful, but the 
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final test of the suitability of a lot must be a thermo- 
electric one. 

In an effort to promote uniformity, steps were 
initiated in 1948 in the Recorder-Controller Section 
of the Scientific Apparatus Makers of America 
(SAMA), representing a number of pyrometer manu- 
facturers, directed toward the adoption of the tem- 
perature-emf relationship known as the “1913” table® 
as a tentative standard for iron-constantan thermo- 
couples over the range, 32° to 1,400° F. This table 
was chosen because it is used more than any other. 
It is apparently based on a single calibration of a 


thermocouple (no longer in existence) that was typi- 


cal of the iron-constantan of forty vears ago. How- 
ever, experience has shown that it is apparently not 
possible to produce a thermocouple, using present- 
day commercial iron, that will correspond precisely 
to the 1913 table. In addition, it was desired that 
the standard table should extend to subzero tempera- 
tures, whereas the original 1913 table did not ex- 
tend below 32° F 

It appeared that experimental data on the char- 


acteristics of thermocouples of the type that repro- 


duce the 1913 table were required in order to (1) 
judge the adequacy of the 1913 table as a permanent 
standard for SAMA, or, if it should prove desirable, 


table, is identical with 


rhis table, sometimes called the Leeds and Northrup 
50, (1926) and in Tech, Pap 


the tables given in International Critical Tables, 1 


Corresponding values of temperature and electromotive force for iron wires versus Pt 





to make adjustments in the table, and (2) to obte 


data at low temperatures. Accordingly, a resea: 
project for this purpose was established at the NJ 
and was supported in part by a grant from SA\M 
and in part by funds of the Bureau. This paper p; 
sents the results of this investigation. 


2. Experimental Procedure 


Each member firm of the SAMA that suppli: 
iron-constantan thermocouples was invited to subm 
samples for investigation at the Bureau. These 
thermocouples were to be of commercial materials 
that were selected to match closely the 1913 table 
over the range 32° to 1,400° F. Of the 8 eligible 
member concerns, 4 responded. Of these, 3  fur- 
nished 1 lot of wire each, while the fourth concern 
furnished 5 different lotsof wire. The samples were 
designated by a letter representing the supplier, fol- 
lowed by a number representing the sample, e. 
A-1. The wire size was approximately No. 8 AWG 
Information received from some of the suppliers 
indicated that the iron samples came from at least 
four different basic sources of ingot iron. 

The thermoelectric measurements were made in 
two parts: (a) From 32° to 1,800° F the wires were 
calibrated in a tube-furnace which has already been 
described [4, 5]. For convenience, the iron wires 


to 


Reference junctions at 32° F 


Thermal electromotive force (absolute millivolts 


BS 14, 306 (1920-21) T170, Table IV, Col. I 
TABLE 1. 
Tempera- 
ture (Int 
1948) VI om D1 
Pr 
319.5 2. 54 2. 4 2. 54 
297. 4 2. 50 2. 60 2. 50 
147. 1 1. 65 1. 72 1. 66 
103. 0 1. 28 1. 33 1. 28 
61.3 0. 91 0. 95 0. 91 
27. 6 0. 58 0. 61 0. 59 
32 0. 00 0. 00 0. 00 
100 +0. 6S 0. 71 0. 68 
200 1. 66 1. 74 1. 67 
300 2. 59 2. 72 2. 60 
100 3. 44 3. 62 3. 46 
500 4. 20 4.42 4. 23 
600 1. 86 5. 18 4. 90 
700 5. 45 5. 76 5. 49 
SOO 5. 99 6. 34 6. 04 
900 6. 52 6. 90 6. 57 
1, 000 7. 06 7. 47 > a 
1, 100 7. 66 8. 09 7.72 
1, 200 &. 34 8. 79 &. 41 
1, 300 9. 14 9. 60 9. 21 
1, 400 10. 07 10. 52 10. 15 
1, 500 11. 06 11. 44 11.17 
1, 600 12. 01 12. 26 12. 19 
1, 700 12. 77 12. 98 12. 94 
1, 800 13. 57 13. 79 13. 75 


for samples 


F-1 F-3 kF-4 F-5 F-6 
2 65 2. 22 2. 43 2. 35 2. 29 
2. 60 2. 19 2. 40 2. 32 2. 26 
1. 70 1. 49 1. 61 1. 56 1. 52 
1. 31 1. 17 1. 25 1. 21 1. 18 
0. 93 0. 83 0. 89 0. 86 0. 84 
0. 59 0. 54 0. 58 0. 55 0. 54 
0. 00 0. 00 0. 00 0. 00 0. 00 
+0. 68 +0. 64 +0. 67 +0. 65 +0. 64 
1. 66 1. 59 1. 65 1. 60 1. 57 
2. 56 2. 50 2. 57 2. 50 2. 46 
3. 38 3. 35 3. 42 3. 32 3. 27 
4.09 1.13 4.18 4. 06 4. 00 
4. 70 4.82 4. 84 4.70 4. 64 
5. 25 5. 45 5. 43 5. 28 5. 21 
5. 74 6. 03 5. 97 5. 81 5. 73 
6. 23 6. 60 6. 49 6. 32 6. 24 
6. 74 7.17 7. 03 6. 86 6. 77 
7. 31 7. 80 7. 63 7. 45 7. 36 
7. 97 8. 51 8. 31 8. 12 8. 04 
8. 76 9. 33 9. 11 8. 92 8. 83 
9. 68 10. 25 10. 04 9. 84 9. 76 
10. 73 11. 21 11. 08 10. 85 10. 76 
11. 83 12. 11 12.19 11. 86 11. 71 
12. 76 12. 89 13. 20 12. 70 12. 61 
13. 61 13. 71 14. 00 13. 51 13. 44 
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were calibrated against a standard iron wire, and the 
copstantan wires against a standard constantan 
wire. The standard iron and the standard con- 
stantan wires had previously been calibrated against 
Pt 27, the platinum thermoelectric standard of the 
NBs, and they were rechecked at frequent intervals 
during each run against the platinum element of the 
platinum versus platinum-10 percent rhodium ther- 
mocouple used to measure the temperature. The 
test wires were calibrated three at a time. The 
measuring junctions of the test wires, the standard 
wire of the same material, and the platinum versus 
platinum-rhodium thermocouple used to measure 
the temperature were welded together. Because of 
the relative instability of the iron and constantan at 
temperatures above 1,400° F, measurements in this 
range were made against a platinum standard only. 
The average rate of rise of temperature in the region 
above 1,400° F was about 8 deg F per minute. All 
measurements were made at 100 deg F intervals on 
rising temperature. After each run, the heated 
section of the iron or constantan standard was cut 
off and discarded. The methods used were similar 
to those described in [4, p. 262-5, 277]. The 
measurements are considered to be accurate within 
the equivalent of ldeg F. (b) From +-32° to— 319.5 

F, the thermal emf of each wire against a copper 
reference wire was determined in stirred liquid baths. 
The copper reference wires were calibrated against 





a platinum standard. Down to —147° F the eryo- 
stat described by Scott [6; 5, p. 206] was used. Ob- 
servations at two lower temperatures were made by 
using baths of liquid oxygen and nitrogen stirred by 
a stream of the same gas. Temperatures in the 
baths were determined with a platinum resistance 
thermometer. Sections of wire that had not been 
previously heated were used for the low-temperature 
measurements. The limit of error in this range is 
considered to be the equivalent of 4 deg F. 

In all the measurements, the reference junctions 
were kept at 32° F in a mixture of ice and water, 
An atmosphere of air surrounded the thermocouples 
at all times. The calibrations of the individual 
wires relative to Pt 27 at the various calibration 
temperatures were compiled and are given in tables 
1 and 2. The data in table 3 for the complete 
thermocouples were obtained by combination of the 
data in tables 1 and 2. 

The deviations (expressed in deg F) of the various 
thermocouples from the 1913 table are shown in 
figure 1. The deviation of the RP1OS80 table is also 
shown. Because the A-1 thermocouple gave the 
nearest over-all match to the 1913 table in the range 
up to 1,400° F, the calibrations of the iron and con- 
stantan elements of this thermocouple were used as 
base lines for figures 2 and 3, respectively, in which 
are shown the differences among the various indi- 
vidual wires. The curves marked “RP1080" repre- 


PaBLe 2. Corresponding values of te mperature and electromotive forces for constantan wires versus Pt-27 
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tron from iron sample A-1. 


F — » 1913 RP 1080 
5 rs r¢ 
F-4 F-5 F , table table 
7. 73 7. 66 7. 60 
7. 44 7. 37 7. 31 7. 83 
4.59 4.55 $51 1. 80 
3. 55 3. 52 3. 49 3. 71 
2. 53 2. 50 2. 48 2. 63 
1. 63 1. 61 1. 60 1. 69 
0. 00 0. 00 0. 00 0. 00 0. 00 
+ 1.95 1. 93 +1. 92 1. 96 2. 00 
1. 93 1.88 1. 86 4. 92 5. 03 
7. 99 7. 91 7. 88 7. 95 8. 15 
11. 09 11. 00 10. 96 11. 02 11. 25 
14. 19 14. OS 14. 04 14. 10 14. 35 
17. 27 17. 15 7. 10 17. 17 17. 43 
20. 35 20). 23 20. 18 20. 24 20. 49 
23. 42 23. 29 23. 24 23. 32 23. 54 
26. 50 26. 37 26. 32 6. 42 26. 60 
29. 63 20. 50 29. 45 20, 55 29. 7] 
$2. 82 32. 70 32. 65 32. 75 32. 89 
36. 11 35. OS 35. 94 36. 05 36. 18 
39. 52 39. 40 39. 36 39. 46 39. 58 
3. O05 12. 92 12. 89 $2. 95 3. 11 
16. 67 16. 51 46. 48 16. 47 16. 73 
dO. 34 50. OS 14. OS 19. 99 V0. SS 
53. 87 3. 44 53. 40 53. 52 53. S6 
7. 16 i. 74 56. 73 57. O04 57. 17 
oO.6 
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Figure 3. Deviation SAMA 
R P1080 constantan from constantan sam ple A-] 


of constantan samples and 


sent the data for iron versus platinum and con- 
stantan versus platinum given in [3]. The 1913 
table was for iron versus constantan only and, conse- 
quently, the corresponding data for the individual 
elements versus platinum are not available for 
inclusion in figures 2 and 3. The constantan marked 
“F-1, F-6” was a single sample intended for use in 
combination with either the F—1 or F—6 iron sample. 
In figure 4 similar plots for the completed thermo- 
couples are given. In preparing these figures, the 
1913 and RP1080 tables were adjusted to the emf 
units and temperature scale used in this investiga- 
tion and currently in use at this Bureau. These are 
the absolute system of electrical units [7] and the 
International Temperature Scale of 1948 [8]. One 
absolute volt is equal to 0.99967 international volt 


TABLE 4, 

Sample Cc Mn P Ss 
A-l 0. 07 0. 25 0. 015 0. 031 
C-1 ‘na . 38 . 005 . O31 
D-1 . 04 . 29 . 009 . 024 
F-1 . 02 . 03 . 004 . 029 
F-3 . OS .43 . O80 . 028 
F-4 . 02 . 23 . O11 . 016 
F-5 . 06 a . 006 . 020 
F-6 . 05 . 22 . 006 . 018 








*Spectographic determination, 
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couples, the 
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mo- 


(U.S... The differences between the International 
Temperature Scales of 1948 and 1927 are not more 
than 0.8 deg F in the range of this investigation [9] 
and so are nearly negligible. 

The differences in shape of the various curves in 
figure 4 are mainly due to the iron elements because 
the deviations of the constantan elements from each 
other are approximately linear. It is well known 
that thermal emf depends both on chemical composi- 
tion and on the degree of strain. In hope of gaining 
some insight into the differences among the various 
thermocouples, the iron samples were analyzed chem- 
ically, and their hardnesses were determined. This 
work was done by the Chemistry and Metallurgy 
Divisions of the Bureau. The results are given in 
tables 4 and 5. 


Chemical and spectrographic analyses of iron samples 


[Values are given in percent 


Cu Si* Ni* Cr* Sn* 
0. 12 <0. 01 0. 044 0. O15 0. 01 
. 054 <. 01 . 030 . 016 . Ol 
.14 . 01 . 042 . 016 Ol 
. 15 . 01 . 039 . 025 . Ol 
. O18 . 02 . O11 01 Ol 
. 023 . 01 . 014 . Ol Ol 
. 022 . Ol . 018 . Ol . O01 
. 022 . 01 .012 Ol . Ol 


Tannin 5 Results of hardness measurements on samples of 


iron wire 


Equivalent 


Sample Vickers . Rockwell B 
number 
number 
Al 99.8 55 
C-1 106 59. 5 
D | 138 76 
| 127 71 
F-3 210 of 
F-4 83. 6 39 
F-5 18 76 
F t 11 73 


* Determinations were made with a 10-kg load applied through a square-based 


diamond-pyramid indenter 


3. Discussion 


The sharp changes in the slopes of the various 
curves in figure | at about 1,650° F are thought to 
be due to the 1913 table having been based in this 
region on an extrapolation that neglected the effect 
of the abrupt change in thermoelectric power of iron 
at the a—y transition. At lower temperatures the 
various curves show a characteristic undulation. It 
is not known to what extent this is due to errors in 
the original 1913 calibration and to what extent to 
differences between the chemical compositions char- 
acteristic of iron used for thermocouples in 1913 and 
at present 

It was hardly to be expected that the curves of 
figure 2 would be fully interpretable, inasmuch as 
the researen was designed to test these commercial 
materials in the “as received” condition, and did not 
provide for that control of the variables of chemical 
composition and heat treatment that would be neces- 
sary for a fundamental understanding of the results. 

The quantum theory of metals leads to the correct 
predictions that iron is thermoelectrically positive to 
platinum and that the effect of small percentages of 
copper on the thermoeleetric power is opposite in 
sign to the effect of manganese. However, a detailed 
explanation of the features shown in figure 2 is be- 
vond the present scope of theorv. Nevertheless, a 
few generalizations may be noted from the data. 

Roeser and Dahl [3] showed that the 1913 table 
could be approximated by using iron wires having 
several tenths percent of Mn, whereas the samples 
of which RP1080 was characteristic contained less 
than 0.1 percent of Mn. The effects of small per- 
centages of the minor elements that are common in 
commercial iron on the thermoelectric power of iron 
have been investigated by Finch [10], who described 
several compositions of iron that can be used in 
thermocouples. Finch showed that the thermo- 
electric power of iron at 500° C (932° F) is increased 
by additions of Cr, 5, or Mn and decreased by Sn, 
Si, P, Ni, and (in percentages greater than 0.10) 
Cu, whereas the effect of carbon is negligible. It 


234 





will be noted that the curve for iron sample F-1 
similar in shape to that for RP1O080 iron, and, indee« 
its manganese content of only 0.03 percent marks 
as an RP 1080-tvpe iron. In the other iron sample 
manganese is the principal minor element, and th 
thermal emf at 932° F (see fig. 2), with minor exce] 
tions, increases with the manganese concentratior 
as would be predicted from the work of Finch. 

The bunching of the curves in the gamma regio; 
above 1,650° F may be the result of differences iy 
electronic band-strue ture in alpha and gamma iron 
However, the stee p rise of F-1 and F—4 is exe eptiona! 
These samples were lowest in carbon (0.02 percent 
whereas the samples showing the opposite behavio: 


C-1 and F-3, were highest and second highest 
respectively, in carbon content (0.11 and 0.08 per 
cent). The effect thus may be connected with th: 


solution of carbon in the gamma-phase iron. A 
study of the phases present initially was not made 
and so no positive statement may be made on this 
point. 

Dahl has shown [11; 5, p. 1238] that thermocouple 
iron decreases in ihennes lectric power on heating at 
temperatures above 1,200° F. However, the anneal- 
ing of high-purity iron increases the thermoelectric 
power, so that it appears that the drift in the thermal! 
emf of the iron elements of thermocouples at high 
temperatures is primarily due to phase changes and 
chemical changes, and that relief of strain is not an 
important factor in this drift. 

The unique shape of curve F-3 may be related to 
any of the following factors: (1) Its manganese 
content was highest of all (0.43 percent); (2) Its 
phosphorus content was very much higher than that 
of the other samples; (3) Its hardness was unusually 
high. 

The curves in figure 3 show the approximately 
linear deviation that is characteristic of constantan [3] 

The maximum differences among the iron samples 
and also among the constantan samples amounted 
to about 0.8 mv at 1,400° F, which is equivalent to 
about 20 deg F in the indication of an iron-constantan 
thermocouple. In addition, there were marked 
differences in the curve shape for the iron samples 
Nevertheless, the calibrations of the completed 
thermocouples fell within a satisfactorily narrow band 
up to 1,500° F. The width of this band on either 
side of the relation for thermocouple A-1 was approxi- 
mately equivalent to commercial tolerances for the 
iron-constantan thermocouple, thus indicating the 
practicality of the A-1 relation for a commercial 
standard. However, the large differences among the 
wires of a given kind illustrate the impracticality 
of establishing standards for the emf of the individual 
elements (against platinum) as has been done for the 
elements of all the other common thermocouples. 

The deviations of the thermocouples below 32° F 
are quite large in terms of temperature, which indi- 
cates the necessity for care in establishing the 
calibrations of thermocouples in this range. How- 
ever, the deviations in emf are nearly linear, so that 
relatively few calibration temperatures need be used 
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4. Reference Table for 1913 type Iron-Con- 
stantan Thermocouples 


As a result of this investigation, the Thermocouple 
Calibration Curves Committee, Recorder-Controller 
Section, SAMA, has recommended adoption of the 
temperature electromotive force relationship of 
thermocouple A-1 as a tentative standard. Because 
of the known instability of iron-constantan thermo- 
couples above 1,600° F, and also because of the large 
increase in deviation of the thermocouples from each 
other in this region, the upper limit of the standard 
has been set at 1,600° F. (Manufacturers of this 
thermocouple do not recommend its use above 1,400° 
to 1,500° F, except in special circumstances.) The 
differences between the A-—1 relationship and_ the 
1913 table, when both are expressed in consistent 
electrical units and temperature seale, are given in 
table 6 


Differences hetween A-1 and 19183 tables for tron- 
constantan thermocouples 


TarRLe 6 


Temper- Emf A-1 table minus emf 
ature 1913 table 


Equivaic nt 


KF Villivolts de grees F 
32 0. 00 0.0 
100 02 7 
200 0! 3 
300 Ol 3 
100 01 3 
500 +. O2 7 
600 01 3 
700 02 7 
800 00 0 
900 02 7 
1, 000 03 9 
1, 100 03 4 
1, 200 04 1.2 
1, 300 . 08 0.0 
1, 400 Ol +.3 
1, 500 06 + 1.7 
1, 600 +. O68 1.7 


In December 1950, preliminary results of this in- | 


vestigation were summarized in a published article 
{12}. In this article it was pointed out that the 
temperature-emf relationship of sample A-1 was 
favored for the proposed new SAMA standard. In 
order to determine whether or not the A-1 relation- 
ship would be acceptable as an industry-wide stand- 


ard, replacing not merely the 1913 table but also the | 
RP1080 table, answers to the following questions | 


were solicited: (1) Is industry-wide standardization 
on a single reference table for iron-constantan desir- 
able? (2) Is the 1913 table, or a close approximation 
to it, acceptable for this purpose? (3) What degree 
of hardship will be occasioned by abandonment of 
the RP1080 table? 

The replies received were uniformly favorable to a 
single industry-wide standard. Users and suppliers 
of 1913-type thermocouples were uniformly favorable 
to the proposed SAMA standard. However, most of 
the replies received from users and suppliers of solely 


| 
| 
| 
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RP1080-type thermocouples indicated that it would 
not be practicable for them to change to the proposed 
SAMA standard. It is important to point out that 
the change from RP1080-type to A-l-type thermo- 
couples would require alterations to most of the 
associated indicating equipment, whereas the differ- 
ence between the 1913 and A-1 relationships is much 
smaller and is less than commercial tolerances for 
this thermocouple. 

Thus it appears that at least two standards, the 
SAMA table (A-1) and the RPIOSO table, will 
coexist in the immediate future. Most commercial 
production of iron-constantan will probably be based 
on the SAMA table, as it has been based on the 1913 
table in the past. On the other hand, one of the 
principal fields in which the RP1O80 standard is 
dominant—namely, in piston-engined military air- 
craft——is declining in importance. 

In order to make the A~1 relationship useful, it 
has been expanded into four tables, numbers 7, 8, 9, 
and 10. Before doing this, additional samples of the 
A lot were calibrated as a check against gross errors. 
The differences between these and the calibration of 
A-1 were within the limit of error of the calibration 
process. In table 7, the desired function for the 
region from —320° to 0° F was obtained at intervals 
of 1 deg by the method of divided differences. For 
the range from 0° to + 1600° F, the LaGrange poly- 
nomial was used. Linear interpolation, either direct 
or inverse, was then performed on this table to obtain 
the remaining tables. The preparation of tables 7 to 
10 was carried out by the Applied Mathematics 
Division of the Bureau. 


The thermoelectric measurements below 32° F 
were made in the Thermometry Laboratory by Ruth 
C. Gallman and Eugenia Thomas. It is a pleasure to 
acknowledge the advice given by Wm. F. Roeser and 
also the cooperation of the various members of the 
Thermocouple Calibration Curves Committee of 
SAMA and their chairman, Edward J. Albert. 
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Table 7. Iron-Constantan Thermocouples—Continued 
(Emf in Absolute Millivolts. Temperatures in Degrees F.” Reference Junctions at 32° F.) 



































°F 0 1 2 3 oe a ae oe 7 8 9 10 J °F 
; Millivolts 
(4)0 | -.89 -.86 -.83 -.80 -.78 |-.75 [-.72 -.70 -.67 -.64 -.61] (+)0 
10 | -.61 -.58 -.56 -.53 -.50 |-.48 |-.45 -.42 -.39 -.36 -.34 10 
20 | -.34 -.31 -.28 -.25 -.22 |-.20 |-.17 -.14 -.11 -.09 -.06 20 
30 | -.06 -.03 -.00 .03 -05 .0o8 | .11 14 17 19 3.22 30 
ko .22 .25 .28 31 34 .36 .39 -42 45 -48 .50 40 
50 50 .53 -56 59 #.62 65 .67 -70 -73 -76 = .79 50 
60 79 ~2-«.82 .84 .87 .90 | .93 .9%6 .99 1.02 1.04 1.07 60 
70 | 1.07 1.10 1.13 1.16 1.19 |1.22 2.25 1.28 1.30 1.33 1.36 70 
80 | 1.36 1.39 1.42 1.45 1.48 11.51 [1.54 1.56 1.59 1.62 - 80 
90 | 1.65 1.68 1.71 1.74 1.77 |1.80 |1.83 1.85 1.88 1.91 1. 90 
- 100 | 1.94% 1.97 2.00 2.03 2.06 |2.09 /2.12 2.14% 2.17 2.20 2.23] 100 
- 110 | 2.23 2.26 2.29 2.32 2-39 2.38 |2.41 2.44 2.47 2.50 2.52] 110 
120 | 2.52 2.55 2.58 2.61 2.64 |2.67 |2.70 2.73 2.76 2.79 2.82] 120 
130 | 2.82 2.85 2.88 2.91 2.9% {2.97 13.00 3.03 3.06 3.08 3.11] 130 
140 | 3.11 3.14 3.17 3.20 3.23 |3.26 13.29 3.32 3.35 3.38 3.41); 140 
150 | 3.41 3.44 3.47 3.50 3.53 13.56 13.59 3.62 3.65 3.68 3.71] 150 
160 | 3.71 3.74 3.77 3.80 3.83 |3.86 |3.89 3.92 3.95 3.98 4.01] 160 
170 | 4.01 4.04% 4.07 4.10 4.13 [4.16 /}4.19 4.22 4.25 4.28 4.31) 170 
180 4.32 4.34 4.37 4.40 4.43 14.46) 4.49 4.52 4.55 4.58 4.61 180 
190 4.61 4.64% 4.67 4.70 4.73 |4.76/4.79 4.82 4.85 4.88 4.91 190 
ie 200 | 4.91 4.94 4.97 5.00 5.03 |5.06/5.09 5.12 5.15 5.18 5.21] 200 
- 210 | 5.21 5.24% 5.27 5.30 5.33 |5.36|5.39 5.42 5.45 5.48 5.51] 210 
220 | 5.51 5.54 2-B) 5.60 5.63 |5.6615.69 5.72 5.75 5.78 5.81] 220 
230 | 5.81 5.84 5. 5.90 5.93 |5.9615.99 6.02 6.05 6.08 6.11] 230 
2ho | 6.11 6.14% 6.17 6.20 6.24 |6.27|6.30 6.33 6.36 6.39 6.42] 240 
) 250 | 6.42 6.45 6.48 6.51 6.54 |6.57|6.60 6.63 6.66 6.69 6.72] 250 
: 260 | 6.72 6.75 6.78 6.81 6.84 |6.87|6.90 6.93 6.96 7.00 7.03] 260 
270 | 7.03 7.06 7.09 7.12 7.15 |7.18|7.21 7.2% 7.27 7.30 7.33 afro 
280 | 7.33 7.36 7.39 7.42 7.45 |7.48|7.51 7.54 7.58 7.61 7.64] 280 
290 | 7.64 7.67 7.70 7.73 7.76 |17.79|7.82 7.85 7.88 7.91 7.94} 290 
300 | 7.94% 7.97 8.00 8.04 8.07 |8.10/8.13 8.16 8.19 8.22 8.25] 300 














jel, bc. Sees 


*Based on the International Temperature Scale of 1948. 






































Table 7, Iron-Constantan Thermocouples—Continued 
(Emf in Absolute Millivolts. Temperatures in Degrees F.” Reference Junctions at 32° F.) 
oF (9) I 2 3 soe oe ‘) TO oF 
Millivolts 

300 | 7.9% 7.97 8.00 8.04 8.07] 8.10] 8.13 8.16 8.19 8.22 8.25] 300 
310 | 8.25 8.28 8.31 8.34 8.20 8.40| 8.44% 8.47 8.50 8.53 8.56) 310 
320 BB 8.59 8.62 8.65 8.68] 8.71| 8.74% 8.77 8.80 8.84 8.87] 320 
330 | 8. 8.90 8.93 8.96 8.99] 9.02| 9.05 9.08 9.11 9.14 9.17 | 330 
340 9-27 9.20 9.24 9. 9.30| 9.33] 9.36 9.39 9.42 9.45 9.48] 340 
350 9.4 2-ae 9-38 2-25 9.61| 9.6 ae! 9.70 9.7 9.76 9.79 350 
360 | 9.79 9. 9.85 9. 9.92] 9.95] 9.98 10.01 10.04 10.07 10.10 | 360 
370 10.10 20-23 10.16 10.19 10.22 |10.25 /10.28 10.32 10.35 10.38 10.41 370 
380 |10.41 10.4 ses! 10.50 10.53 |10.56 /10.60 10.63 10.66 10.69 10.72 | 380 
390 |10.72 10.75 10.78 10.81 10.84 |10.87 |10.90 10.94 10.97 11.00 11.03 | 390 
4OO {11.03 11.06 11.09 11.12 11.15 /11.18 [11.21 11.24 11.28 11.31 11.34 | 400 
410 [11.34 we 11.40 11.43 11.46 [11.49 [11.52 11.55 11.58 11.62 11.65 | 410 
420 |11.65 11.68 11.71 11.74 11.77 [11.80 a1 83 11.86 11.89 11.92 11.96 | 420 
430 /11.96 11.99 12.02 12.05 12.08 /12.11 j12.14 12.17 12.20 12.23 12.26 | 430 
44O 112.26 12.30 12 2 12.36 12.39 |12.42 [12.45 12.48 12.51 12.54 12.37 4hO 
450 12.37 12.60 12. ie! 12.70 cet 12.76 12.79 12.82 12.85 12. 450 
460 |12.88 12.91 12.9% 12.98 13.01 |13.04 |13.07 13.10 13.13 13.16 13.19 | 460 
are 13.19 13.22 13.25 13.28 13.31 [13.34 (13.38 13.41 13.44 13.47 13.50 MA 
4 ap a0 13 «33 ¥ 13.59 13.62 |13.65 [13.68 13.72 13.75 13.78 13.81 | 4 

490 [13.81 13.84 13.87 13.90 13.93 [13.96 [13.99 14.02 14.05 14.08 14.12 | 490 
500 [14.12 14.15 14.18 14.21 14.24 [14.27 [14.30 14.33 14.36 14.39 14.42 | 500 
510 [14.42 14.45 14.48 ie 14.55 BP 38 14.61 14.64 ty | 14.70 14.73 | 510 
520 [14.73 14.76 14.79 14.82 14.85 (14.88 114.91 14.94 14.98 15.01 15.04 | 520 
530 |15.04 15.07 15.10 15.13 15.16 /15.19 [15.22 15.25 15.28 15.31 15.34 | 530 
540 |15.34 15-20 15.40 15.44% 15.47 [15.50 15-35 ¥ 15.59 15.62 15.65 | 540 
550 |15.65 15.68 15.71 15.74 15.77 (15.80 25.84 15.87 15.90 15.93 15.96 | 550 
560 |15.96 15.99 16.02 16.05 16.08 /16.11 116.14 16.17 16.20 16.23 16.26 | 560 
370 16.26 16.30 16.33 16.36 16.39 |16.42 [16.45 16.48 16.51 16.54 16.57 270 
580 16-37 16.60 16.63 16.66 16.69 [16.72 [16:75 16.78 16.82 16.85 16.88 | 580 
590 {16.88 16.91 16.94 16.97 17.00 |17.03 [17.06 17.09 17.12 17.15 17.18 | 590 
600 |17.18 17.21 17.24 17.28 17.31 [17.34 117.37 17.40 17.43 17.46 17.49 | 600 











* Based on the International Temperature Scale of 1948. 











Table 7. 
(Emf in Absolute Millivolts. 


Iron-Constantan Thermocouples—Continued 





Temperatures in Degrees F." Reference Junctions at 32° F.) 






































°F ©) 1 2 3 a ae 2 7 2 + 10 oP 
Millivolts 
600 |17.18 17.21 17.24 17.28 17.31 |17.34 117.37 17.40 17.43 17.46 17.49 | 600 
610 117.49 17.52 17.55 17.58 17.61 |17.64 |17.68 17.71 17.74 ye 17.80 | 610 
620 117.80 17.83 17.86 17.89 17.92 /17.95 17.98 18.01 18.04 18.08 18.11 | 620 
630 |18.11 18.14 18.17 18.20 18.23 /18.26/18.29 18.32 18.35 18.38 18.41 | 630 
640 |18.41 18.44 18.47 18.50 18.54 [18.57 |18.60 18.63 18.66 18.69 18.72 | 640 
650 |18.72 18.75 18.78 18.81 18.84 |18.87 |18.90 18.94 18.97 19.00 19.03 | 650 
660 119.03 19.06 19.09 19.12 19.15 |19.18|/19.21 19.24 19.27 19.30 19.34 | 660 
670 119.34 19.37 19.40 19.43 19.46 |19.49/19.52 19.55 19.58 19.61 19.64 | 670 
680 |19.64 19.67 19.70 19.74 19.77 |19.80/19.83 19.86 19.89 19.92 19.95 | 680 
690 119.95 19.98 20.01 20.04 20.07 [20.10 /20.13 20.16 20.20 20.23 20.26 | 690 
700 |20.26 20.29 20.32 20.35 20.38 /20.41 |20.44 20.47 20.50 20.53 20.56 | 700 
710 120.56 20.59 20.62 20.66 20.69 /20.72 |20.75 20.78 20.81 20.84 20.87 | 710 
720 120.87 20.90 20.93 20.96 20.99 [21.02 /21.05 21.08 21.11 21.14 21.18 | 720 
730 [21.18 21.21 21.24 21.27 21.30 [21.33 }21.36 21.39 21.42 21.45 21.48 | 730 
THO [21.48 21.51 21.54 Sz 21.60 [21.64 |21.67 21.70 21.73 21.76 21.79 | 740 
750 |21.79 21.82 21.85 21.88 21.91 |21.94 |21.97 22.00 22.03 22.06 22.10 | 750 
760 |22.10 22.13 22.16 22.19 22.22 |22.25 |22.28 22.31 22.34 22.37 22.40 | 760 
770 [22.40 22.43 22.46 22.49 22.52 |22.55 |22.58 22.62 22.65 22.68 22.71 | 770 
780 |22.71 22.74 22.77 22.80 22.83 |22.86 |22.89 22.92 22.95 22.98 23.01 | 780 
790 |23.01 23.04 23.08 23.11 23.14 |23.17 |23.20 23.23 23.26 23.29 23.32 | 790 
800 123.32 23.35 23.38 23.41 23.44 |23.47 123.50 23.53 23.56 23.60 23.63 | 800 
810 |23.63 23.66 23.69 23.72 23.75 [23.78 |23.81 23.84 23.87 23.90 23.93 | 810 
820 123.93 23.96 23.99 24.02 24.06 |24.09 |24.12 24.15 24.18 24.21 24.24 | 820 
830 |24.24 24.27 24.30 24.33 24.36 |24.39 |24 4 24.45 24.48 24.52 24.55 | 830 
840 |24.55 24.58 24.61 24.64 24.67 24.70 |24.73 24.76 24.79 24.82 24.85 | 840 
850 [24.85 24.88 24.91 24.94% 24.98 |25.01 |25.04 25.07 25.10 25.13 25.16 | 850 
860 |25.16 25.19 25.22 25.25 25.28 25.32 |25.35 25.38 25.41 25.44 25.47 | 860 
870 125.47 25.50 25.53 25.56 25.59 |25.62 |25.65 25.68 25.72 25.75 25.78 | 870 
880 (25.78 25.81 25.84 25.87 25.90 |25.93 |25.96 25.99 26.02 26.06 26.09 | 880 
890 |26.09 26.12 26.15 26.16 26.21 |26.24 |26.27 26.30 26.33 26.36 26.40 | 890 
900 (26.40 26.43 26.46 26.49 26.52 |26.55 |26.58 26.61 26.64 26.67 26.70 | 900 














* Based on the International Temperature Scale of 1948. 





Table 7. 
(Emf in Absolute Millivolts Temperatures in Degrees F.” Reference Junctions at 32° F.) 


Iron-Constantan Thermocouples—Continued 








°F 


0 


1 


2 


3 


4 


6 


7 


8 


9 


10 





Mil 


= 
llivolts 


°F 





900 


26. 


40 


26.43 


26 


46 


26 


-49 


26 


52 


26. 


55 


26 .58 


26. 


61 26. 


26 


67 


26.70 


900 





910 
920 
930 


940 
950 
960 
970 


980 
990 


26 
27 
27 


27 
2 
2 


28 


«70 
-02 
-33 


64 
38 
.58 


28. 
29. 


89 
21 


26.74 
27 .05 
27 .36 


26 
27 
27 


27 


28. 


28 


ge 28. 
29. 


-TT7 
.08 


39 


-70 
.02 
33 


64 
96 
27 


26. 
-ll 


27 
27 


27 


28 


80 
-42 


-73 
28. 
28. 


05 
36 


-67 
28. 
29. 


99 
30 


26. 


27 
27 


27 


28. 
28. 


28 


29. 
29. 


83 
14 


45 


-76 
08 
39 


-70 
02 


33 


26.86 
ot 
27 ie 
-80 
28. 
28. 


-74 
29. 
29. 


27 


27 


28 


11 
42 


05 
37 


26.89 
27 .20 
27.51 


27 .83 
28.14 
28 45 


28.77 
29.08 
29.40 


26 
27 
27 


eT. 
28. 
28. 


29 


92 26.95 
23 27. 
-54 27. 


86 27. 
-17 28. 
48 28. 


80 28. 
29,11 29. 
43 29. 


26. 
27. 
27. 


27. 
28. 
28. 


29. 
29. 


98 


27 .02 
27 .33 
27 .64 


27-23 
28.26 
28.58 
28.89 


29.21 
29.52 


910 
930 
940 
950 
960 
970 


990 





1000 


29. 


52 


29. 


59 


29. 


62 


29. 


65 


29. 


68 


29.71 


29. 


29. 


29.84 


1000 





1010 
1020 
1030 


1040 
1050 
1060 


1070 
1050 
1090 


31 


29. 
30. 
30. 


30. 
31. 
31. 


-76 
32. 
32. 


84 
16 
48 


80 
12 
yy 


08 
40 


30.19 
30.51 


30.83 
31.15 
31.47 


31.79 
32.11 
32.43 


29. 
30. 


30 


30. 
31. 
31. 


31 


32. 


32 


90 
22 


54 


86 
18 
50 


.82 
14 
-47 


29. 
30. 
30. 


30. 
31. 
31. 


31. 
32. 
-50 


32 


94 
25 
57 


89 
21 
53 


85 
18 


29 


30. 


30 
30 


31. 
31. 


31. 
32. 


32 


"26 


.60 


-92 


24 
56 


88 
21 
-53 


30. 
32 


30 


30. 


00 
64 


30.96 


he 
31. 


31. 
32. 
32 . 


30.03 
30.35 
30.67 


30.99 
31.31 
31.63 


31.95 
32.27 
32.60 


30. 
30. 
31. 
31.66 
31. 
32. 


30. 
30. 
31. 
31. 
32. 
66 32. 


30.16 
30.48 
30.80 


31.12 
31.44 
31.76 


32.08 
32.40 
32.72 


1010 
1020 
1030 


1040 
1050 
1060 


1070 
1050 
1090 





1100 


32 


-72 


32.76 


32 


-79 


32 


.82 


32 


86 


32. 


32.92 


32. 


33. 


33.05 


1100 





1110 
1120 
1130 


1140 
1150 
1160 


1170 
11 
1190 


33 
be) 
33 


34 
34 


35 
D5 
5 


05 
>T 
70 


.03 
36 
34.68 


.O1 
re: 


33.08 
33.41 
33.73 


34.06 


34 39 
34.72 


25-08 
35.3 
35.71 


33 


33. 


33 
34 


34. 


34 
35 
35 
35 


-1ll 
44 
-76 


.09 
4O 


-75 


.08 
-41 
-74 


33 
33 
33 


34 
34 
35 
35 
35 


.15 
-47 
.80 


tg 
"78 
esl 


45 
-78 


33. 
-50 
.83 


33 
33 


34. 
-49 
.82 


el 


34 
34 


35 


35. 
.81 


35 


18 


16 


Mm 


33. 
33. 
33. 


34. 
34. 
34. 


35. 
5 


35.84 


33.24 
33.57 
33.90 


gee 
34. 

3488 
35.21 


35.54 
35. 


33. 
33. 
34, 
34. 
35. 
35. 


33. 
33. 
3h. 
34. 


35.3 


35 


1 
-64 
98 


33.37 
33.70 
34.03 


34.36 
34.68 
35.01 


35.3 
35. 
36.01 


1110 
1120 
1130 


1140 
1150 
1160 


1170 
11 
1190 





1200 





36. 


01 36.05 36 


.08 


36. 


11 


36. 


15 





36. 


18 


36.21 





36. 


36 


dl 


36.35 





1200 





*Based on the International Temperature Scale of 1918. 

















(Emf in Absolute Millivolts. Temperatures in Derrees F.* Reference Junctions at 32° F.) 


Table 7. 





Iron-Constantan The:mocouples—Contanued 





°F 


fe) 


1 


fi 


3 





4 


Ls5] 8 


Le 
f 


8 


9 





10 | 








Millivolts 








1200 


1210 
1220 
1230 


1240 
1250 
1260 


1270 
1280 
1290 


36.01 


36 
36 
37 


35 
.69 
.02 


37 .36 
37-71 
38.05 


38.39 
38.74 
39.08 


36.05 


36 .38 
36.72 
37 .06 


37 - 
37 - 
38. 


38. 
38. 
39. 


40 
74 
08 


43 
77 
12 


36.08 


36.42 
36.75 
37.09 


37 .43 
37.78 
38.12 


38.46 
38.81 
39.15 


36.11 


36.45 
36.79 
Fl oa 


37 47 
37.81 
38.15 


38.50 
38.84 
39.19 


36. 


2 
4 


36 


37 
37 


38 
38 
39 


15 


-48 
.82 
37. 


50 
84 
38. 


16 


19 


~53 
.88 
22 


36.18 


52 
86 
37 .20 


37 .54 
37 .88 
38.22 


38. 
38. 
39. 


36. 
36. 
36. 
37. 
37 - 
37 « 
38. 


36 
36 


38. 
39. 


36.25 


36. 
36.92 
37. 


37. 
27 « 
38. 


38. 
38.¢ 
39.33 


36 


36. 
36. 
37. 
37. 


98 
32 


3 
3 


38. 


9 


Je 


39. 


.28 


62 
96 
30 


64 


67 
02 
36 


36.31 


36 .65 
36.99 
37 .33 


37 .67 
38.02 
38.36 


38.70 
39.05 
39.40 


36. 


36. 
27. 
37 . 


-71 
38. 
38. 


JT 


39. 


39. 


69 
02 
3€ 


05 
39 

















1300 
1310 
1320 
1330 


1340 
1350 
1360 


1370 
1380 
1390 


Ee 


1400 


39.43 


39.47 


39.50 


39.54 


39. 


57 


39. 





39.78 
40.13 
40.48 


83 
19 
54 
-99 


61 





40 
41 
41 
41 
42 


ho. 


82 
17 
52 


87 
22 


58 


39. 
ho, 
40, 


40, 
41, 
41, 


41.93 
42.29 
42.64 


43.00 


39.85 
40.20 
40.55 


40.90 
41.26 
41.61 


41.97 
42.32 
42.68 


39.89 
40.24 
40.59 


40.94 
41.29 
41.65 


42.00 
42 .36 
42.71 


39. 
-27 
.62 


40, 
~33 
.68 


42, 


ko 
4o 


41 
41 


42 
4 


92 


98 


o4 


~39 
-75 


39. 
40, 
40, 


41, 
41, 
41, 


42 0° 
42 43 
4o, 


) ho. 


39. 


39 


-71 


39.75 


39. 








40, 


40.38 


0 ho 


ho, 


40 


3 40 


41, 
41. 
41, 


42 


42 


15 


06 


41 
-76 


12 
4 7 
82 


r 


-53 
89 


40.10 
40.45 
40.80 


41.15 
41.50 
41,86 


42,22 
42.57 
42.92 


ho, 


4 


he 
41, 
41, 


4e. 
42, 


y 


Ge 


40, 


40, 











1340 
135¢ 
1360 
| 1370 


130 
| 1390 





43.03 


43.07 


43 


-10 


43. 


43.25 


43.28 


43. 





1410 
1420 
1430 


1440 |4%.39 44, 


1450 
1460 


1470 
1480 
1490 


———_——+ 


1500 


\4% 32 
143.68 


44.03 


(44.75 
45.10 


45 46 
45.82 
46.18 


43 35 
43.71 
44.07 


ho 
44.78 
45.14 


45.50 
45.85 
46.21 


43.39 
43.75 
44.10 


4 46 
44 82 
45.18 


45.53 
45.89 
46.25 


43 43 
43.78 
4414 


$b -50 
44 85 
45.21 
45.57 


45.92 
46 .28 


43 
43 
yi, 


4y 
4 
45 


45. 
45, 
.32 


46 


-46 
.82 
.18 


-53 
89 
25 


60 
96 


46. 


46. 


43, 
43. 
4h. 
4y 
yy: 
45, 


45. 


53 
0G 
025 


-60 
-96 
32 


.68 
.03 
39 


43, 
43, 
4h. 


KY, 
45, 
45, 


46. 
46. 


43 64 
44 00 
4Y 35 


44.71 
45.07 


45.42 


45.78 
46.14 
46.50 


43.6 
4y , 
44, 


A. 


ys 


5. 
45. 
45. 
46, 


ee 


46.53 


9 | 1430 












| 1400 


| 1410 
| 1420 


| 1440 
| 1450 
1460 


| 1470 
1480 
1490 





+ — 


46 53 


46.57 


46 .60 


46.64 


46 


-67 





46, 74 


46. 


46 85 


46. 





1500 





1510 
1520 
1530 


1540 
1550 
1560 


1570 
1580 
1590 


46.89 
47 24 
47 .60 


46.92 
47 .28 
47. 


63 
47. 


99 
48 34 
48. 


69 
49.04 
49.39 


49.74 


46.96 
47.32 
47 .67 


48, 
48, 
48. 


49,08 
49.43 
49.77 


02 
38 
73 


47.00 
47.35 
47.70 


48.06 
48.41 
48.76 


49.11 
49.46 
49.81 


47 
47 


48 


49, 
49, 
49. 


-03 
47. 
-74 


48, 
48, 


39 


09 
45 


.80 


15 
50 
84 


.10 
46 


.81 


.16 
-52 
87 


22 
-56 
91 


47. 
47. 
47. 


48, 
48, 
48, 


49, 
49, 
49.8 


47. 
47, 
47. 


48. 
48. 


47.21 
47 .56 
47.92 


48.27 
48 62 
48.97 


49.32 
49.67 
50.01 


47.2 
7. 
47 .° 
48, 
48, 
49, 


49 
49 


36 
-70 
50. 


05 


| 1510 











1520 


1530 


1549 


1550 
| 1569 


i570 
1580 
15390 





1600 


50.05 

















1609 





*Based on the International Temperature Scale of 1948. 


241 







(Emf in Absolute 


Table 8. 


Iron-Constantan Thermocouples 
Millivolts. Temperatures in Degrees F". Reference Junctions at 32° F) 






































Millivolts| .00 .10 .20 .30 h0 1-50 [00.70 80 .90 1.00 |Millivolts 
Degrees F 
-7.00 -267 -273 -279 -285 -292 |-298 |-305 -313 -320 -7.00 
-6.00 -212 -217 -222 -228 -233 |-238 |-244 -249 -255 -261 -267 -6.00 
9:00 -165 -169 -174% -178 -183 |-188 |-192 -197 -202 -207 -212 -5.00 
-4.00 -121 -125 -130 -134 -138 |-142 |-147 -151 -156 -160 -165 -4 00 
-3.00 -80 -84 -88 -92 -96 |-100 |-105 -109 -113 -117 -121 -3.00 
-2.00 -42 -45 -49 -53 -57| -61| -65 -69 -72 -76 -8C -2.00 
-1.00 -y -8 -12 -15 -19| -23/] -26 -30 -34 -38 -He -1,00 
(-)0.00 32 29 #2 21 18; 14] 21o 7 3 Oo -4] (-)0.00 
(+)0.00 32 36 39 «644630~=646] 50] 53 57 +60 64 67] (+)0.00 
1.00 67 71 74% 78 =81/ 85| 88 92 9 99 102 1.00 
2.00 102 105 109 112 116/119] 123 126 129 133 136 2.00 
3.00 136 140 143 146 150/ 153/156 160 163 166 170 3.00 
4.00 170 173 177 180 183] 187/190 193 196 200 203 4.00 
5.00 203 206 210 213 216] 220] 223 226 230 233 236 5.00 
6.00 236 240 243 246 249/| 253| 256 259 263 266 269 6.00 
7.00 269 272 276 279 282) 286| 289 292 2 299 302 7.00 
8.00 302 305 308 312 315] 318] 321 325 328 331 334 8.00 
9.00 334 338 341 344 347 | 351/354 357 360 364 367 9.00 
10.00 367 370 373 376 380 | 383| 386 389 393 396 399 10.00 
11.00 99 402 406 4oo 412/415] 419 422 425 428 431 11.00 
12.00 31 435 438 «441 44h | 44 4B 454 457 461 464 12.00 
13.00 464 467 #470 474% 4771 480) 4 4 490 493 496 13.00 
isco | 380 352 Sug Soo Se |Se|Ste G2 See 38 Set] to-c 
: 5 32 535 53 2 . 
16.00 561 565 5 571 574 278 581 23h 5 591 594 16.00 
17.00 594 597 601 604% 607 |610/] 614% 617 620 623 627 17.00 
18.00 627 630 633 636 640 | 643 | 646 649 653 oe 659 18.00 
19.00 659 662 666 669 672 | 6751/1679 682 6 692 19.00 
20.00 692 695 698 701 705 |708|711 714 718 721 724 20.00 











*Based on the International Temperature Scale of 1948. 














Table 8. 


Iron-Constantan Thermocouples—Continued 
(Emf in Absolute Millivolts. Temperatures in Degrees F". Reference Junctions at 32° F) 





























Willivolts ee) PS fo} 220 30 40 | 2 DU i _ 260 2/0 250 * 50 1,00 Millivolts 
Degrees F 
20.00 692 695 698 701 705| 708] 711 714 718 721 724 20.00 
21.00 724% 727 731 73% 737) 7411 784% 787 #750 754 757 21.00 
22.00 757 760 763 767 q7o 4 77 780 q , £20 22.00 
23.00 790 793 796 799 803 09 812 816 819 822 23.00 
24.00 822 825 829 832 835) 839) 842 845 848 852 855 24.00 
25.00 Be 358 861 865 828 871| 874 878 881 Beh 887 25.00 
26.00 890 894% 897 900] 903} 907 910 913 916 919 26.00 
27.00 919 923 926 929 932] 936] 939 942 945 948 2 27-00 
28.00 i ee tee Ses oetl Sablaeee nee ance ute 28.00 
29.00 990 993 996] 999/1002 1006 1009 1012 1015 29.00 
30.00 1015 1018 1021 1024 1028/1031/1034 1037 1040 1043 1046 30.00 
31.00 1046 1050 1055 1056 1059/1062/1065 1068 1071 1074 1078 31.00 
32.00 1078 1081 1084 1087 1090/1093/1096 1099 1102 1105 1108 32.00 
33.00 1108 1112 1115 1118 1121/1124/1127 1130 1133 1136 1139 33.00 
34.00 1139 1142 1145 1148 ey cM 1157 1160 1164 1167 1170 34.00 
35.00 1170 1173 1176 1179 1182/1185/1188 1191 1194 1197 1200 35.00 
36.00 1200 1203 1206 1209 1212/1215|1217 1220 1223 1226 1229 36.00 
37.00 1229 1232 1235 1238 1241/1244/1247 1250 1283 1256 12 37.00 
38.00 1259 1261 1264 1267 1270/1273|1276 1279 1282 1285 12 38.00 
39.00 1288 1290 1293 1296 1299/1302/1305 1308 1311 1313 1316 39.00 
40.00 1316 1319 1322 1325 1328/1331/1333 1336 1339 1342 1345 40.00 
41.00 1345 1348 1350 1353 13 1339 1362 1365 1367 1370 1373 41.00 
42.00 1373 1376 1379 a 1 13 13 1390 1393 1395 1398 1401 42.00 
43.00 1401 1408 1807 1410 1412/1415/1418 1421 1424 1426 1429 43.00 
44.00 1429 1432 1835 1438 1440/1443/1446 1449 a 1454 14 44.00 
45.00 1457 1460 1563 1465 1468/1471/1474 1477 1480 1482 1485 45.00 
46.00 1485 1488 1491 1494 1496/1499/1502 1505 1508 1510 1513 46.00 
47.00 1513 1516 1519 1522 1524/1527/1530 1533 1536 1539 1541 47 .00 
48.00 1541 1544 1547 1550 1553 1596 1558 1561 1564 1567 1570 48.00 
49.00 1570 1573 1576 1578 1581 1584/1587 1590 1593 1596 1599 49.00 
50.00 1599 1602 50.00 

















*Based on the International Temperature Scale of 1948. 
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Table 9. Iron-Constantan Thermocouples—Continued 
(Emf in Absolute Millivolts. Temperatures in Degrees C(Int.1948) Reference Junctions at 0° C) 



































eC 0 I 2 3 a 7 8 7 10 °C 
Millivolt 

100 | 5.27 5.32 5.38 5.43 5.48] 5.5415.59 5.65 5.70 5.76 5.81 | 100 
110 | 5.81 5.86 5.92 5.97 6.03] 6.08] 6.14 6.19 6.25 6.30 6.36 | 110 
120 | 6.36 6.41 6.47 6.52 6.58] 6.63] 6.68 6.74% 6.79 6.85 6.90 | 120 
130 | 6.90 6.96 7.01 7.07 7.12) 7.18] 7.23 7.29 7.34 7.40 7.45 | 130 
140 | 7.45 7.51 7.56 7.62 7.67] 7.73| 7.78 7.84% 7.89 7.95 8.00 | 140 
150 | 8.00 8.06 8.12 8.17 8.23] 8.28) 8.34 8.39 8.45 8.50 8.56 | 150 
160 | 8.56 8.61 8.67 8.72 8.78 8.84 | 8.89 8.95 9.00 9.06 9.11 160 
170 | 9.11 9.17 9.22 9.28 9.33] 9.39| 9.44 9.50 9.56 9.61 9.67 | 170 
180 | 9.67 9.72 9.78 9.83 9.89] 9.95 [10.00 10.06 10.11 10.17 10.22 | 180 
190 |10.22 10.28 10.34 10.39 10.45 /10.50 110.56 10.61 10.67 10.72 10.78 | 190 
200 {10.78 10.84 10.89 10.95 11.00 /11.06 [11.12 11.17 11.23 11.28 11.34 | 200 
210 (11.34 11.39 11.45 11.50 11.56 |11.62 [11.67 11.73 11.78 11.84 11.89 | 210 
220 {11.89 11.95 12.00 12.06 12.12 /12.17 |12.23 12.28 12.34 12.39 12.45 | 220 
230 |12.45 12.50 12.56 12.62 12.67 /12.73 (12.78 12.84 12.89 12.95 13.01 | 230 
240 |13.01 13.06 13.12 13.17 13.23 |13.28 (13.34 13.40 13.45 13.51 13.56 | 240 
250 |13.56 13.62 13.67 13.73 13.78 |13.84 113.89 13.95 14.00 14.06 14.12 | 250 
260 [14.12 14.17 14.23 14.28 14.34 |14.39 14.45 14.50 14.56 14.61 14.67 | 260 
270 (14.67 14.72 14.78 14.83 14.89 |14.94 [15.00 15.06 15.11 15.17 15.22 | 270 
280 [15.22 15.28 15.33 15.39 15.44 /15.50 [15.55 15.61 15.66 15.72 15.77 | 280 
290 |15.77 15.83 15.88 15.94 16.00 /16.05 (16.11 16.16 16.22 16.27 16.33 | 290 
300 /16.33 16.38 16.44% 16.49 16.55 |16.60 |16.66 16.71 16.77 16.82 16.88 | 300 
310 (16.88 16.93 16.99 17.04 17.10 |17.15 |17.21 17.26 17.32 17.37 17.43 | 310 
320 |17.43 17.48 17.54 17.60 17.65 |17.71 |17.76 17.82 17.87 17.93 17.98 | 320 
330 (17.98 18.04 18.09 18.15 18.20 |18.26 /18.32 18.37 18.43 18.48 18.54 | 330 
340 (18.54 18.59 18.65 18.70 18.76 |18.81 |18.87 18.92 18.98 19.03 19.09 | 340 
350 (19.09 19.14 19.20 19.26 19.31 |19.37 /19.42 19.48 19.53 19.59 19.64 | 350 
360 [19.64 19.70 19.75 19.81 19.86 |19.92 /19.97 20.03 20.08 20.14 20.20 | 360 
370 20.20 20.25 20.31 20.36 20.42 [20.47 |20.53 20.58 20.64 20.69 20.75 | 370 
380 20.75 20.80 20.86 20.91 20.97 [21.02 21.08 21.13 21.19 21.24 21.30 | 380 
390 21.30 21.35 21.41 21.46 21.52 |21.57 |21.63 21.68 21.74 21.79 21.85 | 390 
400 (21.85 21.90 21.96 22.02 22.07 [22.13 22.18 22.24% 22.29 22.35 22.40 | 400 
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Table 9. Iron-Constanten Thermocouples—Continued 
(Raf in Absolute Millivolts. Temperatures in Degrees C(Int.1948) Reference Junctions at 0° C) 
*c 0 1 ae 2 2 2 2 7 8 , 3 10; 
Millivolts ne 
400 22. 65 21.90 3 ® 96 22.02 22. 07 2.13 ]22.18 22.24 22.29 22. 35 22. 4o 400 
410 D2 ho 22. 46 22.51 22.57 22.62 22.68 [22.73 22.79 22.84 22.90 22, 95 | 410 
420 |22.95 23.01 23.06 23.12 23.17 (23.23 [23.28 23.34 23.39 23.45 23.50 | 420 
430 /23.50 23.56 23.61 23.67 23.72 [23.78 [23.83 23.89 23.94 24.00 24.06 | 430 
44O 124.06 24.11 24.17 24.22 24.28 (24.33 [24.39 246.44 24.50 24.55 24.61 | buo 
450 |24.61 24.66 24.72 24.77 24.83 [24.88 24.94 25.00 25.05 25.11 25.16 | 450 
460 [25.16 25.22 25.27 25.33 25.38 25.44 125.49 25.55 25.60 25.66 25.72 | 460 
470 (25.72 25.77 25.83 25.88 25.94 125.99 26.05 26.10 26.16 26.22 26.27 470 
480 126.27 26.33 26.38 26.44% 26.49 |26.55 126.61 26.66 26.72 26.77 26.83 | 480 
490 |26.83 26.89 26.94 27. 00 27. -05 P27. 1 27.17 27.22 27.28 27.33 27.39 | 490 
500 [27.39 27.45 27 .50 27. 56 27. 62 27. 67| 27.73 27.78 27.84% 27.90 27.95 | 500 
510 27-95 28.01 28.07 28. 12 28.18 28. 23 128.29 28.35 28.40 28.46 28.52 510 
520 |28.52 28.57 28.63 28.69 28.74 |28.80 [28.86 28.91 28.97 29.02 29.08 | 520 
530 |29.08 29.14 29.20 29.25 29.31 [29.37 129.42 29.48 29.54 29.59 29.65 | 530 
540 |29.65 29.71 29.76 29.82 29.88 129.9% 29.99 30.05 30.11 30.16 30.22 | 540 
550 |30.22 30.28 30.34 30.39 30.45 [30.51 [30.57 30.62 30.68 30.74 30.80 | 550 
560 /30.80 30.85 30.91 30.97 31.02 [31.08 [31.14 31.20 31.26 31.31 31.37 | 560 
570 |31.37 31.43 31.49 31.54 31.60 [31.66 [31.72 31.78 31.83 31.89 31.95 | 570 
580 /31.95 32.01 32.06 32.12 32.18 [32.24 [52.30 32.36 32.41 32.47 32.53 580 
590 /32.53 32.59 32.65 32.71 32.76 |32.82 [32.88 32.94 33.00 33. 06 33.11 590 
600 |33.11 33.17 33.23 33.29 33.35 33.41 33.46 33.52 33.58 33. 6h 33.70 | 600 
610 [33.70 33.76 33.82 33.88 33.94 [33.99 [34.05 34.12 34.17 34.23 34.29 | 610 
620 [34.29 34.35 34.41 34.47 34.53 [34.58 14.64 34.70 34.76 34.82 34.88 | 620 
630 (34.88 34.94 35.00 35.06 35.12 [35.18 [35.24 35.30 35.36 35.42 35.48 | 630 
640 135.48 35.54 35.60 35.66 35.72 [35.78 [35.84 35.90 35.96 36.02 36.08 | 640 
650 |36.08 36.14 36.20 36.26 36.32 [36.38 [36.44 36.50 36.56 36.62 36.69 | 650 
660 |36.69 36.75 36.81 36.87 36.93 [36.99 [37.05 37.11 37.18 37.24 37.30 | 660 
670 37.30 37.36 37.42 37.48 37.54 137.60 4 -66 37.73 37.79 37.85 37.91 | 670 
680 |37.91 37.97 38.04 38.10 38.16 [38.22 28 38.34 38.41 38.47 38.53 | 6380 
690 |38.53 38.59 38.66 38.72 38.78 [38.84 58. 90 38.97 39.03 39.09 39.15 | 690 
700 |39.15 39.22 39.26 39.34 39.40 [39.47 39.53 39.59 39.65 39.72 39.78 | 700 
710 [39.78 39.84 39.91 39.97 40.03 0.10 /40.16 40.22 40.28 40.35 40.41 710 
720 }40.41 40.48 40.54 40.60 40.66 40.73 |40.79 40.86 40.92 40.98 se 720 
730 }41.05 41.11 41.17 41.24 41.30 41.36 141.43 41.49 41.56 41.62 41.6 730 
740 141.68 2-7 41.81 41.87 41.94 2.00 |42.07 42.13 42.19 42.26 42.32 | 740 
750 [42.32 42.38 42.45 42.51 42.58 42.64 |42.70 42.77 42.83 42.90 42.96 | 750 
760 |42.96 43.02 43.09 43.15 43.22 43.28 [43.35 43.41 43.48 43.54 43.60 | 760 
770 |43.60 43.67 43.73 43.80 43.86 43.92 [43.99 44.05 44.12 44.18 44.25 | 770 
780 44.25 44.31 44.38 44.4% 44.50 4.57 [44.63 44.70 44.76 44.82 44.89 | 780 
790 4y. 89 44.95 45.02 45.08 45.15 5.21 145.28 45.34 45.40 45.47 45.53 | 790 
800 45. 53 45.60 45 .66 45.72 45.79 65.85 |45.92 45.98 46.05 46.11 46.18 | 800 
810 /}46.18 46.24 46.30 46.37 46.43 46.50 |46.56 46.62 46.69 46.75 46.82 | 810 
820 }46.82 46.88 46.94 47.01 47.07 47.14 [47.20 47.27 47.33 47.39 47.46 | 820 
830 [47.46 47.52 47.58 47.65 47.71 47.78 [47.84 47.90 47.97 48.03 48.09 | 830 
840 48.09 48.16 48.22 48.28 18-33 8.41 |48.48 48.54 48.60 48.66 48.73 | B40 
850 |48.73 48.79 48.85 48.92 48. 9.04 |49.10 49.17 49.23 49.29 49.36 | 850 
860 [49.36 49.42 49.48 49.54 49.61 49.67 |49.73 49.79 49.86 49.92 49.98 | 860 
870 ho.se 50.04 870 























Table 10, Iron-Constantan Thermocouples 
(Emf in Absolute Millivolts. Temperatures in Degrees C(Int.1948) Reference Junctions at 0° C) 





Millivolts 200 210 220 230 240 | 250 | 200 2/0 280 290 1.00 | Millivolts 












































Degrees C 
-7.00 -166 -169 -173 -176 -180 +184 |-188 -192 -196 -7 .00 
-6.00 -136 -138 -141 -144 -147 +150 |+-153 -156 -159 -163 -166 6.00 
4< -109 -112 -114 -117 | L122 |-125 -127 -130 -133 -136 -5.00 
-4 00 -85 -87 -90 -92 - -97 | -99 -102 -104 -107 -109 -4 00 
-3.00 -62 -65 -67 -69 -71]-741]-76 -78 -80 -83 -85 -3.00 
-2.00 -41 -43 -45 =-47) «-49] -52/| -54 -56 -58 -60 -62 -2.00 
-1.00 -20 -22 -24 -26 -28] -30/ -32 -35 -37 -39 -41 -1.00 
(-)0.00 0 -2 -& -6 -8/-10/-12 -14 -16 -18 -20 |(-)0.00 
(+)0.00 0 2 4 6 8} 10] 12 18% 16 18 20 |(+)0.00 
1.00 20 22 2h 25 e27] 29] 31 33 #35 +37 «+39 {| ~= 1.00 
2.00 39 41 43 45 47] 48} SO 52 5S4& 56 58 2.00 
3.00 58 60 62 63 65] 67! 69 71 #=73 #75 +#«77 3.00 
4.00 77 78 80 82 84 86 88 90 91 93 95 4.00 
5.00 95 97 99 101 102/104/106 108 110 112 113 5.00 
6.00 113 115 117 #%4119 121/)123/124% 126 128 130 132 6.00 
[-0 132 134 #135 137 #139/141/143 144 146 148 150 7.00 
.00 150 152 154 155 157/159/161 163 164 166 168 8.00 
9.00 168 170 172 173 1751177/179 181 182 184 186 9.00 
10.00 186 188 190 191 193|1951197 199 200 202 204 10.00 
11.00 204 206 208 209 211/213 /215 217 218 220 222 11.00 
12.00 222 224% 226 227 229 ana 233 235 236 238 240 12.00 
13.00 240 242 26h 245 247 | 249 | 251 253 254 256 258 13.00 
14.00 258 260 262 263 265] 267/269 271 272 274 276 14.00 
15.00 276 278 280 281 283 | 285 | 287 289 290 292 294 15.00 
16.00 294 296 298 300 301 | 303 | 305 307 309 310 312 16.00 
17.00 312 314% 316 318 319 | 321 [323 325 327 329 330 17.00 
18.00 330 332 334 336 338 | 339) 341 343 345 347 348 18.00 
19.00 348 350 352 354% 356 | 357/359 361 363 365 366 19.00 
20.00 366 368 370 372 374 | 3761377 #379 381 383 385 20.00 
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Table 10, 


Iron-Constantan Thermocouples—Cont inued 
(Emf in Absolute Millivolts. Temperatures in Degrees C(Int.1948) Reference Junctions at 0° C) 









































~Millivolts eo 01D 220 0630 040 | 650 | 60) =—670 S80 Siw 90 sd. 00 | Millivolts 
Degrees C 

20.00 366 368 370 372 374! 376/377 379 381 383 385 20.00 
21.00 385 386 388 390 392] 394| 395 397 399 401 403 21.00 
22.00 hos 405 406 408 410] 412/414 #416 417 +419 421 22.00 
23.00 42] 4423 «42k «426 64428] 430] 432 «434 «435 437 439 23.00 
24.00 439 441 445 444 446) 448) 450 452 453 455 457] 24.00 
25.00 457 459 461 463 464) 466) 468 470 472 473 «475 25.00 
26.00 475 477 479 +481 482) 484/| 486 488 489 491 493 26.00 
27.00 493 495 497 499 500] 502| 504 506 507 509 511 27 .00 
28.00 511 513 514 516 518] 520| 522 523 525 527 529 28.00 
29.00 529 530 532 534 536] 537/539 541 543 5Sh4 546 29.00 
30.00 546 548 550 551 553| 5551557 558 560 562 564 30.00 
31.00 564 565 567 569 571/ 572| 574 576 577 579 581 31.00 
32.00 581 583 584 586 588] 589/591 593 595 596 598 32.00 
33.00 598 600 601 603 605| 607/608 610 612 613 615 33.00 
34.00 615 617 618 620 622] 624| 625 627 629 630 632 34.00 
35.00 632 634 635 637 639| 640/642 644 645 647 649 35.00 
36.00 649 650 652 654 655] 6571659 660 662 664 665 36.00 
37.00 665 667 668 670 672) 6731675 677 678 680 681 37 .00 
38.00 681 683 685 686 688| 690/691 693 694% 696 698 38.00 
39.00 698 699 701 702 704| 706|707 709 710 712 713 39.00 
40.00 713 715 717 718 720] 7211723 725 726 728 729 40.00 
41.00 729 731 732 734% 736) 737|739 7T4O T42 743 745 41.00 
42.00 745 747 748 750 751| 753 |754 756 757 759 761 42.00 
43.00 761 762 764 765 767)! 7681770 772 773 775 776 43.00 
44.00 776 778 779 781 782) 784 |785 787 789 790 792 44 00 
45.00 92 793 795 796 798/ 799/801 803 804 806 807 45.00 
46.00 07 809 810 812 814| 815/817 818 820 821 823 46.00 
47.00 S23 824 826 828 829/ 8311832 834 835 837 839 47.00 
48.00 839 840 842 843 845) 846 | 848 850 851 853 854 48.00 
49.00 854 856 858 859 861) 862 | 864 866 867 869 870 49.00 
50.00 870 872 50.00 








Wasuinoton, March 11, 1953. 
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Epitaxial Deposits of Metals Evaporated on Salt Substrates 
Olive G. Engel 


Experimental data and theoretical considerations relating to the process of epitaxy are 


reviewed, 


The case of epitaxy resulting from vapor phase deposition of metals on rock salt 


is considered in detail and it is concluded that a factor other than best fit determines the 


onrentation 
tion of the metal atoms on the surface. 


the ionization potential of the metal for the ionization state to be expected. 


\ mode! for epitaxy in metal—rock salt pairs is advanced which involves ioniza- 
The epitaxial temperature is found to be related to 


It is found that 


the observed orientations of metals on rock salt are those that would be expected if such an 
ionization process occurred, that is, if the orientation involved charged ions rather than 


neutral atoms 


l. Introduction 


The increasing importance that is being attached 
to thin films has been pointed out by Holland [1],? 
and a number of applications have been cited by 
Lewis [2], who arbitrarily defines “thin” to mean of 
the order of 100 A. Recently the Electron Physics 
Section of the National 
needed films that are not only thin by this criterion 
but that also retain the perfect orientation of a single 
crystal. Films of this kind were required both for 
the electron interferometer of Marton [3] and for 
experiments on the single inelastic scattering of 
electrons. 

For the interferometer a thickness of 100 A is per- 
missible, but the single-crystal structure must be 
perfect enough so that the diffraction pattern will 
consist of sharply defined spots when an electron 
beam is passed through the film. For the scattering 
experiments, films of a single chemical species were 
desired, and to make negligible the possibility of 
plural and multiple inelastic scattering, a thickness 
of less than 100 A was indicated. It appeared that 
single-crystal metal films of this degree of thickness 
would be satisfactory for both applications, and that 
such films might be grown successfully on a substrate 
that would orient the growing deposit. Knowledge 
gained from an analysis of the experimental data on 
the formation of such films, along with new additions 
to the theory, are presented in this review of the 
process of epitaxy. 


2. Analysis of Existing Data 


The term epitaxy has been coined by Royer [4] 
to designate the oriented overgrowth of one crystal- 
line substance on another. Although Royer’s work 


was concerned entirely with growing a salt from its | 


aqueous solution onto a single crystal of a different 
salt, epitaxy may also result when the deposit ma- 
terial is condensed from the vapor phase onto the 
substrate, when the deposit material is electrodepos- 
ited onto the substrate, or when the deposit material 
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is produced as a product of a chemical reaction that 
takes place on the surface of the substrate. 

Results of past experimental work have indicated 
two general conditions that are necessary for the 
formation of oriented single-crystal metal films. 
These are (1) a certain amount of fit between the 
lattice of the substrate and that of the deposit ma- 
terial, and (2) a characteristic substrate temperature 
(epitaxial temperature) associated with every metal- 
substrate pair below which perfect orientation is 
not obtained. 

For use in the interferometer and in the scattering 
experiments, the film must be removed from the 
substrate. Rock salt, which is easily obtained in 
single-crystal form, and which can be removed from 
the metal films by solution in water, is a satisfactory 
substrate from this standpoint. By using substrate 
temperatures that are tabulated in the literature for 
about six metal-rock salt pairs, very thin films of a 
number of these metals have been condensed from 
the vapor phase by J. A. Suddeth and T. F. McCraw 
of the Electron Physics Section of the National 
Bureau of Standards. Not all of these films have 
desirable properties. The gold films are very nicely 
oriented if the rock salt cleavage face is maintained 


at 400° C during the deposition of the metal. They 
are transparent; blue light is transmitted. They 


have, however, the undesirable property that they 
collapse under their own weight at this extraordinary 
thickness because of the low interatomic forces of 
attraction between the gold atoms. Figure | is the 
electron diffraction pattern of one of the few gold 
films of this order of thickness that were successfully 
mounted by Suddeth and McCraw. 

The electron diffraction patterns of silver films 
show that these films are well oriented if the cleavage 
face of the rock salt is maintained at 150° C during 
the deposition of the metal. Silver films of this 
order of thickness are sturdy but are full of small 
holes. See, for example, the electron micrograph 
shown in figure 2, where the angular structure of the 
holes (angles of 45 and 90 degrees) clearly indicates 
the cubic lattice type for silver. Johnson [5] has 
followed the structural changes of thin oriented silver 
films that were condensed in a vacuum on rock-salt 
cleavage faces as a function of temperature. For 
films of 1000 to 2000 A he found that rectangular 





Fioure 1 Electron diffraction pattern of a thin gold film 
eva porated onto the cleavage face of rock salt. 


holes appeared in the films. Their edges were alined 
with the [110] direction of the rock salt He attrib- 
uted the effects to the readjustment of the oriented 
surface under the influence of surface-tension forces. 

To extend experimental investigation to a variety 
of metals for which data are not available in the 
literature, it became desirable to be able to predict 
reasonably well not only the best substrate for maxi- 
mum orientation but also the characteristic epitaxial 
temperature to be expected for the pair. Generali- 
zations comparable to the rules of Rover [4] for epi- 
taxial deposits from solution do not exist in the case of 
condensation of metals from their vapors. The fol- 
lowing approach to the problem is an analysis of 
existing data to find a basis on which such predic 
tions can be made 


2.1. Criterion of Best Fit 


In regard to the first condition, it appears that in 
the case of metal films deposited from the vapor in 
a vacuum too much importance has been given to 
the degree of lattice fit between the deposit and sub- 
strate. This overemphasis of lattice fit seems to 
have grown out of the extensive work of Royer [4], 
Sloat and Menzies [6], and others on the epitaxy of 
salts deposited on various substrates from solution 
under which condition the degree of fit was appar- 
ently rather critical and the dielectric properties of 
the solvent played a role. It was further entrenched 
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Ficure 2 Electron micrograph of a thin silver fiim evaporated 


onto the cleavage face of rock salt 


by the mathematical treatment of Frank and van 
der Merwe i, S| Recently, however, Schulz [9] 
has concluded, even for deposition from aqueous 
solution, that close atomic matching at the inter- 
face is not required for an oriented overgrowth. He 
obtained oriented films by deposition from the vapor 
phase. He then dissolved these films by condensing 
water vapor on them. Finally, he regenerated the 
oriented film by evaporation of the water. In eval- 
uating his conclusion one must bear in mind that in 
the method that Schulz used the ions from an ori- 
ented deposit were just barely brought into solution 
by condensation of water vapor onto the film with- 
out a mechanical action so that the ions of the first 
oriented monolayer on being brought into solution 
probably never strayed far from the site on the sub- 
strate that they had occupied. In this way the 
oriented nucleus in solution continued to. exist. 
Pashley [10] has oriented silver chloride on mica 
from the vapor phase with a misfit of —24 percent 
in spite of the fact that Royer [4] found the misfit 
limit of rock-salt lattice-type salts on mica from 
solution to be about 15 percent. On the basis of 
this evidence Pashley has suggested that the method 
of growth of oriented deposits is an important factor 
to be taken into account. 

While epitaxy that results from crystallizing the 
deposit material from solution may take place under 
fairly critical conditions of fit, results of experimental! 
work on the deposition of metal films from their 
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Figure 3. The | 100}, {110}, and {111} faces of 
a cubic crystal structure. 
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Figure 4. Silver on the {100} face of rock salt at 
150° C. 


@. Silver; X, sodium ior , chloride ion 


vapors have shown that there is a wide tolerance in 
fit in many cases in which epitaxy bas been found, 
and that the depositing vapor often chooses to con- 
dense in an orientation with respect to the substrate 
that from the standpoint of fit would be undesirable. 
To illustrate the point the orientation of silver depos- 
ited on the {100),/110}, and {111 ) faces of rock salt 
will be reviewed (see fig. 3). 

Briick [11] in depositing silver from the vapor in a 
vacuum onto the cube face of rock salt at the char- 
acteristic epitaxial temperature for perfect orienta- 
tion on that face found that the silver single crystal 


251 








© 


= 








FOUND NOT FOUND 
Figure 5. Silver on the | 100} face of rock salt at 450° C. 
@. Silver; X, sodium ion , chloride ion 
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Figure 6. Silver on the |110\ face of rock salt 


@. silver , chloride ion; X, sodium ion 


grows with its (100) face resting on the substrate 
and lying parallel to the cube edge of the rock salt 
with 27.5-percent misfit in spite of the fact that it 
would have a better fit (2.5-percent misfit) if it grew 
turned at 45 degrees to the sodium chloride cell 
(see fig. 4). Gharpurey [12] deposited silver on the 
cleavage (100) face of rock salt at an elevated tem- 
perature (450° C) and found that at this temperature 
the silver deposits with its 110. face resting on the 

100. face of the rock salt. It was inclined at 45 
degrees to the rock-salt cube with misfits of 2.5 
percent and 27.5 percent along the two directions of 
the rectangular parallelogram of the 110) face of the 
silver cube. From figure 5 one sees that an orienta- 
tion parallel to the rock salt cube would have given 
the same percent misfits. It is, however, not 
reported. 

Bré and Gharpurey [13, 14] found that when 
silver is deposited from its vapor onto a 110° face 
of rock salt, the (110) face of the silver rests on and 
is parallel to the (110) face of rock salt. The misfit 
is 27.5 percent in each direction considering the 
chlorine lattice. If the silver were rotated through 
90 degrees, the fit would be nearly perfect as is clear 
from figure 6, but this orientation is not reported. 
Similarly, on depositing silver from the vapor onto 
the 111) face of rock salt, they found that the 111 


face of silver was resting on the 111) face of sodium 
chloride (see fig. 7). The misfit on the (111. plane 
is 27.5 percent for the silver lattice. 

Briick [11], who investigated the orientation of 
nine different metals on the cleavage face of rock 
salt, observed that whereas the cubic face-centered 
metals grew with their (100. face resting on the 

100) face of the rock salt and oriented in the same 
way as the rock salt (i. e., 
parallel to the rock-salt cube), the cubic body- 
centered metals grew on the contact surface turned 
at 45 degrees to the rock-salt cube. In the case of 
the cubic body-centered metals, then, the orientation 
giving best fit is observed. 

From this evidence it appears that best fit is not 
the only criterion that determines whether or not 
epitaxy will occur in metals deposited from their 
vapors, or what orientation will be observed if it does 
occur. There appears to be some other factor 
involved that is more important than fit. An 
interpretation of the orientations that have been 
observed will be discussed at length in section 5, 
and the experimental observations reviewed here 
will be reexamined in the light of that discussion. 

Briick found experimentally (by trial and error) 
the epitaxial temperature to produce single-crystal 
films of Cu, Pd, Au, Ag, Ni, and Al on the {100 
face of rock salt. He tabulated these temperatures 
along with the lattice constant for each of the metals 
as is reproduced in table 1 without advancing an 
explanation for the differences in temperature. The 
most marked of these differences is that between 
silver and gold for which the lattice constants are 
nearly identical. Each of the metals listed has the 
cubic face-centered lattice. Three metals, Fe, Co, 
and Cr, for which Briick was unable to reach the 
epitaxial temperature due to the fact that the rock- 
salt substrate sublimes at 540° C are not included in 
the table. 


with their cube edge | 








Fiaure 7. Silver on the | 111} face of rock salt. 
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2.2. Epitaxial Temperature 


As has already been stated, one of the conditions 
for laying down an oriented deposit of a metal from 
the vapor phase onto a substrate is that there exists 
a minimum substrate temperature, the epitaxial 
temperature, below which perfect orientation is not 
observed. 

Uyeda [15] has extended the investigation to a 
variety of substrates. 

The data of Briick and those of Uyeda for the 
deposition of silver from tabulations given by Uyeda 
are listed directly in table 2. 

From this set of data it appears that the epitaxial 
temperature for orienting a given face parallel to the 
substrate (for example, the (100) face of silver 
parallel to the {100} face of galena, iron pyrite, and 
rock salt) is not very different for different substrates. 
The temperature for laying the {100}, {110}, or 
the {111} face of a given material parallel to a given 


TABLE | substrate, however, is different. Briick observed 
that, in general, the most densely populated plane 
Deposit Lattice Epitaxial tends to go parallel to the substrate at a relatively 
ww | Swen | eee low temperature; then comes the next most densely 
- populated at a higher temperature, and soon. The 
Ni 3. 52 70 most densely populated plane for cubic face-centered 
Pa 3 88 230 metals is the {111} plane; for cubic body-centered 
Al oo - metals, it is the {110} plane. Lau 
Ag 108 10 The question of what determines the epitaxial 
temperature for a given element remains as yet 
TABLE 2 
Substrate . — Orientation Fn ance 
°C 
Molybdenite C7 (0001) srosy!|(11 1) ae; (1010) mos, [110 ag <20 
Zincblende B3, B4 (110) zes!!(110) ae; [110]za5 [110]ae 80~ 100 
Galena. Bl (100) pps! (100) ag; [100}p.s [100]. 100~150 
Iron pyrite C2 (100) res, ||(100)4e; [100] re8,!|[100]r. ~150 
Rock salt_. Bl (100) wact || (100) ag; 


[100}~wact [100 ]a¢ 150 
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unanswered. Smollett and Blackman [16] 
stated that is is hard to understand the temperature 
dependence of orientation in terms of critical misfit 
alone. Although it 
epitaxial temperature may serve the purpose of 
expanding the lattice, such expansion can hardly 
play an important role since Eppelsheimer and 
Penman [17] have shown that the copper lattice 
constant inereases only slightly over 1 percent 


have | 


has been suggested that the | 


between 18° and 770° C, which is negligible in com- 


parison with observed misfit in epitaxy. Neuhaus 
18] has associated the epitaxial temperature with 
the lattice energy of the deposit and has concluded 
that the oriented deposits are not the result of 
adsorption (partial energy exchange) but are the 
result of a true two-dimensional reaction (complete 
energy exchange). 


3. Process of Epitaxy 


Since the first experimental work on oriented 
deposits was done, attempts have been made to 
explain the observed orientations. At the end of a 
very extensive experimental investigation, Royer [4] 
summarized his observations in three empirical 
rules. These were (a) that within a certain tolerance 
the geometric distribution of lattice sites in the plane 
of the substrate and in that of the deposit should be 
identical or the distribution in one should be a 
simple multiple of that in the other, (b) that the 
deposit ions, which take positions that substrate 
ions would have taken had the substrate continued 
to grow, should have the same polarity that the 
substrate ions would have had, and (c) that the type 
of binding between elements of the deposit should be 
identical with that between elements of the sub- 
strate. He wrote the third rule with hesitancy since 
he was not certain that the evidence that was avail- 
able really proved it. Since then, Sloat and Menzies 
(6] found orientation of ammonium bromide and of 
cesium chloride deposited from aqueous solution onto 
a single crystal of silver. They interpreted the 
observed orientation by assigning positive charges of 
not necessarily effectively unit value to the atoms of 
the silver crystal. They assumed that these charges 
alined the negative ions that form lattice planes in 
the salts. Johnson [19] found epitaxy between 
sodium chloride deposited from solution onto ori- 
ented films of silver, which had the {100} and the 

111} planes parallel to the surface, and on poly- 
crystalline silver sheet. He concluded that Royer’s 
third rule should be modified or abandoned. 

One of the earliest of the theoretical treatments of 
epitaxy was that of Dixit [20]. He considers that 
the first layer of adsorbed atoms in random lateral 
motion on an unattractive (amorphous) solid surface 
may be represented as a two-dimensional gas obey- 
ing essentially the ideal gas equation 


xA=RT. (1) 


Here x is the two-dimensional pressure, A is the 
surface area per adsorbed molecule, F is the gas 
constant, and 7’ is the temperature of the surface. 
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The two-dimensional pressure 7 is given by 
7=Y0—-7; (2) 


where yo. is a constant and ¥ is the surface tension. 
By substituting the value of y given by the Ramsay- 
Shields equation, and by substituting successively 
the area per molecule in the (100), (110), and (111) 
plane, he succeeded in showing that the temperature 
at which any of these orientations appear is directly 
proportional to the area A in this plane. By means 
of additional assumptions, he was able to calculate 
appearance temperatures for these planes that com- 
pared very favorably with those observed. Although 
his treatment applies to an unattractive substrate, 
the difference in the appearance temperatures of 
these planes should be the same on any substrate 
if the assumption of a two-dimensional gas can be 
applied to the diffusion process on the substrate 
surface. 

Briick was aware of Dixit’s work at the time that 
he discussed his own experimental results, but he 
regarded the treatment as inadequate although it 
was in line with his own observation that the densest 
plane of the deposit tends to go parallel to the sub- 
strate first. He noticed the seeming preference for 
misfit in the orientations of nine metals on the 
cleavage face of rock salt that he observed experi- 
mentally, and attempted to explain the results by 
saying that the orientation appeared to be such as 
to make the sum of the distances of the substrate 
ions of one species from the deposit atoms a minimum 
inacell. He says that the arrangement takes place 
in such a manner that for the range of an elementary 
cell the sum of the distances of the atoms of the 
uppermost layer of the substrate from those of the 
deposit lattice will be as small as possible. Briick’s 
rule has been supported by Uyeda, as we shall see 
in a moment, but it has been criticized both by 
Bri and Gharpurey [14] and by Thomson and 
Cochrane [21]. 

Thomson and Cochrane point out that what would 
be a favorable orientation in the first cell if Briick’s 
rule were obeyed would be an unfavorable one later 
on. This is, of course, true, and in this respect 
Briick’s rule is inadequate although it points in the 
right direction. If Briick’s rule is extended, how- 
ever, it can be used to define the most favorable 
orientation in terms of attractive bonding energy to 
the substrate for the largest area of monolayer 
deposit lattice that is capable of forming over the 
substrate in question. This will be discussed in 
more detail in section 5. Bra and Gharpurey dis- 
miss Briick’s explanation as unconvincing in that 
among ions of the same polarity the distances have 
rather to be large than small [13]. However, it 
does not appear to be necessary to minimize a 
distance-sum over ions of the same polarity. It is, 
for example, not necessary to conceive of the (111 
face of rock salt as being made up of sodium ions. 
To be absolutely realistic, the best {111} face that 
one could cut experimentally would undoubtedly 
house both ionic species. It appears that the com- 
ment of Bri and Gharpurey must have been directed 


toward the fact that Briick, in considering the cleav- 
age face of rock salt, imagined the sodium ions de- 
leted from the lattice, so that a simple face-centered 
cubic structure of chloride ions remained. ‘t was 
in a cell of this modified lattice that his rule applied. 
An extension of Briick’s rule to include antibonding 
as well as bonding energy will be discussed at 
length in section 5. 

Thomson [22], and Thomson and Cochrane [21] 
thought that steps and cracks on the surface are 
often a determining factor in producing epitaxy. 
Wilman [23] rejects the explanation of Thomson 
and Cochrane. He states that rotational-slip results 
show that the existence of angles or edges of cleavage 
steps, whether these are macroscopic or only mona- 
tomic in height is quite unnecessary for epitaxial 
orientation to occur, but that if such steps or cracks 
are present, they may often determine the location 
or habit of the deposit crystals. 

Uveda [15] regards the cracks and steps as con- 
tributing only a disturbance to the ideal case. He 
regards the preferred orientation as being induced 
by the periodic field of surface foree, and has ad- 
vanced the following two assumptions in terms of 
which to discuss the oriented growth of crystals: 
(1) The configuration of atoms in a crystal bud 
formed on a base is almost the same as that formed 
in a free space, and (2) the interaction between the 
base and the crystal decides the orientation of 
crystals. He concludes that the orientation should 
be fixed while several atoms are gathering around a 
nucleus (compare Briick’s orientation in an ele- 
mentary cell) and states that this point of view 
supports Briick’s rule of orientation and explains a 
few phenomena concerning orientation qualitatively. 
He has, however, not made clear what he means by 
interaction between the base and crystal. Kainuma 
(24) has found evidence in favor of Uveda’s con- 
clusion. 

Kubo and Miyake [25] suppose that the vaporized 
metal atoms, after wandering on the surface, finally 
arrive at the foot of a nucleus and it is from here 
the crystal growth commences. 

The mathematical treatment of Frank and van der 
Merwe [7, 8] has already been referred to. Their 
theory is based on the properties of a one-dimensional 
dislocation model. It consists of a row of identical 
balls connected by identical springs having force 
constant uw. The balls are acted on by a force that 
varies periodically with the position on the substrate. 
They let a be the wavelength of the potential due to 
the substrate and W/2 be its amplitude. The 
potential energy is represented by the first harmonic 
term (amplitude W’/2) in a Fourier series. A differ- 
ence may exist between the natural spacing 6 of 
the balls and the wavelength a of the substrate field. 
From the theory for the model it turns out that 
when the natural spacing 6 of the balls differs from 
the wavelength a of the substrate field, the lowest 
energy state of the system is one in which there are 
no dislocations (small regions in which the balls 
would be spaced more closely or less closely than 
normal) up to a certain critical value of the misfit 





1/ Py, where 


1/Py)=(b/a 2/ri, 


1) critical 


and 
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w-(Si7) : 


By using Lennard-Jones forces, and by assuming 
interactions with other deposit atoms and with 
substrate atoms are similar, they find the value of 
/, to be about 7 and hence the critical value of th 
misfit to be about 9 percent. Below this critica! 
misfit the monolayer in its lowest energy state is 
deformed into exact fit with the substrate. Above 
the critical misfit it is only slightly deformed on thy 
average, but it contains many dislocations. Up to 
the larger critical misfit of 14 percent the monolayer 
can be deposited metastably in exact fit on the 
substrate at temperatures that are sufficiently low 
Completely oriented overgrowth is to be expected 


only if the first monolayer can form over the entire 


surface under subcritical conditions. 

Smollett and Blackman [16] have investigated the 
stability of a strained two-dimensional layer. They 
evaluated the lattice constants of a two-dimensional 
alkali halide layer strained to fit on an alkali halide 
crystal substrate. For a series of cases in which 
oriented overgrowth is found experimentally they 
found instability in the initial deposited layer if it is 
strained to fit the substrate. They conclude that 
experimentally observed orientations cannot result 
from an initially uniformly strained layer such as is 
indicated by the Frank-van der Merwe _ theory. 
They further conclude that the initial lattice must 
be distorted in the sense that displacements occur, 
but that the lattice constant of the deposit material 
remains more or less unchanged. 

In analyzing the forces that operate in epitaxy 
from solution, Sloat and Menzies |6] state that the 
most important are the mutual attraction between 
the solute ions and the substrate ions, and the mutual 
attraction between oppositely charged solute ions 
Following Lennard-Jones and Dent [26] they further 
analyze the first force mentioned into four parts 
(1) Electrostatic attraction of the ionic charge to 
the valency charges of the substrate; (2) the force 
of attraction between substrate ions and the dipole 
which results from the polarization of the attracted 
ion by the surface field; (3) the attractive force 
resulting from the polarization of the surface tons 
by the attracted ion; and (4) the van der Waals 
attraction between neutral atoms. The van der 
Waals attraction is operative at longer distances, but 
at short range the electrostatic force is much greater 
in magnitude. They conclude that it is first the 
van der Waals force which is operative; then the 
electrostatic. They ascribe the final capture of an 
ion, and the fixing of the ion into the space lattice 
of the substrate, to the electrostatic forces. 

It appears that a similar mechanism must operat: 
in the epitaxy of metals deposited from their vapors 
although Rhodin [27, 28] has found rather surprisin: 
agreement on the assumption of van der Waal- 
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horces. 
he deposition of a monolayer was conductivity across 
the film. It is hard to believe that a continous 
nonolaver of metal is ever produced because of the 
strong ‘cohesive forces of the deposit. Rather it 
vould seem that aggregates of a number of layers in 
depth eventually grow together. One can then 
hardly talk about the binding force of the first laver 
from a conductivity measurement. 


4. A Model for Metal—Rock Salt Pairs 


In order to establish a basis on which prediction 
of metal-substrate pairs for vapor-phase deposition 
and of their epitaxial temperatures can be made, it 
is necessary to know the process involved. If this 
process can be understood for one substrate, it can 
he generalized to others. With this view in mind the 
special case of the epitaxy of metals on a rock-salt 
substrate has been considered [28a] because a larger 
amount of experimental data exists for this substrate. 

We shall start with the following assumptions: 

|) Orientation must result from the interplay of 
forces, and (2) the preferred orientation must be the 
configuration of lowest energy. If we assume that 
the first layer of deposit-metal exists as atoms, the 
forces, in addition to van der Waals forces, that are 
involved are electrostatic forces acting between the 
substrate ions and the dipoles induced by the sub- 
strate ions in the neutral metal atoms. If, on the 
other hand, we consider that the atoms of the first 
laver of deposit-metal become ionized, they are 
electrostatic forces acting between the ions of the 
substrate and deposit-metal ions on the surface. 
When neutral metal atoms arrive on the surface of 
the substrate, forces of the first kind are operative. 
If, after arrival on the surface, the deposit-metal 
atoms permanently lose one or more electrons, forces 
of the second kind exist. A true two-dimensional 
reaction can then occur between the deposit-metal 
ions and the anions of the substrate so that at the 
interface the three lavers—substrate salt, deposit- 
metal salt monolayer, deposit-metal—exist. The 
deposit-metal salt that forms as a result of this re- 
action must be sufficiently mobile to diffuse on the 
surface possibly in the form of symmetrical units, or 
lattice cell fragments, such as discussed by Zener [29] 
until larger static aggregates are formed. It would 
appear that the epitaxial temperature should be re- 
lated not only to the mobility of this salt, which in 
the form of larger aggregates is eventually firmly 
anchored to the substrate by true electrovalent bonds, 
but that the epitaxial temperature is also a function 
of the ionization potential of the respective metals. 
Thermal ionization at the expitaxial temperature is, 
of course, out of the question because of the vanish- 
ingly small value of the Boltzmann factor. But the 
combined effect of thermal energy and the availabil- 
ity of electron traps on the surface could be sufficient 
as we shall see in a moment. If one plots the first 
ionization potential of each of these metals against 
the epitaxial temperature, there is no relation be- 
tween them. If, on the other hand, one plots the 
ionization potential of the most commonly encoun- 


However, his criterion of the completion of | 


tered ionization state, the trend is along a line (see 
fig. 8, and reference [32]). For these metals the only 
valence states used that are higher than the lowest 
ionization state possible for the atom in question 
under usual conditions are Cuu, Au mt, and Ni 1. 
We shall return shortly to a discussion of the reason- 
ableness of assuming that these will be the most 
stable valence states to be expected for these metals, 
and we shall consider first whether or not the process 
that has been advanced is energetically feasible. 

In considering what could constitute trapping sites 
for the electrons, one might at first sight fix upon the 
sodium ions. That this is not a logical choice, how- 


| ever, is evident from the faet that if rock salt is 





heated in sodium-metal vapor, F-centers are formed 
in preference to the alternative that the lattice should 
house sodium metal interstitially. The F-centers 
(Farbzentren, or color centers) are vacant chloride 
ion sites, which due to their effective positive charge, 
have trapped an electron whose orbit then extends 
over the six surrounding sodium metal ions. They 
may be regarded as a vacant negative ion site at 
which one of the nearest neighbor positive ions has 
gained the electron and become a neutral metal atom. 
If we write .M/ for the deposit-metal, S for a vacant 
negative ion site, and F for a color center, then we are 
concerned with the energy for the two-step process, 


M°—ne->M"** (5) 


Step (1): 


Step (2): nS+-nenF. (6) 
The second step in sodium chloride involves about 
2.7 ev. This is the peak of the characteristic absorp- 
tion band of the F-centers [30]. The metal atom, 
let us consider it to be silver, is first adsorbed from 
the vapor onto the cleavage face of the rock salt. 
After a certain period of time, during which it is 
adsorbed on the surface, loss of an electron takes 
place. Since the silver atom is on the surface of the 
sodium chloride, the energy required is less than the 
gas phase ionization potential of silver. That is, the 
ionization potential is reduced because of the dielec- 
tric constant of the salt. Because a thermal rather 
than an optical process is involved, the static dielec- 
tric constant A is in question rather than K optical 
arising from the distortion of the electron clouds. 
The whole value of K would be effective only if the 
process took place below the surface of the crystal, 
that is, if the adsorbed atom were completely dis- 
solved in the rock salt. If the cleavage plane were 
atomically smooth, one might reasonably assume that 
the dielectric constant should be reduced by some- 
what less than half. Actually the factor 0.44 has 
been used as this gives the lowest positive value for 
AF in the case of silver for which the data are known 
most accurately. In addition to the energy available 
from eq (6) there is also available the energy of plac- 
ing a Ag* ien over a Cl~ ion on the cleavage eae. 
If the cleavage plane were atomically smooth, this 
energy would be given by [31] 


0.066 é 
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where ¢ is the electronic charge, and where r is the 
Ag*—Cl- distance of separation. The first ionization 


potential of silver [32] is 7.574 ev, the dielectric | 


constant of rock salt [30] is 5.62, and if we take the 


Ag*—Cl> distance to be the sum of the ionic radii | 
(Goldschmidt radii are Ag 1.13 A, Cl 1.81 A), 
then 
, 7.574 0.066 67 
) -—2.7- 0.04 ev. (S) 
a | i r fa 


In the case where a doubly charged ion is produced, 
as in Pd**, an additional energy is involved due to 
the polarization of the surface ions by the double 
positive charge. Wright [33] gives the potential 
when an electron is added at a positive ion site in 
sodium chloride as —2.5 ev. It would appear that 
on the surface, although the Pd** ion has only one 
nearest neighbor Cl~ ion, it sees all of the ions in the 
plane on which it is standing, so that the polariza- 
tion energy is a large fraction of that which results 
when an electron is given to a positive ion in the 
interior of the salt. A reasonable guess might be 
0.7 of this value or 1.75 ev. In the case of forming a 
doubly charged ion, two F-centers are involved and 
the energy for placing the ion on a Cl~ ion in the 
plane is doubled. Because the second ionization 
potential of palladium is 19.42 ev [34], the energy for 
producing the Pd** ion is given by 


19.42 oo 2) 


~— 2(2. 0.01. 
044K ~' 


AE 2(0.066 e2/r)—0.7(2.5) 


(9) 


There is little point in estimating the energy change 
for a triply charged ion such as Au*** because of the 
uncertainty both in the value of the ionization 
potential and in the polarization energy. 

This value of AZ, however, does not determine the 
equilibrium constant. It is really the thermo- 
dynamic AF} or energy change for the reaction at the 
absolute zero of temperature. The equilibrium 
constant A is given by 


—InK=AFY/RT, (10) 


where # is the gas constant, 7 is the absolute tem- 
perature, and AF” is the standard free-energy change. 


AF’ =AE} +A ( a) T. (11) 


For the ionization process postulated, the second 
term on the right for silver is 


(PES FE; 
r| ( T ) +f * gg iy 


F-center 
F°—E° F°— Ee 
(77*)_-(77*)_.. } 


This information is not available and consequently 
the equilibrium constant cannot be calculated. 








If it is conceded that ionization of the depo 
atoms does occur on the surface of the substrate, » 
the evaluation of the energies involved appears 
indicate that such a process is feasible, then it mu 
| also be conceded that the epitaxial deposits consi 


of the layers: rock salt, deposit-metal salt, deposi 
metal, of which the deposit-metal salt layer is 
monolayer. If this is indeed true, we must expe; 
in every case in which epitaxy is observed that ther 
is epitaxy between rock salt and the deposit-meta 
salt as well as epitaxy between the deposit-meta 
salt and the deposit-metal, and we must consider th, 
lattice fit between the salt lattices, and that between 
the deposit-metal salt lattice and the deposit-meta! 
lattice rather than the fit between the metal lattic: 
and the rock-salt lattice as has been done so far. 

Following this point of view, when the silver meta! 
itself deposits, it does so on an oriented surface of 
silver ions. Here, three cell lengths of the silver 
chloride lattice (16.62 A) are closely equivalent to 
four cell lengths of bulk silver (16.32 A). The ease 
with which the films can be floated from the salt 
cakes seems to indicate that they are not bonded to 
the substrate cell by cell such as an additional layer 
of sodium chloride would be. Rather, there must 
literally be tunnels through which the water can 
flow under the film. 

Some very interesting experimental work has been 
done recently by Johnson [35] in regard to epitaxy 
of salts on metal surfaces. He found epitaxy be- 
tween rock salt and the following metals: Cu, Ag, 
Au, Zn, Pb, Sb, and Bi, when he evaporated a salt 
solution on the metal. ‘Two metals known to exhibit 
epitaxy when deposited from the vapor onto a rock- 
salt cleavage plane were negative to Johnson's test. 
These were nickel and aluminum. While one might 
ascribe the negative result in the case of nickel to 
the large misfit between the nickel lattice and the 
rock salt lattice, one could hardly advance this argu- 
ment in the case of aluminum where the lattice con- 
stant is very close to that of silver and gold. The 
answer may then possibly be that it is not necessary 


| that there should be epitaxy between rock salt and 


the metal itself, but that it is only necessary that 
there should be epitaxy between rock salt and the 
metal salt; and again, epitaxy between the metal 
salt and the metal as has just been indicated. John- 
son explained his results by assuming the presence of 
an oriented oxide layer on the metal. Where the 
oxide layer was known to be amorphous, orientation 
did not occur. In the case of silver and gold no 
oxide layer was considered to be present. In this 
connection Johnson states that it is possible then 
that the oriented growth of sodium chloride on silver 
did involve an intermediate thin layer of silver 
chloride. Silver produced a much better orientation 
of the sodium chloride than did the gold. In regard 
to this, Johnson states that the difference in behavior 
of silver and gold must be attributed to differences 
in their chemical activity. Briick [11] also found 
that silver was a more efficient substrate materia! 
than gold as is evident from the epitaxial tempers 

tures that he found for deposit-substrate pairs give: 

in table 3. 
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Substrate Deposit 
Ag Au 
« ( 
NaCl 1%) 400 
Ag me oon 
Au 2m) 


Sloat and Menzies [6] have shown that sodium 
chloride can orient silver chloride and Wilman [36] 
has shown that silver can orient silver chloride. 
Consequently, orientation is possible through the 
three layers: NaCl, AgCl, and Ag. In the list of 
positive and negative cases of epitaxy given by van 
der Merwe [8] it is indicated that copper does not 
orient cuprous chloride. It remains to be shown 
that copper does orient cupric chloride since accord- 
ing to the plot of epitaxial temperature against 
ionization potential shown in figure 8 one is led to 
believe that the divalent salt forms. Only by obtain- 
ing similar data for the other metals for which the 
epitaxial temperature on rock salt is known, can it 
be shown that the valence state indicated by the 
epitaxial temperature is the only ionization state for 
which the metal salt will orient both the substrate 
salt and the metal deposit. 

Of the metals so far considered (data of Briick [11]) 
the only valence states used that are higher than the 
lowest state possible for the elements in question 
under usual conditions were Cu u, Au ru, and Ni m1. 
In determining whether or not chlorides in these 
valence states would be stable at the epitaxial tem- 
peratures required, one must consider both the 
stabilizing effect of the sodium chloride and the fact 
that a monolayer is in question. Morris [37], for 
example, has found that auric chloride solution 
contaiming sodium chloride is stable to 450° C but 
auric chloride by itself decomposes at a much lower 
temperature. What effect stabilization due to the 
sodium chloride may have in determining the valence 
state that should be observed can only be surmised. 
The chemistry of thin layers, let alone monolayers, 
has been investigated very little. In regard to the 
use of the higher oxidation state for gold, it is note- 
worthy that aurous chloride undergoes dismutation 
with ease to produce auric chloride and gold. 

The limitations of the graph shown in figure 8 
must be clearly borne in mind. Each of the six 
metals on it has the cubic face-centered lattice type. 
The epitaxial temperature for each is for deposition 
of the {100} plane parallel to the {100} plane of 
rock salt. It is to be noted that for face-centered 
cubie lattices, the orientation with the {111} plane 
parallel to the substrate should occur at a lower 
temperature than that for orientation of the {100} 
plane parallel to the substrate. However, the 
orientation with the {111} plane parallel to the sub- 
strate in many cases never appears. Briick, for 
example, observed that this was true for copper, 
palladium, nickel, and cobalt. He found that the 


111} orientation parallel to the substrate did appear | 
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for gold and aluminum and that there was some 
indication of its appearance for silver. The fact 
that it does appear in aluminum, and must be forced 
out by increase in temperature, may account for the 
poorer fit of aluminum on the graph. _ If the epitaxial 
temperature for laying down the {111} plane of 
aluminum parallel to the substrate (which may be 
taken to be about 350° C since above this value the 
(100} plane starts to deposit) were used for alumi- 
num, this metal would have a much better fit on the 
graph. This, too, would be in line with the assump- 
tion that the metal must be ionized to be oriented at 
all. Gold would also be displaced to the left on this 
reasoning. 

In the discussion of Dixit’s investigation in section 
3, it was pointed out that although his treatment 
applies to an unattractive substrate, the differences 
in his appearance temperatures should be correct for 
any substrate if the assumption of a two-dimensional 
gas can be applied to the diffusion process on the 
substrate ‘aon For silver, Dixit finds the {111} 
appearance temperature to be 659° C and that of the 
{100}* face 801° C. The difference is 142 degrees. 
Brick reports that silver can be oriented fairly well 
at room temperature. At 150° C he obtained single 
crystal orientation with the {100} face parallel to the 
substrate, but sometimes, in the foils nearly oriented 
to a single crystal, there was a light superimposition 
with the {111} structure. Assuming that the first 
observed orientation was with the {111} plane parallel 
to the rock salt, the temperature range between the 
two orientations is about 130 degrees. For alumi- 
num on an unattractive substrate Dixit found for 
the {111} appearance temperature 264° C, for the 
{100} appearance temperature 346° C, and for the 
{110} appearance temperature 602° C. The differ- 
ence between the appearance temperatures for the 
{100} and {111} planes is 82 degrees. In regard to 
aluminum Brick states that above 350° C the single 
crystal drives out gradually the {111} structure until 





Dixit, who was considering interplanar spacings, used the notations {200} and 
{222} instead of {100} and {110} 
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it only exists at 440° C. The aluminum film at 
440° C has the {100} plane parallel to the rock-salt 
cleavage plane. This is a temperature range of 90 
degrees, which again compares favorably with the 
difference in Dixit’s calculated values. This agree- 
ment appears to indicate that the same assumption 
of a two-dimensional gas may be applied to attractive 
substrates. 

There appear then to be two temperature depend- 
ent processes: (1) the ionization process and (2) the 
process by which the {111}, {100), and {110} planes 
orient parallel to the substrate. The second process 
is a complicating factor in the first. The answer to 
the question as to why orientation with the {111} 
plane parallel to the substrate does not make its 
appearance for every metal is not clear. It must be 
energetically forbidden. The answer may be simply 
that the temperature required to produce ions that 
can be oriented at all in sufficient number is already 
in excess of that at which the {111} plane would 
orient parallel to the substrate if ions had been avail- 
able. Then the next most densely populated plane, 
the {100} plane for the face-centered cubic structure, 
would show up. 


5. Retrospect 


In this section we return to the experimental re- 
sults that have been cited for the deposition of metals 
from their vapors onto a rock salt substrate to find 
if the ionization process that has been advanced can 
explain the orientations which have been found. At 
its epitaxial temperature, Briick found that silver 
deposits on the cleavage face of rock salt in the way 
indicated in figure 4, although, if there were a relative 
rotation of 45 degrees, the orientation designated as 
not found in figure 4, the misfit would be very much 
less. It is interesting to observe from the two- 


dimensional sketches that the misfit orientation that. 


is actually found places a chloride ion in closest prox- 
imity to a silver atom (or silver ion, in terms of the 
ionization process discussed) and leaves the sodium 
ions with no silver close by. If the silver cube 
diagonal did lie along the rock-salt cube edge (the 
orientation involving best fit, but which is not found) 
a sodium ion would make a close approach to a silver 
atom. If the silver atom is really a silver ion, this 
would involve repulsion. If silver ions, rather than 
silver atoms, are involved, the orientation observed 
is in accord with Royer’s second rule, which says 
essentially that you will not get epitaxy if ions of the 
wrong polarity have to be brought together. This 
rule is for the orientation of charged ions depositing 
out of solution. The fact that it indicates as impos- 
sible the orientation that is actually not found, but 
which one would expect to find on the basis of best 
fit, lends support to the idea that charged ions are 
involved. The orientation that is found appears to 
try to establish as closely as ible a monolayer 
lattice of silver chloride using the anions of the sub- 
strate. The lattice constant of silver chloride 
(a=5.545 A) is nearly the same as the lattice con- 
stant of sodium chloride (a=5.628 A). A close 





approach to this hypothesis was made by Dral 
[38] who concluded that the epitaxy is such that | 
nearest possible coordination fit between the sil. 
chloride (sodium chloride type structure) and sil 
structures occurs across the interface. The conclus; 
of Drabble is cited by Bré and Gharpurey [13, 

in regard to their own experimental results. 

In the deposition of silver on the {100} face of roc k 
salt at 450° C one sees again in the orientation not 
observed that sodium ions would be brought close jo 
silver ions and the situation is the same as that just 
discussed (see fig. 5). It was pointed out in regard 
to the deposition of silver on the {110} face of rock 
salt that if the silver were rotated through 90 degrees, 
the fit would be nearly perfect, but that this orienta- 
tion is not found. The orientation of best fit which 
is not found experimentally would again involve a 
close approach of the silver ions to the sodium ions 
(see fig. 6). 

As was also pointed out previously, Briick found 
that whereas the cubic face-centered metals grew 
with the {100} plane resting on the (100) of rock 
salt and oriented with their cube edges parallel to the 
rock-salt cube, the cubic body-centered metals grew 
on the contact surface turned at 45 degrees to the 
rock salt cube. 
case of a cubic body-centered metal, there would be 
no interference from repulsion of sodium ions in de- 
positing such a metal with its cube diagonal along 
the sodium chloride cube edge because there would be 
no deposit-metal ions in these positions. 

We see that although Briick’s rule points in the 
right direction, it needs to be extended. If we con- 
sider that electrostatic forces are involved, we must 
say that the arrangement is such as to give maxi- 
mum bonding to the substrate for the largest aggre- 
gate that is capable of forming. One must then 
sum the attractive bonding energy, and the repul- 
sive antibonding energy between each silver ion and 
its four nearest neighbor substrate ions of each spe- 
cies over the area that is defined by the points where 
exact fit is obtained between the substrate salt lat- 
tice and the deposit-metal salt lattice. If the sum 
of the attractive energies (—) added to the sum of 
the repulsive energies (+-) is negative, the arrange- 
ment is stable, but if it is positive, the arrangement 
in question should not exist. It should be possible 
to explain all of the orientations that have been 
found by determining these sums graphically both 
for these orientations and for other possible orienta- 
tions which are not observed. It should also be pos- 
sible to explain the laying down of the {100} plane 
and the {110} plane of the face-centered cubic met- 
als parallel to the substrate at higher temperatures 
as higher energy configurations which are neverthe- 
less stable (i. e., the sign of the sum is negative). 

A small square region (15 cell lengths by 1) cell 
lengths) of silver-chloride lattice imposed on the 
rock-salt lattice has been mapped. A silver ion was 
placed directly over a chloride ion at the origin and 
actual measurement of Ag*t—Cl~ distances, and 
Ag*—Na* distances was carried out in the first quad- 
rant only. At a cell having coordinates of seven 





From figure 4 it is clear that in the: 
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cell lengths from the origin along both the ordinate 
and abscissa the misfit has become about 10 per- 
cent. The cell at this distance from the origin is 
shown in figure 9. Each silver ion in the area con- 
sidered has four chloride ion neighbors and four 
sodium ion neighbors in the plane below the silver 
ion monolayer. Planar distances, which were meas- 
ured, were converted to actual distances by means 
of a graph of planar distance versus actual distance 
for which the actual distances were computed by 
using a 20-inch log-logvector slide rule. The bond- 
ing energy is the sum of e*/r, where the distances, r, 
are those of the silver ions from their chloride ion 
neighbors. The antibonding energy is the sum of 
Jr. where the distances, r, are those of the silver 
ions from their sodium ion neighbors. In the count- 
ing process, each distance is counted twice except 
for the Ag*—Cl~ distance at the center of the cube 
face. Where distances were counted twice, half the 
sum was taken. The result for the first row along 
the abscissa for seven and a half cells was 7.092 &? 
for the bonding energy, and 4.821 e? for the anti- 
bonding energy. The net bonding energy was 2.271 e’. 
For the seventh row above the abscissa where et 
the maximum misfit reached about 10 percent, the 


NeT bonding energy had decreased only to 2.259 e?. 


At this rate of decrease a very large area could be 
covered before the antibonding energy became equal 
to the bonding energy. In this connection, Hass [39] 
has found that very thin silver films deposited on 
collodion exist as isolated granules 100 A to 150 A 
in diameter, but that thin silver films deposited on 
rock salt do not contract as strongly into isolated 
granules. Measurement of distances for only the 
first cell at the origin with the silver chloride turned 
at 45 degrees with respect to the rock salt cube gave 
an energy already slightly antibonding. Hence at 
a nucleus consisting of the first cell, this orientation 
became impossible. 

To find if this basis for predicting epitaxial tem- 
peratures has any real value, Suddeth and McCraw 
have carried out evaporations of metals in a vacuum 
on rock salt at temperatures which their common 
ionization potential would indicate to be about 
correct on the graph of epitaxial temperature 
against ionization potential (see fig. 8). In the case 
of lead, which according to its ionization potential 
should have oriented at about 200° C, a slight orien- 
tation was found at 225° However, the results 
with lead were not reproducible even among separate 
deposits on the same cleavage surface which were 
made simultaneously and in no case was complete 
single crystal orientation obtained. It is possible 
that the misfit was too large and that better results 
would have been obtained on potassium chloride, but 
there was no potassium chloride single crystal on 
hand at the time to use for the purpose. Lead is 
cubic face-centered with a=4.93. Lead chloride is 
orthorhombic. The cell dimensions are: a=9.030, 
b=7.608, and e=4.525. 

According to the graph of epitaxial temperature 
against ionization potential, antimony, with a com- 
mon valence of three, should orient at about 380° C. 
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Ficgure 9. Misfitted silver ion lattice on rock salt. 


@. Silver ion; X, sodium ion; ©, chloride ion 


Suddeth and McCraw have found partial orientation 
for antimony from 250° to 410° C. They found 
appearances in the diffraction patterns similar to 
those observed by Briick for aluminum and gold, 
which he interpreted as the transition from deposition 
with the { 111) face parallel to the substrate to deposi- 
tion with the {100} face parallel to the substrate. 
Suddeth and McCraw have not as vet obtained a 
perfect single crystal of antimony having only the 
100} face parallel to the rock salt. Antimony has 
the arsenic-type lattice. For the antimony cell 


| a=4.50, and the nearest neighbor distance d is 2.87. 


Data for antimony trichloride do not seem to be 
available. For SbI; in a bimolecular cell, a) =8.186 
with a=54°14’ and in a hexagonal pseudo cell 
a, = 7.466 and c= 20.89. 

According to the graph of epitaxial temperature 
against ionization potential, platinum for an ioniza- 
tion state of two should orient at about 250° C, 
Thirsk [40], by employing evaporation from a hot 
platinum wire, has obtained orientation of Pt on 
potassium chloride at the epitaxial temperature of 
250° C with the Pt-cube edge parallel to the cube 
edge of the substrate. The | 100) face of potassium 
chloride was used. The epitaxial temperature on 
rock salt will not be very different judging from the 
data of Briick where he used potassium chloride as 
substrate as well as sodium chloride. The data for 
the three metals for which he found the epitaxial 
temperature on both substrates are given in table 4. 


TABLE 4 
Gashe . Metal Epitaxial 
Substrate deposit temperature 
°c 
NaCl Ag 1%) 
ACi... Ag 150 
NaCl Pd 2%) 
KCl Pd 25) 
NaCl Ni 370 
KCl Ni 410 
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Briick pointed out in regard to these data that the 
temperature varied in the case of nickel where the 
misfit was large. 

Consequently, it seems fairly certain that when an 
adequate method is used for depositing platinum on 
rock salt, it will be found to orient at a temperature 
close to that predicted by the graph of epitaxial 
temperature against ionization potential. 

Suddeth and MeCraw made an evaporation of tin 
at the epitaxial temperature which the graph would 
predict for the valence state of two. The result was 
negative with pure tin but showed considerable orien- 
tation in an impure sample which probably contained 
lead 


6. Summary 


From an analysis of experimental data on the 
formation of oriented deposits of metals from their 
vapors in a vacuum onto salt substrates, it appears 
reasonable to suppose that a surface ionization proc- 
ess may explain the existence of metal ions that are 
brought into preferred orientations by attractive 
electrostatic forces that act between them and the 
subtrate ions. A relation between ionization poten- 
tial and the epitaxial temperature at which the 

100, plane deposits parallel to the (100) plane of 
the rock-salt substrate has been found which may 
be useful as a basis for making predictions as to the 
temperature at which orientation will occur. The 
results appear to indicate that the sequence of 
layers —substrate salt, deposit-metal salt monolayer, 
deposit-metal— exists. 

In making a prediction of a substrate-metal pair 
and of the temperature at which orientation of the 
100) plane of the deposit parallel to the {100 
plane of the substrate should occur: (a) It is neces- 


sary to consider the lattice fit between the substrate 


salt and the deposit-metal salt, and the lattice fit 
between the metal salt and the metal. (Structure 
data for many of the metal salts in the anhydrous 
state are at present unavailable.) (b) The epi- 
taxial temperature can be roughly estimated from 
the accompanying graph by using the most common 
ionization potential of the metal in question. But 
if the metal salt corresponding to this ionization 
potential is unstable, or if this salt cannot orient on 
the substrate salt, a higher ionization state (and cor- 
respondingly higher epitaxial temperature) may be 
required. 
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Torsion of Anisotropic Elastic Cylinders by 
Forces Applied on the Lateral Surface’ 


Harold Luxenberg’ 


The classical Saint-Venant theory of torsion presents a method for the determination 
of the elastic behavior of an orthotropic elastic beam of uniform arbitrary cross section 


twisted by forces applied on the end sections, with the lateral surface free of stress 


In the 


present paper the theory is extended to include a more generally anisotropic beam (possessing 
only a single plane of elastic symmetry, rather than three mutually orthogonal planes of 


elastic symmetry, as in the case of orthotropic materials) 


built in at one end and twisted by 


forces uniformly distributed along the lateral surface 
The stresses in a beam of elliptic cross section twisted by constant tangential traction 


are obtained as an application of the theory 


l. Introduction 


The classical Saint-Venant theory of torsion 
presents a method for the determination of the 
elastic behavior of an orthotropic elastic beam of 
iniform arbitrary cross section twisted by forces 
applied on the end sections, with the lateral surface 
free of stress. In the present paper the theory is 
extended to include a more generally anisotropic 
cylinder (possessing only a single plane of elastic 
symmetry, rather than three mutually orthogonal 
planes of elastic symmetry, as in orthotropic ma- 
terials) built in at one end and twisted by forces 
uniformly distributed along the lateral surface. 

The extension of the Samt-Venant theory of tor- 
sion of cylinders by forces applied on the end sections 
to anisotropic cylinders was made by Lechnitzky 
13, 4]° in 1939; 

The torsion of isotropic cylinders byfforces applied 
on the lateral surface has been considered by Filon 
(1] in 1902, Timpe [11] in 1912, and Zwolinsky and 
Riz [12] in 1939. 

Filon and Timpe showed that the stresses and dis- 
placements in an isotropic circular cylinder twisted 
by forces applied on the lateral surface are derivable 
from a biharmonie stress function. Zwolinsky and 
Riz treated the problem of the isotropic cylinder of 
arbitrary cross section twisted by constant tangential 
forces, and showed that, in this case, the stresses and 
displacements are derivable from two stress func- 
tions, one biharmonic, the other harmonic. For a 
circular cylinder the harmonic stress function van- 
ishes, and the Filon-Timpe solution is obtained as a 
special case of the Zwolinsky-Riz solution. 

The results obtained here are more general than 
those given by Zwolinsky and Riz in two important 
respects. First, the applied forces assumed in the 
present treatment are not necessarily constant or 
tangential to the lateral surface: and second, the 
results are applicable to a wider range of elastic 
materials, including wood, which is orthotropic, and 
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certain crystalline media [5, p. 155 to 160], which 
possess a single plane of elastic symmetry. 

It will be noted that no mention has been made of 
the distribution of stresses over the built-in end of 
the cylinder. This distribution will, in general, 
depend on the mode of clamping. If these stresses 
are prescribed the exact solution of the problem 
presents considerable difficulty because of the com- 
plicated boundary conditions. 

In 1855 Saint-Venant stated as a principle of 
elasticity that two statically equivalent systems of 
forces applied to a portion of the surface of a medium 
will produce essentially the same stresses at point 
of the body not too near the region of application of 
the forces. 

If this principle is accepted, and if it is borne in 
mind that in most applications only the resultant of 
the stresses on the end of the cylinder is known. it is 
seen that a solution of practical usefulness may be 
obtained by “relaxing” the boundary conditions, and 
specifying only the resultant stress on the end section, 

The extra freedom allowed by the relaxation of the 
boundary conditions permits us to obtain a solution 
by the so-called semi-inverse method of Saint- 
Venant, which consists in making certain assump- 
tions about the distributions, while still 
leaving enough arbitrariness in the assumed stresses 
to satisfy all the conditions of the problem. 

By the use of this method it is shown that the 
stresses in the twisted cylinder are derivable from 
two stress functions, ¢(z,y) and ¥(z,y), which satisfy 
the differential equations. 


stress 
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in the region 2 corresponding to the interior of the 
evlinder, and the conditions 
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on the bounding curve C of each cross section. The 


parameters /, m, n, 6,, and by appearing in these 
equations are functions of the 13 distie constants of 
the material, NV and 7 are the normal and tangential 
components of the applied surface tractions per unit 
area, v designates the unit normal vector to C, and 
a is given by 


f[er-vxe 





i er (m se +n 3 —y(i 88 +m e¢)] 3° (3) 


In eq (3), VY and Y are the z and y components of 
the applied surface traction, and the contour integral 
is the applied twisting moment per unit length. The 
surface integral is proportional to the torsional rigid- 
ity of the cylinder. 

In the case of an isotropic beam the differential 
equations (1), simplify to 


vV’o=0, V‘y=0, (4) 


the boundary conditions (2) become 
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and a@ is given by 
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It is shown, from energy considerations, that (| 
boundary-value problems posed above possess « 
unique solution, and several methods of clatien are 
discussed. 

As an illustration of the method, the stresses in a 
cylinder of elliptic cross section twisted by constant 
tangential traction are obtained. 


2. General theory 


The first fundamental boundary-value problem of 
elasticity may be stated as follows [8, p. 72]: 

Determine the distribution of stress and the dis- 
placements in the interior of an elastic body in 
equilibrium when the body forces are prescribed and 
the distribution of the forces acting on the surface of 
the body is known. 

The second fundamental problem is similar, except 
that in this case the surface displacements rather than 
the forces acting on the surface are prescribed, and in 
the third, or mixed, boundary-value problem the 
forces are prescribed over a portion of the surface and 
the displacements over the remainder. 

The problem of the cylinder that we are considering 
is of the first kind, where body forces are absent. 

In the following paragraphs we shall review certain 
basic elements of the mathematical theory of elastic- 
ity and present the mathematical formulation of the 
first fundamental boundary-value problem. 

The state of stress in an elastic medium is char- 
acterized by the six components Tzz, Tyy, Tzzy Tyz) Tex 
Try, of the stress tensor. The components of stress, 
NX,, Y,, Z,, on an element of surface whose normal 
has the direction ¥ are given in terms of the stress 
tensor by the relations [8, p. 40] 


ds Trr COS(L,v) + Try COS(Y,v) + Tzz COS(Z,v) 


Y,=r,, cos(r,v)+17,, cos(y,v)+1,. cos(z,r) 7) 


Zy=Tez COS(Z,v) + 7,2 COS(Y,v) +722 COS(Z,r). 


The relations (7) may be regarded as defining the 
six components, Tzz, Tyy, -, Try, Of the stress tensor. 
If we let the vector ¥ take on the directions of the 
three coordinate axes in eq (7), we obtain the usual 
definition of the stress tensor, 


NX,=Tyzz, Y,=Try, Z.=T2 
X, = Try) ) y=Tyy LZ, = Ty 5 
is = Tory ) = T yz Z, = Te 
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From considerations of static equilibrium it can be 
hown that 


(9) 


nd that the defining relations (8) are self-consistent. 
For any material in static equilibrium under a 
rescribed system of surface loading .X,, Y,, Z,, and 
ee of body forces, the stress tensor must satisfy the 
quilibrium equations 
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n the interior of the medium, and the conditions (7) 
on the surface, where we now regard X,, Y,, Z, as the 
prescribed surface tractions. 

The strain tensor, €,,, €yy, €z2, € €rr, Cry, IS defined 
in terms of u, rv, and w, the x, y, and 2 components 
of the displacement vector, by the relations 
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It may be shown [5, p. 48] that a necessary and 
sufficient condition that a given tensor may be con- 
sidered as a strain tensor arising from actual physical 
displacements in an elastic material is that the 
tensor satisfy the conditions 
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Equations (12) are known as the compatibility 
equations and are the integrability conditions for 
the system of (11), when the displacements u, v, w 
are to be determined from a given strain tensor. 

The stress and strain tensors in an elastic medium 
are related by the generalized Hooke’s law. If the 
elastic constants of the medium are denoted by 
ey (ij=1,2, . . . ,6) Hooke’s law for a medium in 
which the zy-plane is a plane of elastic symmetry 4 





may be written in the form [8, p. 62]: 
€rr Curr CiaT yy Ci3T 22 CVeT ry 7 
C yy =Car Tez tT C22T yy + C237 C26T ry 
¢ Cai T rr Cg2Ty O€33T2 Cagl y 
an (13) 
¢ C4aT yz T C457 
( C547 Cs:T 
‘ CeT CegoT Ce3T CeeTry J 
where Che Ca;, Cy Cay. Cry Cei, Cag Cao, Cre Ca2, 


C36 = Cozy Ca5— C5t- 


For isotropic materials, the number of independ- 
ent elastic constants reduces to two, and Hooke’s 
law has the simpler form 
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where # is Young’s modulus and @ is Poisson’s ratio 
for the material. 

We may now reformulate the first fundamental 
boundary-value problem of elasticity in mathemat- 
ical terms: 

Determine a set of six functions, the components 
of the stress tensor, which satisfies the equations (10) 
in the interior of the elastic medium and the relations 
(7) on the surface of the medium, where X,, Y,, and 
Z, are prescribed functions, and such that the strain 
tensor calculated from these functions by the use of 
(13) satisfies the compatibility conditions (12). 


3. Formulation of the problem 


Let the coordinate system be chosen as shown in 
figure 1, with the built-in end of the cylinder in the 


‘ The most general elastic medium is characterized by 21 independent elastic 
constants. Here we consider a medium whose elastic properties are identical in 
any two directions symmetric with respect to the ry-plane. This type of elastic 
symmetry reduces the number of independent elastic constants to 13 





plane z=0, and the free end in the plane z=A. The | 


boundary curve of any cross section will be denoted | 
by C; the region enclosed by C and the unit normal | 


to C will be denoted by FR and ». 














Figure 1. 


The prescribed surface traction 7 (force per unit 
area) is assumed to be independent of 2, directed 
parallel to the (x,y) plane, and such that the result- 
ant on any cross section is a twisting moment of 
magnitude M (torque per unit length). 

This condition requires that the components 
and Y of the applied surface traction satisfy the | 
relations 


| Xas 0, | Yds—0, 
Cc JC 


In order that the first fundamental boundary 
value problem of elasticity may be formulated in 
such a way that the solution is unique, it is necessary 
that the stresses be prescribed over the entire sur- 
face, subject, of course, to the conditions of static 
equilibrium. This would require that Z, and Z,, the 
components of the clamping stresses on the end 
z=0, be prescribed, subject to the conditions 


(2¥ yX ds M. 


(14) 


(1Z,—yZ,dS=Mh. 
KR 
(15) | 


| 4d8=0, | as 0. 
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In this case the boundary conditions on the end | 
sections would be of the form: 





(a) on the end z=0: 7 | 

T2230, Tes=Zr, T yz sy: | 

- (16) | 

(b) on the end z=A: 
T,.2=0, T2r=0, Te=0. J 


Inasmuch as Z, and Z, are unknown (and depend 
upon the mode of clamping), we shall replace the 
exact boundary conditions (16) by the relaxed 
conditions, 








(a) on the end z=0: 


| r,.dS | rr.dS | y7..dS=0 
JR R JR 


| 72248 r,,dS =0 
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(b) on the end z=h: 





Tye 0. J 


These are conditions on the resultant stresses, rathe: 
than on the stress distributions. 

The boundary conditions on the lateral surface 
that is, on the curve C, are 


Trr COS (7,v)+7,, COS (Y,v) =X 


Tr, COS (r,v)+7,, cos (y,v)= } 17 


Ter COS (7,v)+T,, COS (Y,v) = 0. 

Our problem, then, is to determine a system of 
Stresses Tyr, Tyy, . . ., Try Which satisfies the equilibrium 
eq (10), the relaxed boundary conditions (16a) and 
(17), and for which the strain system calculated by 
the use of Hooke’s law (13) satisfies the compatibility 
equation (12). 

The solution of the problem as thus formulated is, 
of course, not unique, but, according to the principle 
of Saint-Venant [8, p. 95], the stresses given by any 
two solutions will be essentially the same at points 


_ of the cylinder not too near the end sections. 


We shall make use of the Saint-Venant semi-inverse 
method of solution [8, p. 100]. This consists in 
making certain assumptions about the stress distri- 
bution, and still leaving sufficient freedom to permit 
the satisfaction of the differential equations and the 
boundary conditions. If these assumptions are 
carefully made, a unique solution (unique under 
these assumptions) exists. 

It is apparent that the torsional couple trans- 
mitted across each cross section is a linear function 
of the distance from that cross section to the free 
end of the cylinder. Hence, it is natural to assume 
that the stresses in the cylinder are given, in part at 
least, by a Saint-Venant torsion function. 

There is a type of elastic deformation known as 
plane strain [7, p. 231], in which the displacements 
are of the form 
w=0. 


u=u(z,y), v=v(z,y), 
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lhe stresses in this type of deformation are derivable 
rom an Airy stress function and can be chosen to 
satisfy very general boundary conditions of the type 
riven by (17). 

We shall assume that the stress system in the 
wisted eylinder may be derived, in part, from a 
uperposition of the torsion solution and the plane 
strain type solution. 

For a deformation to be of plane strain type, it is 
wecessary that 7,,.40. Because a variable tension 
on the end sections may produce a bending or exten- 
sion of the evlinder, we shall further superpose the 
bending and tension type solutions. 

Thus, if we denote by 7, any one of the six com- 
ponents of the stress tensor, we may write 


Tr, Tis Thi 9 (18) 


where @ is a constant to be determined by conditions 
of static equilibrium, 7;! is derived from an Airy 
stress function, 7;? is derived from the Saint-Venant 
torsion function, and 7;* is required to satisfy the 
condition that the end surfaces be free of stresses 
that would produce bending moments and extensions 
The stress system 7,! may be written [7, p. 231] 
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where ¥[=¥(z,y)] is the Airy stress function and p, 
q, and r, are functions of the elastic constants. We 
note that the stress system 7;}) satisfies the equilib- 
rium equation (10) identically. 

The stress system 7,? has the form [8, p. 101, 108] 


te ~A+Br+Cy, te = TO = 19 =r P= OM =), 
(20) 


where A, B, C are determined by the three relations 


f, t,dS=0, _ [_ rr,.dS=0, 


f y7t,,dS=0, 
R 
(21) 


which state that the resultant bending couples and 
tension vanish on the end sections. 

_ The system 7{? is a modification of the Saint- 
Venant stress system 7/,. Sokolnikoff [8, p. 214] 
gives the stresses 7/, in an orthotropic beam twisted 

by end couples as | a 
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where 1, m, n are functions of the elastic constants 
and @|=¢(2,y)] is the torsion function. Because 
we wish the torsional couple to vanish on the face 
2>=h, we shall set 
: = Og ] 
Po (h =) (ay? )+m(5"—y )¢ 
(23) 
2 _\4,(9°° O¢, \l 
1 (h— 2) U( Gry )+m( site )- 
The third of the equilibrium equations (10) will be 
satisfied if @(z,y) satisfies the condition 
0’o 0’o 0’ 


l dr? 2m + 
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In order that the first two equilibrium equations 
in the set (10) be satisfied identically, we set 


) Mm «a , m ' 
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To simplify the satisfaction of the compatibility 
conditions (12), we choose 


' 8, ft. 
Tr = —OT5r—5Y, (26) 


with s and ¢ as vet undetermined. 
Our assumed stress system, then, is the following: 
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with the constants a, /, 
undetermined. 

We next determine p, 4, and ¢, 
that ¢,, will be independent of ¥(z,y) and of the 
terms in 2 and 7°. The insertion of the stresses 
given by (27) into the third of (13) leads to the 
following expression for ¢,, 


m,n, p,q, 7, 8, 
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the last five terms in the right-hand member of (28) 
will vanish. 
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If the values for p, q, 7 and ¢ given in (29) are 
substituted in the stress equations and the 
strains are then calculated from Hooke’'s law (13), 
we have the following expressions for the components 
of the strain tensor: 


(27), 
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If the components of the train tensor given by 
are inserted in the compatibility equations (12 


| is seen that all but the fourth of the compatibility 


equations are satisfied by the assumed stresses if /, 
m, and n can be chosen to satisfy the six equations 
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It is easily seen that only three of these six equa- 
tions are linearly independent and that the compati- 
bility equations are satisfied if we set 


l C33€ 44/8, m C33€s5/ B, N= C33€35/ B, 


where 


B= C44 (C55 + C31) + C55 (Cas + Ca2) — Cap (245 + Cye). 


For convenience in calculation we rewrite the ex- 
pressions for s and ¢ given by (29) in terms of the 
elastic constants, using the relations (33): 


$= (€31€45 — Cae€s5)/ B, t= (C3245 — Cae ag)/B. (34 


The fourth of the compatibility equations (12), 
the only one that is — on ¥(z,y), requires 


that ¥(z,y) satisfy the relation 
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In the next section we shall consider the boundary 





(Cs4n T C55m (Se+z)) 


conditions on the functions ¢(z,y) and y(z,y). 
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4. Boundary conditions 


It has been shown in the previous section that the 
.ystem of stresses given by 
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will satisfy the equilibrium equations (10) if $(z,y) 
satisfies 
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0, (37) 
and will satisfy the compatibility conditions (12), 
provided that 
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(38) 
C45(2C45 + Cag), 


B= C4q( C55 4 Cai) 4 C55 (Cag + Ca2 


and ¥(z,y) satisfies the differential equation 


‘ 4 4 
b Ov _ 35, o'v + 2( bi. + Deg) oy — 
or ~ Or Oy : Or*oy" 
; o'y oy Ooo o*o 
35,, b 6. —— 25, 
6 Szdyst " Oyt 2 Oz? 8 Ody > 
0’ ee 
b, dy? t Po» (39) 
where 
bi j=b 4 ci i, j=1,2,6, 
and 
bh, =Cylt+ceon—e Cam, i=1,2,3,6 








Since both of the functions ¢(z,y) and y(z,y) 
appear in (39), it is apparent that we must first 
determine ¢(z,y). Fortunately the third boundary 
condition of the set (17) is independent of y(z,y), 
so that the two boundary value problems are 
separable. 

If the values of r,, and r,, given by (36) are 
inserted into the third of (17), we see that this 
boundary condition will be satisfied if 


y! SF tm set cos (x,2)+}m aaa Se cos (y,) 
(ly —mx) cos (xr,v)+(my—nar)cos(y,v) (41) 
on ©. It is shown by Sokolnikoff [8, p. 215] that the 


change of variables 


> m 
f= yln—m’r, n=yY- ] I 

"(et 9)—vin—m? 2 
9 (5,7 j o(t,n) (42) 
will reduce the systems (37) and (41) to the simpler 

| system 

0°¢’(E,n) , O°@’(E, n) 
FAS) 8) 4 OS NS 99 (43) 
of on 
d 





jp? (&) n cos (£,v)— £ cos (n,v)=0 on C”. 
av 


The problem of determining ¢(7,y) is thus reducible 
to a problem of Neumann. Since the right hand 
member of the second of equations (43) is a perfect 
differential, the determination of $’(£,n) is possible, 
and is unique to within an additive constant. 

It should be pointed out that the quantity 
(In—m?*) is positive for all elastic media, since from 
the defining relations (33), we observe that 


» Cy - 
lIn—m?= BF (€C44€55— Cis) (44) 


and the positiveness of (¢y¢4;—cis) is guaranteed by 
the positive definite nature of the strain energy form.’ 
Once the function ¢(z,y) has been determined the 
right-hand member of (39) is a known function. 

Let us express the surface tractions Y and Y in 
terms of the tangential traction 7 and the normal 
traction \, by the relations 


X= N cos (z,v)—T cos (y,v), 


Y=WN cos (y,v)+T cos (z,v). (45 


5 See appendix 





stresses 


If we make use of the relations (45), the 
given by (36), and the relations 


dy dz 
cos (r,v)=-> cos (y¥,v)=->” 
7 ds 


ds 


the first two equations of (17) may be rewritten in 


the form 


4 (Y)_(ig4™ 2 


N\ 
dy - ) cos (x, t 


a 


-mor Jeos (y, »), 


a 


: , (47) 
-) (Re Sy — a )cos zt 


d (ow 
~ de ( Ou 





n A 
( no—= y’—— )eos (y, »), 
o~ a J 


on (. Let us denote the right-hand members of 
these equations by F,(z,y) and F\(z,y), respectively. 
If we can find a pastiouler (single -valued) solution, 
¥*(z,y) of (39), we must then solve the homogeneous 
system: 


o*y 
2( b10-+ Das) . — 
+ 2(O12-+ 06 dzty? 


3b, 2° ory 


3 Oi. droyit OM dy 


oy 
Or’ Oy 


b Ov 


25,1 3 bre 


ra(34) -- (2 oy ~)4+ F(z, y) 


Oy d 
-# (2) ds (x 


It is easily deduced from the results given by 
Sokolnikof [7, sections 66 and 95] in his discussion of 
the theory of anisotropic plates that the system (48) 
yey an essentially unique, single-valued solution 
¥ (z, y), provided that (a) the roots of the character- 
istic equation 





) + Fiz, y). J 


boo" 3bor° T 2 (bi. t Deg)” _ Bbi—d T by, =(0) (49) 


are all complex, and (b) 


| F.ds=0, | Fide=( 
Cc ¢ 


It was shown by Lechnitzky* [3] and [4], through 
energy considerations that the roots of eq (49) are 
always complex. Thus condition (a) is fulfilled. 

We shall test the functions F(z, y) and F;(z, y), 
in the first two equations of (50). These equations 
become 


| (xF',—yF,ds=0 
c 
(50) 


* Lechnitzky considered a similar equation, arising in the ocmty of anisotropic 
plates, whose roots are the negatives of the roots of eq (49 


| J (o+3 


(46) 


| we may rewrite the expressions (53 


| Salar 





. i ta 5 
z ) dy —( > in y- 


1¢ , an 
= | N cos (x, v)—T cos (y, v) } ds, 
JC 


a 


| (= a 4 
r*—— y* 
JC 2 2 Y 


Re os mn 
N cos (y, »)4 T cos (x, v)ds. 
aJwc 4 


T mo) dx 


-m@ ) dy —( no— = y? ) dx 





The right-hand members, by the use of eq (45) are 
seen to represent the components of the resultant 
force per unit length. 

By making use of Green’s theorem, 


e “(ON OM 
Ne 4 = ( ; - 
|, ro 208s LG& r OY 


we rewrite the left-hand members of (51) in the forms 
" §, (0@ 0d, .\? ac 
opel i My ’ 
Jim f_y)+ m (3° 


Since ¢(z,y) satisfies the differential equation (37 


) ds, 


Oo oo oo 
l -+2 0, 
drt m he dy? 


| tei he Site 
Oo _y)4-2m (32+2)t+ 
2 {em (2¢—y)+-2n (3842)} ]as 

54) 

f. [= i! yn (So+2)+um (22-»)}+ (54 


2 ym (S$+z)4 yt (5? -. yt] ds. 


Once more we use Green’s theorem (52), this time on 
the expressions (54). These now have the forms 


' = ae :, 4 
foz[¢ 5rt™ = ly ma} cos (x, v)4 


jm aaa 38 my+nz} cos (y,») ds j 


(55) 
Joy [i OF +m Sey +mz} cos (x, v)+ 


jm Sin $e my+nz} cos (2,») | ds. 
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The expressions in brackets vanish, by (41), and 
hence the first two of our integrability equation 
50) become, in this case: 


Yds 0 


( 


NVds=0, 


a ( 


These are the conditions which the prescribed 
omponents YX and Y of the applied stress were 
issumed to satisfy in section 3, hence the first two 
ntegrability conditions (50) are satisfied. 

Let us now consider the third equation of (50). 
We express N and T in terms of X and Y and see 
that we must have 


| (rY¥—yNX)ds 


Pee s l ; m a2lg 
el? J I~5 Yormyo rigor 5 ry (¢ J 


in l ; ] m 
rye? — 5 ty rmrge—y e—- 5 


ya? dy. (56) 


We apply Green’s theorem to the right-hand member 
of (56), to obtain 


F : oO re) 
|, [ nz?—2m, Y ly? + di ( m sete >y) 


O¢ , Og . nin 
‘ y (ls +m 3°) dS. (oa) 


Now, the left-hand member is proportional to the 
twisting moment per unit length, M, that is, 


Let us denote by D’ the expression 7 


, 4 = , , ; ~ ne 
D =|, [mz —2mrytly’+a ( 7 oat oy a 


o¢ de oe 
y(J >t ™ 32) Jas. (59) | 


Then, if @ is taken as 
a= M/D’, 


the third equation of (50) is satisfied, and a single- 
valued solution of (48) exists. 

We have satisfied the boundary conditions on the 
lateral surface; next we shall consider the relaxed 
boundary conditions on the end sections. 

On both end sections we require that there be no 
resultant forces. Hence we require that 


_' The quantity D’ is proportional to the Saint-Venant coefficient of torsional 
rigidity D, which arises in the theory of torsion by end couples. 


| constants 7, m, 


(60) | 





( T4dS 0. 
R 


The first of these conditions can be satisfied by 
adjusting the parameters A, B, and C. The second 
and third are easily shown to hold by considering 
(53) through (55). 

On both ends we require that the bending couples, 


M, and M,, vanish. These requirements can be 
met by adjusting the parameters A, B and C, so that 


M, — ( rr.dS 0. (O61) 
JR 


On the end 0, we require that M,—A-A and on 


the end z=hA, we must have \/,=0. The second 
condition is easily seen to be satisfied[since + 
Ty2=0 for z=h. 


Since M, for -=0 is given by 


M. hea [_ I [ ( > +7 ) +m ( s -¥ )] . 


y [: ( 5H y +m (5° +2 ) Jas, (62) 


and the integrand is equal to the integrand of (57), 
it follows that 
M.=haD’=hM. (63) 
Thus the stresses given by (36) can be made to 
satisfy the boundary conditions (16a) and (17), the 
conditions of static equilibrium (10) and the com- 
patibility conditions (12) by proper choice of the 
and n, and the functions $(7,y) 
and y¥(r,y). 
In the next section we shall specialize the equations 
to the isotropic case, and in section 6 we shall 
illustrate the use of the method by calculating the 


| stresses in an elliptic cylinder twisted by constant 


tangential traction. 


5. Isotropic beam 


If the beam is isotropic the differential equations 
and boundary conditions are considerably simplified. 
In this case the elastic constants c,;; may be expressed 
in terms of F, the Young’s modulus, and o, the 
Poisson ratio, as follows. 


Cy = Co2=— C3=—FZ 


; 
k 
y 


o 
Cog = Cag Cq1 = Cp. = C713— Cy = — FH 


; 
kf 
4 


Cu C55 Cen 





C16 = Co = C2g= Con Cag= Cog = C45 = Cy = O. 
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The elastic constants /, m, n, 6,, 6, are given by 


j=a_=} ) | 
m0 
oe ey, | 
beg = (1 +-0)/F | 
. (65) | 


bi. buy it | 


6, =6, h l+o)/E 


6,.=0. ! 


The function ¢(z,y) must satisfy the conditions 





Vols, y 0 ) 
de . (66) 
= y cos (4,v)—s cos (y, v), on ¢ j 
and ¥(z,y) must satisfy 
Viv(r,y)=0 5 


f= ($—* Joss (| += eos ty 


—dalSE Aa aeons H(G—Zooeino | 


(67) 





The torsional rigidity 1)’ is calculated from the 
defining equation 


,1f 2 2, , Ov Oo . 0) | 
D 9 Jf ty? ta dy! oe | aedy (6S) 


and is related to the Saint-Venant coefficient of 
torsional rigidity D by the relation 


D’=D(i+oyV/E. (69) | 


The stress system in this case has the simpler form, 


l l+o ote Oy 1l—o’ Oy ao _ 
a 9E°" EF dz* E oy? EA Be+Cy 

l l+o l—o’ Oy o+ec’ Oy o 
a” 9E*— EF dztt EF dy EAT Bt ey 
ve, TP Oty (A+ Br+Cy) 

1 ite, (da, f 
a’ oF (h- (a, +2) 

1 ite, (da 

l l+oa Oy 4! tors, 
la” E dyor' 4E™  * 


1» », 
hy ta YT Sry’ 


and the strains are given by 





6. Stresses in an elliptic cylinder 


As an illustration of the theory developed in the 
preceding sections, we shall calculate the stresses in 
an anisotropic cylinder of elliptic cross section 
twisted by a constant tangential traction 7. 

Let the equation of the elliptic boundary, C, of the 
cross section of the cylinder be given in the paramet- 
ric form 


r=a cos 6, y=b sin @. (72 
Then, 
cos (4, ¥):cos (y,”)=6 cos 6:4 sin @. (73 
Let us set 
or, y)= >) Ayr'y’, +, 7=0, 1,2. (74) 
| and 
Vz,y)= >> Byzr'y’, i,j=0,1,2,3,4, (75 
inc 


l oo OY ) 


T : 5 
a a oy" 


dy 


T > 
a ” 2 or 


: Tez —<aV*y—o L(4+Br- Cy) 
a 
. (70) 
l i — Og 
a “- 3 (h =)(5,+7) 














and seek to determine the coefficients A, and B,, in 
such a way that (37), (39), (41) and (47) are satisfied 
We first calculate \/ from the defining relation (58), 


M=| (x¥ —yX)ds. (76 
a 


If we insert in (76) the expressions for Y and X given 
by (45), we find that M is given by * 


* Although N is assumed to be zero in this discussion, we include it in eq 7 
to show that a constant normal traction produces no twisting moment. 











ir 
in 
mn 








V/ | (—yN t rT) cos (x, v) t (yT t rN) cos (y,v) ds. 
oe ( 


(77) 


By making use of Green’s theorem (52), we may re- | 


write (77) in the form 


M=2T | dS=2rabT. (78) | 
JR 


We insert the expression (74) into (37) and (41) 
and, after some manipulation obtain 


1 , , . , 
ol, y)=-s5 = —bh?mr? t (b*l—a*n)ry ramy 
. a*n- 6? 7 ' 


(79) 


as the solution of the systems (37) and (41). 
Using the value of ¢(z, y) given by (79) in the defin- 
ing equation (59) we find that D’ is given by 
ra*bh*(nl—m?*) 


Dp’ é . . (SO) 
(b?/ an) 


Henee. from (60), it follows that 


2(b7/4-a?n)T 


7 = © (Sl) 
a*bh? (nl— m*) 


a 


We next insert the expression (75) for ¥(z,y) into 
the two equations (47) and after setting r=a cos @, 





y=6 sin @, and making use of (73), we obtain a pair 
of polynomials in sin @ and cos 6. These are expanded 

| in Fourier series and the coefficients of the various 
terms are equated to zero. We thus find that, in 
order to satisfy the boundary conditions, y(r,y) 
must be given by the expression 


24a*b* (a?n + bl) -W(z,y) = 6' Blmb? — 5Smna? + y)2*- 
4a*b* (b?m? —a?n®) xy + 20°*b* yy 


+ 40°? (6°? —a?m?) ry’ + at (Slmb?— 3mna?+ y)y' 
2a°b* (3lmb?—3mna?+- y)2° 
2a*h? (3lmb?—3mna?+y)y?, (82 
where y is a constant, as vet arbitrary. 

If the value of ¥(z,y) given by (82) is inserted in 
the differential equation (39), it is found that y must 
satisfy the relation 

3b.ob* +- 2075? (by + bya) +- 3,04 } +7 

3b..b*(Smna?— 3lmb?) + 9b.,a°b? (b?m? — a?n*) + 

9b,,a7b? (6°? — a?m?) +-36,,a4(3mna*— dlmb?) — 

6b,a7b'm + 66,076? (a2n — 671) + 6b,a*b?m +- 
30°76" (a?n + 671) (bogn —bygl). = (83 


Inserting the values of (x,y), a, and ¥(z,y) given 
above into (36) we obtain the stress system 








ter (3mna?—3lmb?—y)(b?2?4+-307y?—a*b?)—6a7*b*(lIn—m*)ry ) 

- 3a?b4(/n—m?*) 

tr, (3mna?—3lmb?—y) (3672? + a?y?—a*b*) + 60° b*(In—m*)ry 

7 3a*h*(In—m?) 

ty 3ln—m?*)(b?2?—a*®y*?)—2(3mna*?—3lm b?—y)ry 

7 3a? b*(In—m?*) 

> (S4 

Ty  42(h—z) 

f- a’ 

Tez 4y(h—z) 

T b? 

arp’ r?+q’rytr'y?+A4 Br+Cy 

1 J 
where 

, 367e3(3lm b? —5mna?+ y)—3a%eq¢( b?m?—a*n*?)+aeny ¢ym- Caglt b?m > 

} - a — > > 
' 6a*es9( 6°71 + a?n) C33 (b*/+-a?n) 

»  3x( b?m?— a?n?) —cagy +3¢y( 672? —a?m?) + 3¢,;(a?n — bl) : (ns 
. Ta ) 
. 3C39( b* + a*n) 
wins Bey — 3 b7e56( 670? — a*m?*)+-3a7¢4,(5lm b?—3mna*4 7) 4 Cum ; Cpl a*m : 

6 b7e3,( 671+ a?n) C33 b*l+-a?n : 
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To evaluate the parameters A, B, C, 
that the conditions 


in order 


r,dS=0, rr..dS=0, y7,,dS=0 
R JR ? 


Jf 


. 
are satisfied we write 


. 5 


(px? +q'ry+r’y+At+Br~Cydrdy=0 
JR 


(pr +q' rytr'rty t+ Ar+ Br? +Cruy)drdy=0 > 
R 


(p’xtytq' ry tr’y 
JR 





Ay+ Bry + y*)dady 0. | 


(S86) 


If we denote by R, J,, and J, the area and the 

moments of inertia of the elliptic section, the system 
86) can be written in the form 

p'l,+r'I,+AR=0, BI,=—0, CI,=0. (87) 


For the ellipse, 


rh a ra 4 
R rab, I ’ . (S8) 
1 t 
Hence the parameters A, B, and C are given by 
a*p’ + hb?’ . 
A=- u » B=0, C 0, (89) 


4 


and the proposed problem has been completely solved. 
For the isotropic elliptic cylinder the solution may 

be simplified by using the relations given in section 5. 
The solution for this case is given by: 





b?—a . 
Or, Y) ; , ry 
P a’?+- 5? “* 
rb® 
Dass 
2(a?+- b*) 
_— 4(a?+-b*)T 
a’ bh’ 
y=0 
b?ry'—a'*r®y 
x,y): ©, = 
uz, y) 12(a*+ 6?) 
Tez 2ry 
es 
> (90) 
Tw_ 2ty 
zz ¢€ 
Try _ b*2z?@—a*y?| 
r ees: 





ty. 42(h—2z) 

T a- 

tr _4y(h—2) 
Zz b? Ci 
tz: 20(a*—b*)ry 

T a*b? J 


For an isotropic circular cylinder of radius a 
results are given by 





dor, y)=—0 7 
a 
pb’ —* 
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ST 
a " 
a~ 
ry°—ax°y 
vr, y) 24 
Tos ry 
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Ton 2ry | 91) 
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in “-y 
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Ty: 42(h—2z) 
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7. Methods for obtaining solutions of the 
boundary-value problems 


The determination of the stresses and displace- 

ments in the twisted cylinder has been made to 

| depend upon the solution of the two linear partial 
| differential equations, 





o*¢ ‘ oo 0*d = (Q” 
l 5x2 t 2m Sady** dy?» (92) 
and 
ory oy ory oy 
bx rt 32 dxidy t 2012 + bee) dz?dy? 8 dzdy'* 
ot" 0*¢ 0*¢ 0*o ” 
bi; Sy 5x2 246 doy to dye t Om (93) 
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ogether with the appropriate boundary conditions, | 


t1) and (47). 

It was shown in section 4 that (92) can, by a simple 

hange of variables, be transformed into the form 

Vo(r,y) 0, (94) 
nd the various methods of potential theory are avail- 
ible for the solution of this problem of Neumann. 
in elegant method of solution of the Neumann 
woblem, due to Muschelisvili, is discussed in 
Sokolnikoff [8, p. 170]. 

The determination of the function ¥(z,y) is con- 
siderably more difficult. An equation similar to (93) 
arises in the theory of anisotropic plates and a great 
deal of work has been done on the problem of solving 
this equation. A series of papers by Green and 
lavlor [2] gives a method of solution of the system 
93)-(47) for the special case (corresponding to the 
orthotropic beam), where 

b6..=6,,.=—0, (95) 

A perturbation method of solving the system (93)- 
17), which consists in expressing the solution as a 
series of biharmonie functions, is described by 
Sokolnikoff [9]. Another, more elegant method, 
essentially due to MuscheliSvili, depending on con- 
formal mapping is also described by Sokolnikoff 
7, p. 375] and [10]. 

If the beam is isotropic, (93) has the simple form 


Viviz,y) =—0, (96) 


and the solution is easily obtained in the form 


¥(zr,y) =R[f(z) + Zg(2)], (97) 
where f(z) and g(2) are analytic functions of the com- 
plex variable z{[=2-+-iy], which are chosen so that the 
function ¥(z,y) given by eq (97) will satisfy the 


boundary conditions (47) [10]. 
8. Summary 


The determination of the state of stress and the 
displacements in a built-in anisotropic cylinder 
twisted by forces uniformly applied along the lateral 
surface has been made to depend on the solution of 
two boundary-value problems. 

The existence of the solutions of the boundary- 
value problems has been demonstrated, and several 
methods of solution of the boundary-value problems 
are indicated. 

As an illustration of the method, the stresses in an 
anistropic beam of elliptic cross section twisted by 
constant tangential traction were obtained by assum- 
ing that the stress functions were second- and fourth- 
degree polynomials. 


In a succeeding paper it is planned to extend the 
theory to the case where the applied surface tractions 
vary along the length of the cylinder, that is, where 


the boundary conditions are given in the form.’ 


NX =f(z).X(s 
g(z)¥(s) on C@. 
Z=0. 

It is hoped that the theory can be further extended 
to the problems of extension, bending, and flexure of 
anisotropic cylindrical beams by forces distributed 
on the lateral surface.” 


The boundary conditions for these problems are: 
(a) Extension: 


0 on’ 
0 onc, 


and Z is prescribed on C subject to the conditions 


( Zds+0, rZds=0. 
JC J ¢ 


| yZds=0. (99 
JC 


(b) Bending: 
0 ond 
“0 onC, 


and Z is prescribed on C subject to the conditions 


| Zds=0, rZds #0, [ uZds 0. (100) 
Cc JC 


Cc 


(c) Flexure: 


7=0 on, 
and VY and Y are prescribed on C subject to the 
conditions 


| Xds <0, | Yds+0, | (2Y 
Cc Cc Ce 


It is apparent that a superposition of solutions 
corresponding to the boundary conditions (98) 
through (101) can be made to satisfy quite general 
boundary conditions on the lateral surface. 


—yNX)ds=0. 
(101 


* Zwolinsky and Riz [12] treated the case (for isotropic cylinders) where / 
N=0, T=constant on C. 

'* The problems of extension, bending, and flexure of isotropic (and orthotropic 
cylinders by forces applied on the end sections have been solved by Saint-Venant 


[8 p. 100). 
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10. Appendix. Strain Energy and 
Elastic Symmetry 


The strain energy, W, of a strained elastic medium is 


defined by [8, p. 88] 


2 Text F aCe y of T sylyz- 102 


re‘ee 


linear 
to be 


In the theory of elasticity the strain energy is 
assumed given also by a positive definite quadratic 
form in the components of the strain tensor by 
W = Ibje?, + by00 


: 2 rz€ yy 


7 ! b, sf 


T Dis aet an” DisO esl ys T Dis€r2€ ex T Dist relzy 


. I + + 2 . p > 
Shove vytee bos yylve bos yyl ez bot yylry 


1 2 
’ bb are bag se€ye bast s2€ ex 4 base arlzy 


1 
Mayes, bes ysl ex t bys ysl zy 


, . 
+ bbsse2, + bse sexy 





i 2 
; beet * 


(103) 


rhe positive definite nature of the form (103) imposes certain 
requirements on the signs of the elements of the matrix 


k » te be & ; &* 
he Os Os fy fe be 
b, 2 ; 5 b 6 


bi 











An elastic medium is said to possess elastic symmetry j 
plane if a reflection of coordinates in that plane leaves + 
strain energy form (103) invariant. Let us consider elas 
symmetry in the zy-plane. The form (103) must be invaria 
under the transformation 


The components of the 
to the relations 


and if the form (103) is to be invariant 
tion, we must have 
Bg = 015 = bog = 25 = bas = Dig 5 , 107 


The matrix (104) then has the form 








KN big b . b ry 


Since the matrix (104a) is obtained from a positive definit: 
quadratic form it follows that 


Bas - bes! 
(10s) 
es te 


By differentiating (102) and (103) with respect to e,., and 
comparing, Hooke’s law is obtained, expressing the stresses 
in terms of the strains 

If this expression of Hooke’s law is solved for the strains 
in terms of stresses and the coefficients are denoted by 


Cu, « « «» Cee &8 in section 2, it will be seen that the determinant 


! ' 
Cu (C45 

(10Sa 
% C5 


is a positive multiple of the determinant in (108), and hence 
is positive, as asserted in section 4. 


Los ANGELES, May 26, 1951. 
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On Methods for Obtaining Solutions of Fixed End-Point 
Problems in the Calculus of Variations’ 


Marvin L. Stein? 


Two methods for constructing solutions to the problem of minimizing an integral in a 


certain class of ares joining a pair of fixed points are proposed. 
a generalization of Newton’s method, while the other is a ‘“‘gradient”’ 


One of these procedures is 
method. Conditions 


for convergence to a strong relative minimum are given in both cases. 


l. Introduction 


In this paper we consider methods for effectively 
constructing solutions to the problem of minimizing 
an integral 


T(y) f(r,y,y dx (1) 


in a certain class of ares y;(x), (a<r<b; i=1,2, 
, nh) joining the two fixed points (r=a, y(a) 

o) (4=b, y(b)=e). As is well known, the minimum 

must necessarily satisfy the Euler equations 


d 


fAr.y.y’ =f, (r,y.y’) 
dy Y be vi Js 


Thus in any particular case the solution of these 
equations becomes of prime importance. The two 
methods to be proposed below enable us to compute 
ares that satisfy the Euler equations subject to 
boundary conditions at two distinct points. For 
both methods an initial estimate of the solution is 
made. From this estimate a new estimate is ob- 
tained in a definite way, and so forth. These pro- 
cedures have the advantage of giving estimates that 
satisfy the boundary conditions at each step, and 
their iterative nature makes them adaptable to auto- 
matic computation, 

The first method of obtaining new estimates that 
will be discussed is a generalization of Newton’s 
method for functions of a single real variable. This 
method is described in section 2 where it is actually 
applied to a more general class of systems of differen- 
tial equations than the Euler equations, namely 
those systems of the form 


d 


= gir, y,y =hdz, yy’) 


To construct the convergence proof, we make the 
class of functions under consideration into a Banach 
space and regard the system of differential equations 


rhe preparation jof this paper was sponsored (in part) by the Office of 
Naval Research. 

* Present address, Consolidated-Vultee Aircraft Corporation, San Diego, 
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as an operator 7’ on the space into itself. The neces- 
sary Banach space tools are developed in sections 3, 
4, and 5. Sections 6, 7, and 8 are devoted to an 
analysis of the properties of 7 and of its first varia- 
tion. The results obtained in these sections are used 
in section 9 to show that under simple conditions on 
the initial estimate the sequence of estimates derived 
from the first method will converge quadratically to 
an admissible function that satisfies the svstem of 
differential equations and the two-point boundary 
conditions. In fact, it is demonstrated in section 9 
that convergence will occur if the initial estimate has 
nonconjugate end points (see definition 5) and is con- 
tained in some neighborhood on which the norm of 7 
is bounded by a suitable constant. These conditions 
are shown to be satisfied if a solution exists and has 
nonconjugate end points. In this case there will 
exist a neighborhood of the solution such that the 
sequence of estimates will converge for any initial 
estimate chosen from this neighborhood. In_ the 
special case of the Euler equations, considered in 
section 11, it is shown that the limit function not only 
satisfies the differential equations but under the 
proper conditions also affords a strong relative mini- 
mum to J(y). Section 10 concerns itself with an 
interesting lemma about conjugate points, that 
proves useful in section 11. 

While not having as large a range of applicability 
as the first method, the second of the two methods of 
choosing new estimates that will be discussed has the 
advantage of being simpler to handle at each itera- 
tion. This method is described in section 12, where 
it is shown to be particularly useful for the solution 
of the Euler equations arising from an integral whose 
integrand is quadratic in y and y’. In this case the 
integral itself, instead of the differential equations, is 
considered; and it is shown that under suitable hy- 
pothesis the sequence of estimates converges to an 
are that affords a strong relative minimum to /(y). 
In fact, convergence will occur if J has a lower 
bound on the class of admissible functions and the 
second variation is positive. 

Neither of the methods to be described below is 
presented as a quick computational procedure for a 
desk computer. However, they are processes that 
might well be carried out on automatic computing 
machinery. Numerical examples are to be found 
in section 13. 








2. An Iterative Method for Solving Certain 
Systems of Second Order Differential 


Equations 


In the following pages we shall be concerned with 
constructing solutions for systems of the form 


d 
(gdz,y.y’))=hidz,y,y’) 
dg . . ‘ ‘ ‘ 


y¥(a)=p,, ylb)=a;. 

Here y= (%.%,. - Ya) denotes an are defined on the 
interval a<x<b and p, ¢ denote arbitrary constant 
vectors. By formal integration of (2) it is found 
that an are y that satisfies these equations must 
also satisfy 


T2439) [ ata’ hissyyds | at 
ot * 1 
aft? - 
| gly.) | his, y.y’)ds |dt=0 
ada! 4! a Y,! 


(3) 
y(a)=py,, yilb)=oe jae BD, . . op) 
Conversely, any are without corners that satisfies 
3) must also satisfy (2). For reasons that will be 
clarified below, we prefer to make use of the differ- 
ential eq (2) in their integral form (3). Accordingly, 
let us examine some of the properties of the function 
T (xy). 
That 
T (a; y) =7T,(b; y) =0 (j==1,2,...,2) (4) 
independent of the are y is immediately apparent. 
* One recalls that a function of one or more vari- 
ables is said to be of class C*, n>0, in its arguments 
on its domain of definition if it is continuous on that 
domain together with all of its (partial) derivatives 
up to and including those of order n. An are 
Y= (%,Yo,- . .Ya) is of class C* if all of its compo- 
nents, y;, have this property. From further inspec- 
tion of the definition of 7;(z;y) it is evident that if 
g and A are of class C*~' and y is of class C* for some 
k on O<k<n, then 7,(x;y), considered as a function 
of z, is also of class C* that is, 7(z;y) belongs to the 
same class of arcs as y, provided g and h have proper 
differentiability properties. Let a primed variable 
signify differentiation of that variable with respect 
to z. Then it will be noted that the are 7’’(z;y), 
which when set equal to zero merely presents the 
differential eq (2) in their usual form, does not pos- 
sess the property of belonging to the same class to 
which y belongs. In fact, even if y is of class C°, 
T’’ (x;y) may only belong to class C’. 
Before describing an iterative procedure for con- 
structing ares that satisfy (3), one must first introduce 


| y have sufficient diflerentiability. 


_ regardless of the choice of y and 7. 





the functions 


69 (2, 4; 0) =P, (2.4 Uns t Galt. YY Nn; 


bhdz,y;n=hy nythyan (t,9==1,2,... 


and 


t 
bh ds, Yn as | dt 


6T (zr, y;n)= | [ batt. ys 9) 
~e@ 


Ja 


r—a (° "t 
-- ; | [oy - | ahd | dt (gaz I, 2, 
b—aJal ° Je 


In (5), (6), and elsewhere in this paper, unless other- 


0) 


| wise noted, the tensor convention of summing on 
| repeated indices is utilized. 


The are 7 is taken to 

be such that it vanishes at the end points of a<r</ 
The function 67;(z,y;n) can be obtained from 

T (x;y) as 

TiKz;y- 


- 1 
67 (2,439) = A Q@N) a 
. da 


' and is often referred to as the “variation” of 7, 


It has properties similar to those of 7;. For example, 


| 67,, considered as a function of x belongs to the 


same class as the function 9, provided that g, A, and 
This fact follows 
from a small bit of reflection on (6) as does the fact 
that 

57 (a,y3n) =8T (by) =0 (i - ee a 
It is seen that 
these properties do not hold for the variation of the 


| differential equations (2) taken in their usual form. 


In the iteration process to be discussed here, it is 
assumed that the 7“ estimate y“ of an are satisfying 
(3) has been obtained. An “improved estimate” 
y“*» is then found as follows. Determine a ‘‘varia- 
tion” » by solving the system 


T(x3y) +67, (2,y® jn) =0 
nj(a) =n,(b) =0. 


The function y“*" is then chosen to be 
ys) y? ae 


y is replaced by y“*”, and the steps just described 
are repeated. 

This method may be justified heuristically by 
observing that upon expanding T7,(r;y+ 7) by 
Taylor’s theorem one gets 


T (x3y+n) =T(z5y) +6Ti (ayn) +. « 





sus it is reasonable to assume that a good approxi- | 
ition to a solution of 


T (xiy+n) =0 (omc) 7... sft) 
a known function y might be obtained by solving 


T (xzy) +67 (2,y3n) =0 G=13....2) 6 
yn. The analogy between this procedure and 
wton’s method for a function of a single real 
ariable is clear. Consequently, it will henceforth 
referred to as Newton’s me thod. 

\s one sees by setting g,(z,y,y’) =y; ((=1,2,.. .,n), 
he differential equations (2) contain the important 
lass of equations 

yi =hilayyy’) (i=1,2 (11) 


* 


as a special case. By setting 


glaysy’) =f, (r,y,y') 
hilryy’) =f, (ruy’) 


”) can also be made to assume the form 
d f 
dr on 
(13) 
Yl) =p, yi(b) =o. 
These are the Euler equations associated with the 
integral (1), and their solution is the chief applica- 
tion of Newton’s method which we wish to make in 
this paper. The functions 7, and 67, associated 


with (13) are, respectively, 


ry) ( [ys fds | at - 
a [ys 
5 Je i. 


-| wn ds | dt— 
“a, t 
: [« — [ Wy ds | dt 
~@ Ja ae ; 


“J 


and 


Tru; 9)= | [ 


where w,; 
form 
fy; ni 0, + Suzy; nj 


20( 2,459) Sow (2yy’) nes 4 


«ylt). 


(GQj=1,2,. . 
\ deeper insight into Newton’s method as applied 
to (13) may be acquired by 
theorem to the integrand f(x;y+7,y’+1’). 
expansion results in 


The 


F(ayy T 


| equivalent to (10). 


| are equivalent for y>0. 


and w,, are derivatives of the quadratic | 


applying Taylor’s | 


ny’ +a) —flauy)+Sy (uy at flew mn 


T + Avy mins +. T how; i 0; +4 hy, ¥; mn, T 


Sayy) thy yy thy + o(z,y,n) 4 


_ Let us drop the bigher order terms and consider the 


problem of minimizing 


6 
J(y+n) f, (f(z,y,y’) 


thyme tSyn + o(2,y;n) |dz, 
where y is a known function and 7 is variable. The 
Euler equations with respect to » of this integral are 


d d 
dz" Sy, +9 r 


wy, =0 
This last system of differential equations is formally 
Thus Newton’s method con- 
sists of repeatedly replacing f(z,y,y’) with an inte- 
grand that is quadratic and minimizing the new 
integral so obtained. Consequently, if the original 
integrand is quadratic, Newton's method should 
vield a solution in a single step independently of the 
choice of the initial estimate. That this actually 
occurs is readily verifiable. 

Finally, we note that although the Euler equations 


_ often have an equivalent formulation as a system 


of type (11), Newton’s method may not necessarily 
be the same for both forms of the problem. This is 
due to the fact that in general 67 is not invariant 
under a transformation of 7 to an equivalent form. 
For example, the two differential equations 


,2 


d 
—— (yy’) =1+-21 
dz YY Y >» 


In the first case 


I t : 
k[y*(x)— y*(a)| — ( [ f. (1+ 2y"\ds |at 
/@ /@ 
—;— « ) bly*(6)—ya)|— i [ fia L2y’ ‘ds Jar} ’ 


T(z,y)= 


| and 


z t 
57 (2x,y;n)= y(r)n(x)—4 ( lf y'n'ds |at 
42- a 
+f [. ’ y'n ‘ds Jat 


1 z *¢ 1 j y" 
T(x;y)= y(r)— y(a)— [ [| ( a )as Ja 
r—a ( ” ° ees 36 oy 
ag » ¥() ya)— fi [ f. = )as Jatt. 


| In the second case 





Ml (Zuri, Jas Ja 
/@ o 7 7] y 

] a bd / °, oy’ . 4 

h | lf. ( . o : / = )as Jar 


It is seen that the two variations are not the same. 
Since the conve rgvence theory to be deve loped below 
will be general enough to include either (11) or (13), 

one will be free to use whichever system of differe ntial 


equations seems to be simpler in any given problem 


3. Preliminaries 


There are several definitions and known results 
of general analysis to which there will be frequent 
occasion to refer It will be convenient to have 
these collected in one section Accordingly, the next 
few pages are devoted to this end. The reader who 
wishes more detail about what is to follow in this 
section is advised to consult a paper of Graves and 
Hildebrandt [1].2 The treatment he will find there 
is more general and couched in somewhat different 
terminology. 

Functions defined on a linear subspace l” of 
Banach space Y into itself are dealt with here. 
These functions may depend, in addition, upon a 
parameter chosen from Y or a region of Y. When 
the arguments are indicated, the parameter will 
always appear first and be denoted by a lower case 
y perhaps bearing a superscript. It will be set off 
from the independent variables, which will be denoted 
by Greek letters ", ¢, & etc., by a semicolon 

Dertnrrion |. F(y;m) is said to be continuous in 
y at y uniformly relative to » if given e>0 there 
exists 6>0 such that 


Fiy;n)— Fy ;n en (14) 


for all y in (y™)s and every » in (7. Here and else- 
where (y),; and similar notations denote the 6-neigh- 
borhood of y° 

The extension of definition 1 to the case of a func- 
tion F(y;n, ¢,. 
is clear. The only essential change is made in re- 
placing eq (14) by 


Fly; 9,¢,...,8)—Fly 


fre 


+ eer 


All of the definitions that follow may be generalized 
to functions of several variables in a similar way. 

Derinition 2. Let Y, be a region of Y. Fly; n) 
is said to be uniformly continuous in y on Yo uni- 
formly with re spect to n if for every €>0 there exists 
a 6 >0 such that 


Fey 39) — Fly® 50) elim 


Figures in brackets indicate the literature references at the end of this paper. 


., £) of several independent variables | 





that satisfy | y y |) <i 


for all y, y® in Y, 
all n in Ll’. 

Derinition 3. F(y;n) is said to be 
y" uniformly with respect to n if there exists a \ 
stant W>0 such that 


» hound: 


Fy sn) <M in 

for all nel ; If some \J>0 holds for all y belong 
to some neighborhood | eet F is said to be 
formly bounded on (y)q uniformly with respect t 

For the sake of brevity the words ‘“‘uniformly wit! 
respect to 7’, that appear in the above definitions 
and that should appear in the concluding definit 
(4) of this section, will henceforth be omitted. Thy 
reader, however, should supply them for himself in 
all of the necessary places. 

The following lemma is an eas\ 
the definitions. 

Lemma l. Jf F (ysn)tol 
and bounded at ae” there exists a ne ighborhood of 
such that F is uniformly bounded on this ne ighborhood 

The conclusion follows from the inequality 


conmseq uenhe { ol 


is continuous in y at ¥ 


Fiysn)) <) Foy sn F(ysn)— Fy sy 
DeFINITION 4. F(yzn) to U’ is said to be of class 
C”’ in y uniformly on Yo in case there exists a function 
dF (y;n.¢) on Y,UU to U having the following proper- 
ties: 

dF is uniformly continuous in y on Y, 
(2) dF is linear in ¢ and uniformly bounded on Y 
(3) the function R(y™ ;y® yy) defined by 


F(y™ jn) — Fly 3n) —dFy® jny® —y 


Rey sy jn) yO —y ||, yy 


is such that given e>0 there exists 6>0 such that 


Rey sy 39 €7 


for all y" .y® in Y, satisfying | y?—y™ <6 and all » 
in L’. 

Here and elsewhere 
space Y. 

LemMa 2. If F is of class Cin yun iformly on Yo, 
there exist non negative constants a and M such that for 
every y and y® of Yo for which \y—y™ | <aitw 
true that 


refers to the zero of the 


F(y™ :9)—Fly® :9)|| SM yy —y®)\\-\\n 


for all nin U. 


The conclusion of this lemma follows from an 
application of (2) and (3) of definition 4 to eq (15) 

Finally, we state a lemma that corresponds to 
lemma 16.3 of Graves [1]. 

Lemma 3. Let F (y;n) to U’ be bounded at y and 
linear in for each y of the region Yo. In addition let 
y” be a point of Y, and let F have a bounded inverse 











for 
the 


for 


at 
for 
is j 
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rer 
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| y’ sn) with respect to n and be of class C’ in y 
iniformly on Yy Then there exists a constant a such 
that for each y in (y 


(ony un iformly on (y . 


a I Yon) €2 ists and is of class 


4. A Convergence Theorem 


In section 2 an iterative method for constructing 
solutions of the Euler equations was described 
However, no mention was made of conditions of 
convergence. In this section we will describe and 
vive conditions for the convergence of a= similar 
procedure that can be applied to a large class of 
funetions. Later the process of section 2 will be 
considered as a special case of that to be given here, 
and it will be shown that under simple assumptions 
the special process satisfies the hypothesis of the 
convergence theorem of this section. 

Let S(y) be a continuous function not identically 
zero defined on a region ), of a Banach space Y to a 
subspace (7. The problem of imterest is that of 
obtaiming solutions of the equation 


S y ‘a 16 


Toward this end one introduces a linear operator 
Kiys) on L to UL, which depends in addition on 
some parameter taken from Y. It is assumed that 
Kiysy) is nonsingular in yn. Let the 7th estimate 
i of the solution to (16) be known. The 7+ Ist 
estimate is constructed as follows. A function 7 
is determined by solving 


S y K "se HW 17 


and the function 7] is chosen to be 


y"" yOtn®, IS 


The process is then repeated. 

In order to simplify the statement of the con- 
vergence theorem, several conditions will now be 
introduced on the function S(y) and on A™~'(y;n), the 
inverse with respect to » of A(y;n). 

Ll). The are y® of Yo, the constant y>0, and the 
operator K (yn) are such that there exists a positive 
constant M for which 


KO (yn) — KO (yn) <M ly y? |e! in 19 


for all n in U and all yy in ye dee Furthermore, 
there exists B>O such that 
K(y3n) <Bon (20 


for all y in (y™), and all yn in U. 


From (19) it is clear that A™! is continuous in y 


at y°. Therefore by lemma 1 one sees that in order | 


for (20) to hold on some neighborhood of y®, which 
is also a subneighborhood of (y),, it is sufficient to 
assume that A™'(y;n) is bounded at y®. The 
remaining part of (1) is, of course, undisturbed on a 
subneighborhood of (y), 

In the next two conditions certain non-negative 
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constants /7, L, and N appear. These are taken to 
be such that thev satisfy 


H—B-L~M-N<1. 21 


(Il ; For all 2 in y ' the inequality 


Ky? sy 


is true, 

This condition asserts that A is an approximate 
differential of S on the neighborhood (y a 

Ill). There exists a constant C such that 0< C-: c' 
and 
INC B, 23 
while 


for all other | in Y . 

THEORI Ml. Let the initial estimate 7] and the 
positive constant + be such that (1), (1), UIT) are 
valid, Then the se quence of ¢ stimates determined by 
(18) conve rges to a unique element y of (Yy a The 
are 7] is the unique solution of eq' 16) in (y° y: 

First we prove that 


y yt I , TT’ l Hl €%. y in y : ) >I, 


In order to construct the proof by induction con- 
sider the inequality 


y y y | K ~ aad S y ) y , 
<B Sly <(1—-MC<4 
This last follows from (23) and establishes (25) for 
i=1. Assume that (25) is true for all << m—1. 


By this assumption in conjunction with (19), (22), 
and (24) one may write 


yom —y" yi" +A (y("-? -— Sly 
yin Ky"? ;—S(y" 
y 7] Kk-'(y SY } 
A y Sly ko '(y 
S(y\™ Ko '(y Sty 
Sk“ KG; 
y"" y Ss(y™ Sy" 
+ M||y'*-? —y*-||-|| Sy 
SB|K(y* 5y"-? —y®) —Sy"*"") + Sy" 


+ MM) lyo"—) — y™ 





This establishes the first assertion of (25) for 7 m. 
The rest of the proof of (25) follows from the in- 
equality 


m m 
y'—y® ayy) Ss Qa lly —y" 
i=l i=l 


<1 C= (1-H) C< 9. 
i=l 
(26) 


If a sequence {y‘° } in a Banach space is such that 
the series Sy" —y""" | is convergent, it has a unique 
limit. Therefore, by (25) there exists a unique ele- 
ment 7 of Y that satisfies 


lim | y?—¥)) 0. (27) 


+a 


It follows from (26) and (27), that 7 must belong to 
(yy: 


Upon passing to the limit in the expression 
y m yi™ l K ‘iy m—1 :Sty m )), 
‘and S are continuous, 


and recalling the fact that A 
one obtains 


AK a; SQ)) ” 


This implies S(7/) =@. 

To see that 7 is the unique solution of (16) in 
(vy), let y* be another element of this neighbor- 
hood satisfying (16). Then by (20) and (22) 

'(y; S(y* 


¥—y*||=ll¥-y*-K'YG;SY))+K 


KUNG; AG —y*) -S@M+Sly* 
<B KG;y¥—y*)—SG+Sly* 
<BL\\j—y* 


Since (21) implies that BL<1,7 =y* 


It is possible to modify the method of obtaining 
iterated estimates by always holding y in K(y;n) 
fixed at the initial estimate y®. This is equivalent 
to originally introducing an operator K(n), which 
depends only on functions from (7. In this case, of 
course, the convergence theorem may be simplified. 
Condition (1) cle anty no longer applies, and condition 
(111) may be modified so that (24) is eliminated. 
The resulting theorem may be stated in the following 
way. 

Tueorem 2. Let K(n) be a nonsingular linear 
operator. Let K~'(n) be bounded and let the positive 
constant L be such that \|\K~'\|\L<1. Then if the 
initial estimate y and the positive constant a are such 
that 

(1) for all y™, yin Y™)a 

K(y™—y™) 


S(y) + Sy)! < Lily —y® 


the sequence of estimates {y } 
| element 7 of (ya 


(2) for some C, 0<CU<a 


(l—||AK-"!|-L)C 
S(y") e 
y _ 
will conve rgetoa un 
The are 7 will satisfy eq 
uniquely in (y). 
Thus when K is independent of y, only an initial 


| approximation to the solution at one point and 


approximate differential of S are required to inswre 
convergence. Consequently, this last theorem jis 
seen to be just the classical implicit function theorem 
in general analysis. One notes that in this case the 
continuity of S is not required. 

This section will be completed by establishing 
conditions under which the sequence of estimates 
will converge quadratically. 

(IV). Let y® be an element of Y and let y and L’ 


be constants such that 


| of some ne ighborhood (iy 


K(y™; y—y™) —Sy)+Sy™) 
<L’|\\y—y ||? (28 


for ally, y in (y),. 

We may now state the following theorem. 

TueoreM 3. Let the initial estimate y and th 
positive constant y be such that (1), (111), and UV 
are valid. Then the sequence of estimates determined 
by (18) converges quadratically to a unique element 
la, ay. The are ¥ is the 
unique solution of eq (16) in (y™) 

To see that (IV) implies (II) on A a a<y 
merely = a small enough so that L’, —y |1<L 
for all y“ in (y) As (1) and dip are undis- 
turbed on the smaller ne ighborhood, all conclusions 
but that of quadratic convergence follow at once 
from theorem 1. The quadratic nature of the con- 


vergence can be established by the following calcula- 





tions. From (17) and (18) it follows that 
y?—yY§=K"'(y ;—S(_y)). at) 
Therefore, 
K(y; yo—y) Sy"). 


Thus by (28) 
S(y™) || <L’||y—y? 
Making use of this last result in (29) one obtains 
Pr Sei 


This establishes theorem 3. 


5. Some Properties of 67'(x, y;») 


In this section some properties of the function 
57. which was defined in section 2, are developed 
; ; ; 
with an eye toward applying the considerations 0! 











ii 





> — 








section 4 to Newton’s method. First, however, the 
space Y to which the ares y belong must be more 
precisely defined. The space } is chosen to be the 
collection (y;(r)) of all n-tuples of real valued fune- 
tions that are of class C’ on the interval a<z<b. 
Addition of elements in ) and multiplication of 
elements of Y by real scalars are defined in the usual 
manner. Under these definitions of addition and 
multiplication ) forms a linear vector space over the 
real numbers 


For anv y in ¥ let 


Yr) ( > Y(t) . ). 
i=] 


Then a topology may be introduced into Y" by 
defining a norm 


y sup |y(2)|\V sup |y’(r) (30) 
a<r<b 


a<r<h 


The fact that (30) satisfies the properties of a norm 
is easily verified as is the fact that  y“’—y), has all 
of the properties of a metric. It is also not difficult 
to establish that any Cauchy sequence of elements 
of Y must converge in the sense of the norm to an 
element of Y. In faet, convergence in )" implies 
uniform convergence in the ordinary sense of the 
sequence of first derivatives. Thus in the topology 
specified by (30) Y is a Banach space. 

The subelass of ares in ) that have the property 

y(a y h 0 (31 
will be denoted by the symbol 0” It is clear that 
is a linear subspace of }. 

In order to eliminate repetition when discussing 
certain properties of the norms of elements of Y,, it is 
convenient to have the following simple lemma 
available for reference. 

Lemma 4. Let y,n, and Q(x,y;n) belong to Y. Let 
there corre spond to each component Q, of Q a non-nega- 
tire number C™ such that 


sup QGilryyin) SC ln 


z 


Also let there corre spond to each component () of (’, the 


first derivative with re spect to x of (), a non-negative 


number C'’ such that 
sup Qi (x,yin <¢c* e n. 
rz 
Then there exists C>0O such that 
Clry3n) SCin 


By (30) either 


a 4 
x,y; n)\|=sup (9 /@?)< 
z t=1 


nn 
24 Sup ap <> C® \n 


or 
Signy 
i Q4,y5 0) sup (>> @, +) < 
5 a t=1 
Sy sup [Qi < Ce" |y 
i=] z i=] 


Obviously © can be chosen to be the greater of the 


n 
numbers 5) C'? and 53 CO’, 


t=1 t 

In addition to defining the space of functions under 
consideration, it is desirable at this point to state 
some assumptions that will hold throughout the 
remainder of the paper. Unless otherwise noted, it 
is assumed that the functions g(z,y,y’) and h(x,y,y’) 
are of class C? ina region 2? of Euclidean 2n + 1 —space. 
Any are of Y, all of whose elements (z,y(x7),y'(x)) are 
in 2 will be termed admissible. At least one admis- 
sible are is required to exist. This implies the exist- 
ence of a subregion Y5, all of whose elements are 
admissible. It is also assumed that the determinant 
Diy: (LYN) never vanishes in ?. In the case of the 
Euler equations, that is, when g and A are given by 
(12), the above assumptions imply that the deriva- 
tives f, and f,, of the integrand f(z,y,y’) must be of 
class C, in RP. 

In eq (6), which defines the function 67\(2,y;y), 
let the parameter y range over Y, and the variable 
n range over U7. The following result is then valid. 

Lemma 5. For each y belonging to Yo the transfor- 
mation 6T(x,y;n) = (6T;, . «> 67y) 0 on YU te ; 
and is linear in n and bounded at y. 

The fact that 67 is an element of l’ follows from 
the considerations of section 2. The linearity in 9 
is a consequence of the corresponding property of eq 
(5) and of the operation of integration. 

To simplify the proof of the boundedness we 
introduce a lemma that will be found useful in what 
follows. 


LemMa 6. Let q(x;y) be given by 


qQ\t,y A \D)Y iL) 4 bi(x)y; (a 


where a; j(Z), b,(x) (] de % mn) are continuous on 
as r<b and yu in 2 The n there exists a constant 


C>0 such that for each y in Y 


’ 


sup gir;y)' <C y 32 
I 
and 


sup [ “attsy)dt<(o ayy). 33) 
Z a 


On recalling that sup y;|<|.y|) andsup|y; < (y| for 


z z 
| all(¢=1,2,...,n) one quickly obtains the inequality 


| sup| q(25y) <(Qusup a;\)\\y 4 (SUsup by ily) =Cly 


J 





The second statement follows from the first) by 
means of the inequality 


. a . 
qeyy dt. 
1 a 


sup qty dt< sup ity dt 


| pon noting that 6g and 6/4, are of the same type 


as q(zjy) and applying lemma (6) one obtains the 


inequality 


6 ry Ssup 6 7, 


sup | [9 bh ds | 
; : [oy bh ds Jat 


< 2 64 | bh ds dt 


= 2 [ sup 64g, +sup bh ds ]« 


<Cc* n , By > 


A similar result follows for 67° Thus the conelu- 


sion follows from lemma 4 


6. An Existence Lemma 


In section 7 the existence of an are yn in UU that 
satisfies 


67 (7, y3n) | [és ‘bh ds |ar 
; “| [9 shud | at 
I— AJ « Ja 


L(z), (=1,2,...,n) (34) 


for an arbitrary function d(z) in LU’ will be discussed 
To facilitate this discussion it is convenient to have 
available the following lemma. 

Lemma 7. For each admissible are y and every 
ne=tuple of functions m,(x) continuous on ac 2<b 
the system 


bq (2, uj;n)— bh (t, yjn)dt mr) (7=1,2,...,n) 
1 


(35) 


has a unique solution mir) which is of class C’ and 
satisfies n(a@) 0. 

An alternate form for (35) to which known 
methods of proof used in constructing standard 
existence theorems can be directly applied will be 
introduced Let 


n> Dat N (s=- 1.2, ....%). (36) 


) ony 
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Because by hypothesis the determinant 


never zero, eq (36) can be solved for 9. The s 
tion is 
n,m (r,n,c p= 3.2. nl 


is a linear and homogeneous function of 


where x 
and ¢ By use of (37) 6h, may be expressed as 


bh h n Lh fr lr“ne? 


a linear and homogeneous function in 9, and 
Therefore, any pair n, 9 satisfving (35) are tray 
formed by (36) into a pair y, ¢ satisfving 


¢ | lh yy t h ¥ riit.n jdt mAr) (7 @ 3 
n | r(t.n, Cjdt—0 (7 =e | 38 
a 


The last ~ equations of (38) are obtained by int 
grating (37) and recalling that »,(@)—0 

Conversely, assume 7,¢ satisfy (38) the second 
condition of that set implies (37). Inverting this last 
system of equations, one obtains (36), and by substi- 
tution in the first condition of (38) eq (35) result 
Thus instead of seeking a solution of (35) which is of 
one may equivalently seek a continuous 
solution of (38). 

The system (38) has the form 


class (’ 


y(r) A, (Oy (tdt+-™m (x) (7 


where A,,, m, are all continuous on a<ax<b. From 
well known methods of proof which have been used 
bv Bliss [2, p. 274-278] and others for systems of 
differential equations which can be obtained by 
formally differentiating (39) it can be shown that 
(39) has a unique continuous solution satisfying 


y,(a4)=™,(a). Since 
Y¥n41(2) =n, (a)=—0 G=1 2. ...,8) 


the lemma is proved 


7. Existence of an Inverse of 57(x,y,») 


In this section the existence of a linear bounded 
inverse in 7 for the function 67(2,y¥:) is discussed 
An idea that plays a leading role in this discussion 
is that of conjugate points. 

Derinition 5. A point 7; is said to be conjugate 
to the point P?, on an are joining P?, and P, if ther: 
exists a solution of 


d 


-bg lr, u3n)—bhdz, y3n)=0 (¢=°1,23.-. 
i ia Shale ’ 


~ 








ng 
of 
of 


( ‘ol 











Jong this are, which vanishes at P, and P, but is not 
dentically zero between P; and P? 
\ consequence of definition 5 that will be apropos 


is stated in the following 


» our present purpose 


emia 
LEMMA S Let n i.) A _n) be solutions of 
L n 6a bh dt Ons, n a 0 
. ae 
co eee $1) 
Then 2b is conjugate to r—a if,and only if. the de- 


ferminant on hb)! is zero 

By lemma 7 the »{” (2) described in the hypothesis 
xist and are of class ¢ In addition to satisfy- 
ng (41), the »{” are n linearly independent solutions 
of (40). Thus the n-parameter family of solutions 
of (40) which vanish at z—a is given by 


“ah AT HAR er i ” J 


Setting »)(b)=0 in (42), one obtains the system of 


linear algebraic equations 
n,” ANC 0 i a ee (43) 


a, the system 
This is pos- 


In order for z=6 to be conjugate to z 
3) must have a nontrivial solution. 
sible if, and only if, n (b) 0). 

LemMMA 9. For each admissible y in Y 
points are not conjugate 67 (x,y;n) is al 
/ onto itself. 

By lemma 5, 67(z,y;) is an element of LU’ for each 7 
in 0’ and any admissible function y. It remains to 
be demonstrated that the system (34) has a unique 
solution 7 in Ul for each L(x) in 7 provided that the 
end points of y are not conjugate. By lemma 7, there 
exists a unique (2) of class C’ which satisfies 


whose end 
] mapping 


bq bh, dt 


a 


I’ (xr), na ) (a 


. 


Let »{? be the functions mentioned in lemma 8& 
Then by lemma 8 there always exists a unique set 


C,} G@=1,2,...,n) such that 
n(b) +92? (b)C,=0 1,2, n 
Consequently, by (41) 
ni(J ny (d n.?’ (r)¢ =1,2, ’ 14 
is a solution of 
1 
bqgi— bh dt=1 + ¢ t= 1.2,...,%) (45) 
Ja 


aand r=b. That n(x) is of 
That n(x) is unique follows 


which vanishes at 2 
class C’ is apparent. 


from the condition of nonconjugate end points. For 
if 7»? and »® both satisfy (45) and the boundary 
conditions, their difference must satisfy (40) and 
vanish at both end points of a<z<h, that is, their 
difference must be identically zero 


By integration of (45) it follows that (7) must 
satisfy 
| 6s. . | 6h ds ja C(r—a)=/(2r) 
7a /@ 
i - re |S tt) 
Letting z—b4 in (46) one sees that 
%, 7 
| [0 bh ds Jat C(b—a)=—1(b)=—0 
Ja Ja 
(a=1,2,...,m 
Therefore, 
: l Js , Ja , : 
( i 0qdi- bh ds dt (7 ee , ft) 
4 dite . Ja 


Consequently, n(2) satisfies (34). 

We conclude with a lemma that states the result 
we have been pointing at in this and the preceding 
section 

Lemma 10. For each admissible y in Y whose end 
points are not conjugate the inverse of 6T (x,y; ) exists 
and is linear and bounded at y. 

By lemma 5, and lemma 9, 67'(2,y;n) is a bounded 
linear transformation mapping a Banach space in a 
1—1 manner onto itself. As is shown in Hille [3, 
p. 28-29], a transformation having the above prop- 
erties must possess a bounded inverse. 

The inverse of a linear transformation is linear 


8. Differentiability Properties of 67'(x,1;7) 


In this section two lemmas that will be useful in 
establishing the fact that 67 is one of the class of 
operators that was denoted by K(y;n) in section 4 
are proved. The first lemma will help to show that 
67 satisfies condition (1) of section 4, whereas the 
relates to conditions (I1) and (LV) of the 
same section. First, however, we present a lemma 
that will considerably simplify the statement of the 
proofs of our two principal ones 

LemMa 11. Let plryey’n.&) be of the form 


second 


Playsy’ sn,€ ay (ryy nts by (yy ne, 


ey (ry née .g=1,. + of), (47 


where a,,. bs. ¢:, are all continuous functions on the 
sae region R on which gq and h are defined, y is an 
admissible element of Y, and n, &— are elements of U. 
Let i be any admissible element of Y and 6->0 be 
such that the closure of (y)s contains only admissible 
arcs. Then give ne 0 a follows that there exists d~>0 
such that 


(48 


ftr 


Play? yor’ n,&) — play? yo’ m,&) Ce in 
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m (y ”),, which satisfy y" yr ii" d 


for all y, y° 
Furthermore, there exists a constant 


and all n,& on / 
('>0 such that 


sup plryyy’ ne) SC on &, (49) 
and 
sup Axryy’ neds <(b—a' in + é (50) 
Ja 
for youn the closure of (y , nbn / 


Statements (49) and (50) follow from a simple 
extension of lemma 6 
In order to establish the remaining part of the 


lemma, let 
xe.y”’ 9 .Be 
at het head © 3» 


Then ¢ is a continuous function on the closed and 
bounded set S:a<a<bh, y y™ in the closure of 
y Ys, on é 1. That it is sufficient to merely 
consider 9, — on the unit sphere follows from the 


s 


bilinearity of p. It is clear that g=0 on the subset 
of S for which y? =y Thus ¢ <¢ on a neighbor- 
hood of this set, that is, on a<zr<h, gy”. im 
ys, ly—y ||<d, |\n £ | for d sufficiently 
small. This proves lemma 11 


each admissible y in Y there 
0 such that 6 T(x, 450) is of class (" 
in y uniformly on (y)s 

Let 6°g. (2, y3n, &) and &A, r.y;n, &) be 


LemMaA 12. For 
erists a constant 6 


defined by 


Pd =i (LU Unk Diy ole UW ngk 

q n,§ d n,t (51) 
and 
Pho sh r, uu’ )nyt h ni 


Negty. mh 
where y is some admissible element of Y and n,& are 
ares in Ul. Then the function d(67')) mentioned in 
definition 4 can be defined as 


. bd 
TT (r.uin, 8 | [ ea - ihads Jat 
7 — | | ea shade |at i = “2 on 
b—aJa , Ja 


By reasoning similar to that used in considering 67 
it follows that #7 is on YUU’ to l That &T is 
bilinear in » and € is obvious. To complete the 
proof of the lemma conditions (1), (2), and (3) of 
definition 4 must be established. 

(1) For every admissible y® in Y and every finite 


d 57’,) 


5>0 such that the closure of (y)s contains only 


admissible ares #& T(z, y;n,£) 8 continuous in y : 


0) 


formly on (y)s. 
Let F(z,y,y' 3n,€) be any function defined on } 
LU’. Then we define 


AF (1,2) =F(az,y™ yo” jn) — Play ye’ snk 
Keeping this convention in mind, one may writ: 


sup 6°7) (2, y'"; n,£) 


sup | [ a0 q (1,2) | Ad*h (1,2)ds ] dt 


—— [ 36 q,(1, 2) | Ash 1, 2)ds | ae 
hb ade , Ja 


< 2| Aé"4, (1, 2) | a5%,(1, 2)ds dt 


. 2| Ad*4 (1. 2)idt+-2(6b a) | Ab7A, (1. 2) dt 
Ja Ja 


Examination of (51) and (52) reveals that both &y 
and 5h, are of the type p(z,y,y'.n,£) defined in 
lemma 11. Thus application of lemma 11 to (53 
vields the conclusion 


sup aT, watt PT (2, y -n,t en-:eé 


ls Which satisfy a condition of 


for all y, y® in (y° 
The same result can be 


the form y?—y <d 


established when #7, is replaced by &7;.  There- 
fore, (1) follows at once from lemma 4 

(2) For every admissible y an Y there exists 
5>O such that &T(xz,y3n,£) is uniformly bounded or 
(y? Ds. 


From an application of (49) and (50) of lemma 11 
the inequality 


” 7 


sup 87, (z,y™; 9,8) : 2| bq; bh.ds dt 
! 1 Ja 


q 2| [ sup dq, + sup | bh.ds qa 


<C|\ yi! -lell, C’>O (54 


From analogous considerations the 


-O such that 


is obtained 
existence of ¢ 


sup #7; (4,439, &)) < CO" \n 


fre 


can be demonstrated. Because of lemma 4, in- 
equalities (54) and (55) imply that &7 is bounded 
inyaty®. By (1) &T7 is continuous in y at y 
Thus by lemma | there exists 6>0 such that #7 is 
uniformly bounded in y on (y)s. 

(3) For each admissible y® in Y and each finite 
6 >O which is such that the closure of (y)s contains 
only admissible ares, the function R (yy sn 
fined by replacing F by 6T and dF by &#T in | 











wh 


D 





n 








such that given «>0 there exists d>O for which 


Rly :, Yy' sn) € 7 


wv all -, in al satisfying y y <i 
ind all nin U 
Consider 


[Agi (1,2)n 


t 
(Ah yy (1,2)n) + Shy, (1,2)n/)ds]dt 


J a 


h a. 
Ahy, (1,2)n) 4+ Ahy (1,2)n)ds\dt. (56) 


\pplying the law of the mean to a typical A-symbo!l 
n (56), one obtains 


where 


Y y~ M\Y v] ; O< pu a 


Thus by applying the law of the mean to the remain- 
ng A-svmbols and taking note ef (51) and (52), we 


may write 
Pry sn) 6T (ay? 3) FT (ay ny? —y™) 
“D r,yin,g)dt—* ” Dat, 
. b a, 

where 

D=Aq 3,1) ly y Ag Ln yy y 
Ag aly —y Aq 3a (yy y,”) 

Meiy y (3,0 a(ye? —y Meiy (3,0) 05 

; y Ys, 
Mey (3, Daye? —y Ah 3,1)9/ 


oN ye Ids. (57) 


Of course, the 3 appearing in the arguments of the 
various A-symbols does not always refer to the same 
value of yu 


The function gy.95(M yo’), that is typical of 
the funetions to which the A operation is applied in 
57), is of the type considered in lemma 11. Thus 
| readily follows from that lemma that, 


sup 6 ¥ rey? sn)—s Ti(x.y =) —i 7 (r,y so 


y) <2{"|D J, Das\ae. eliaii- iy —y~ 


for all y®?, y™ in (y)s such that |\y y*\|<d. A 
similar result holds if 67, and 6°7, are replaced by 
their derivatives. Therefore, by lemma 4 it is true 


that 

Ry youn -y yc 
6T (ry 39) —8T (ry 59) — PT (ayy 3 ny — yy" 
en: ¥y y . (58) 


forall y?, y@ in (y), that are close enough together. 
Statement (3) follows at once from (58) 

By choosing the 6’s of (1) and (3) to be equal to 
each other and less than the 6 mentioned in (2). 
one obtains a 6 that is suitable to be used in lemma 
12 

LemMMA 13 kor any admissable yO in VY there 
erist constants a >0O, L’>0 such that 


6T (x,y sy y" T(rsy™) + Tay? 
for all y, y° 
Consider 


T Aa sy Ta atl 


[ a9 (2,1) | Ah (2,)ds Jat 
U ada Ja 


By the law of the mean 


Agi (2,1) = Gy, (z,y"’,y ) (Yi? —Y 
Ii," (2,y"" y 7] 
where 


y y p(y™—y™), Ocul. 


As a similar result holds for Ah,, we may write 


bT (x,y sy y" T(xsy)+T (rsy? 
a . J a "0 
Alz,uin, dz Zdr, 
Ja : 4 a, , 
where 
Z=AgiA1,3)yl? —y)4+- Aga: 3)(y$ y 
[Ah 1,3)(y)' ys Ah ,,7(1.3 y " lds, 


Of course, the 3 appearing in the arguments of Ag 
and Ah, does not always refer to the same value of ug. 
Application of the law of the mean to a typical term 


of Z yields 





where 


y y u*(y y“), O< pu l. 
It will be noted that Agy/(1.5) is of the same form 
as the funetions considered in lemma 11. Let @ be 
chosen small enough so that for any y“’, vy in (ya, 
y lies in some fixed neighborhood (y/), whose 
closure contains only admissible functions. Then 
because of our previous hypothesis on the class of 
h, and g;, we may apply lemma 11 from which the 


existence of C >0Osuch that 
sup bT (ry sy y~ T(ayy T (arcy” 

( y y 
follows. As similar considerations are valid when 


67,7 are replaced by 67” .7 . respectively, the con- 
clusion of lemma 13 follows from lemma 4 


9. Convergence of the Sequence of 
Estimates Derived by Newton's Method 


The main theorem on the convergence of the se- 
quence of estimates derived from Newton’s method 
is the following 

THrorem 4 Let y be an admissible element of Y 
having nonconjugate ‘ nd pownts The " uf there i ists a 
suficv ntly emall constant [)>0 such that T(x.y ) DD) 
the sequence of estimates derived from Newton’s method 
will converge quadratically to a unique are 7 of class C 
which lies within a prede termined ne ighborhood of y 
and satishe S T(r; “ y a 7] a Y b y (h 

By lemma 12, there exists a constant a >0 such 
that 67 (2,y;n) is of class C’ in y uniformly on (y) 4, 
By lemma 10, 677 '(z,y";) exists, and is linear and 
bounded at By lemma 5, 67 is on (vy), U to U 
and is bounded in y and linear in ». Therefore, by 
lemma 3, there exists a such that 677! exists on 
(vy), and is of class C’ in y uniformly on (y).,, 
Thus by lemma 2, there exist nonnegative constants 


a, and VW such that for y™, y® in (yy). and 
" y™) <,. it follows that 

67 (x,y 39) —8T (ny 59) SM oy yous» 
Consequently, if we choose a, = ay/2, a part of condi- 


tion (1) of section 4 Is satisfied on (y "4 
By an argument similar to that applied in section 
t, it can be shown that there exist nonnegative 


constants a;, B such that 
6T'(r7,yin) <Bon 59 


for all y in (y®’)., and all yin U. If ag is taken to be 
less than or equal inf (ay,a;), and such that (y®),, 
contains only admissible arcs, the part of condition 
(1) previously proved and inequality (59) are still valid 
on (y) ag, that is, condition (I) holds on (y’),,. 
The existence of nonnegative constants 17,“L, N, 
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which along with the nonnegative constants /J 
already determined in this section satisfy rel; 


(21) is clear. By lemma 13 there also exist a, an 
such that 


bT (ry sy y Tryy V(ryy < 

ky i y 
for all vy, v® in (y 4, That is, condition (1 
of section 4 is valid on (y),.. Thus we may cho 


a, <a; such that conditions (1) and (IV) both hy 


on (y¥)a.. By continuity there exists ag such t) 
Tixz,y)|<D> on y), Let C be such ¢] 


‘ ‘ 9 


O< C<inf (ag.ay Then if the constant D is | 
than or equal to the smaller of N and (1—/H/)¢ 
condition (II1) is” satisfied on (y), Con 


7) 


quently, if the positive constant y in theorem 3 


defined to be the inf (ag,a 9), we may conclude by tha 


theorem that Newton’s method converges quadratic- 


ally to a unique element 7 of (y),, a<y, which 
is such that (3) is satisfied. That 
yla y'(a), Wb y bh 
follows at once from the fact that all the variations 
7 vanish at r=—a and r—b 
The are 7 is of class C’ by virtue of being in Y 
To see that 7 is actually of class C° consider th 
equations in (2,7 
G r.0) adr, yr), v) h(t, yd, 7'(O)dt 
| [ a oA 
oldr. ular) yota 
4 a Ja . } / 
Ads, 8), 7'(s ds \|dt ) ; ey oat H0 
Ja 
By what has just been proved the point (x, 7’(r)) is 
a particular solution of (60) for each « on a<xe</ 
Since for each such x, (7,9(4),7’(x)) is an interior 
point of the region ? in which g and h are of class ( 
G(x) is of the same class in a neighborhood of each 
point (7, 7’ (2) By hypothesis the Jacobian 
G,, g 1) 


| 


in the region ?. Thus the implicit function theore: 
(see Bliss [2, p. 269]) is applicable. This theore: 
tells us that 7’(r) is of the same class as G(s 
that is, 7(x) is of class C°. 

The following theorem gives further conditix 
under which the conclusions of theorem 4 will 
valid 

Tueorem 5. Let y* be an admissible are with 1 
conjugate ¢ nd points such that T(x;y*)=6. Then th 
erists a constant a >0 such that for any initial esti 
y” taken from y*), the conclusions of theorem 4 


valid. 


’ 
; 
‘ 


‘ 




















to that employed in the 
proof of theorem 4 there exists 8 >0 such that condi- 
tions (1) and (IT) and (IV) are satisfied on (y*), 
By continuity 8 may be made to have the additional 
property that  7(y)| <N for all y in (y*) Let a 
be such that 0<a< 8/2. Then for any y belonging 
to (y*)q it follows that (y)s,. is contained in (y*)s 
Thus (1), (11), (IV), and the inequality | 7T(y)) <N 
hold on (y)g;2. Let C be such that 0<C< 86/2 
Then by continuity a may be taken so small that for 


By reasoning similar 


lly? im (y*)a Ty) ||) <0-IDC B. This says con 
dition (111) also holds on the (8/2)-neighborhood of 
any y taken from (y* We now invoke theorem 3 


and the remainder of the proof follows as in theorem 
} 


If one wishes, he may hold fixed in 67 (2,4; 
throughout the entire process. Theorem 2 is ap- 


plicable in this simplified case, and with its help all 
of the previous results can be shown to be true for 
the simplified Newton’s method. 


10. A Lemma Concerning Conjugate Points 


It is easily seen that the equation 
57 (x,y 7 [és t.yin bh (s,y5 nds | at - 
. i ea 
r—a {(* *t 
g bh ds Jat 0 
h a Ja Ja 


) 


(i ea (61) 


and eq (40) of definition 5 are equivalent. Thus in 
the considerations of this section it will be assumed 
that (40) has been re placed by (61) in the definition 
of conjugate points. 

If y is an are such that 67~'(2,y;n), the inverse 
with respect to » of 67(x,y;n), exists for all y in U, 
then y has nonconjugate end points, because the 
unique element of U' which 6 7T(z,y;n) maps into @ is 6 
itself. In lemma 9 it was shown that if y has non- 
conjugate end points, 67~'(z,y;n) exists. Thus a 
necessary and sufficient condition for 67™'(z,y;n) to 
exist for all » in €’ is that y have nonconjugate end 
points. However, by lemma 3, if 67~' exists at y 
it must exist in some neighborhood of y®. There- 
fore, the following lemma may be looked upon as a 
corollary of lemma 3. 

Lemma 14. Let y® be an admissible element of Y 
having nonconjugate end points. Then there exists 
56 >0 such that no elements of (y)s have conjugate end 
pornts 

In this section we wish to give an alternate proof 
of lemma 14, which depends directly on the proper- 
ties of functions with Lebesgue square integrable 
derivatives rather than on the Banach space theory 
of section 3. 

To construct a proof of lemma 14 by contradiction, 
assume that there exists no neighborhood of y™ in 
which all elements have nonconjugate end points. 
Then in every neighborhood of y there exists at least 
one function y suc ch that (61) has a nontrivial solution 
along y which vanishes atz=aandr=6. Consequent- 


ly, from consideration of the //q-neighborhoods of y“ 
for q : 


1,2, ... one obtains the sequences (yy), (9), 
which are such that 


uniformly on a<s<b, 


ye (r)ynn"(z) 
q 2 oe. ) 
(2) 
y@"(r)—y'"(2) uniformly on a<ar<b, 
gq=1,2 
57 i 0 q=1,2, 
(63) 
n a n h 4 q 1.2, 


— all ares » under consideration belong to the 
they must certainly belong to the class of the 


spac c 
n,)-space defined by 


rete of ares in (x, m, . 
a set of n real valued functions 
nr a<zr<b; i1=—1,2, — 

that are absolutely continuous and have integrable 

square derivatives n’(z) on a<r<b. Henceforth, 

we will think of the y’s as belonging to this class. 
Without loss of generality one can choose 7” such 

that 


(no sum on q). 


(64) 


Equations (64) imply the existence of a subsequence 
n%’ >, which converges weakly in L? to some element 
7’. We can assume that this subsequence is the 
same as our original sequence. The weak conver- 
gence in L? of the sequence | * implies 


n> uniformlyona<zr<b (i=1,2, n) 
(65) 
7? 9," weakly in L? (i=1,2,...,n). 
It is also true that 
nen? dx <lim inf | 2’? ‘de (nosumon q). 
Je reas 


(66) 


A proof of this last fact has been given by MeShane 
[4, p. 355}. 

Inequality (66) will be used in the proof of the 
following lemma 

Lemma 15. Let (n®’. be a sequence in L? such 
that (64) and (65) are satisfied and let Aj? (x), Bi (zx) 
(¢,j=1,2, . . . , n) be continuous functions for almost 
all q such that 
Ay A$ 


uniformly on a<r<b 


(G7) 





BS —Be uniformly on a<2r <b 

D A n) 
Then if 
W(x) (AS? 9“ + B® »\ lt q a See 


no sum on q), 
if follows that 


WYows uniformly on a<.s: b 
te 4 n) 
(OS) 
W?'’-W? weakly in L 
] 2, vi) 


Consider 


1Wer) W (zx) 


A® 9 )dt\ 4 A 9 —A®_®)dt + 


*, . 


{ 
| (Bane Bo n/ lt + BY n° Bone dt- 


The first of the four terms on the right of the sign of 
inequality satisfies 


| (Alene Ai? ni? dt < | A® —A®|.\9@ {dt 


(b—a) max |Ajy —A,° - |}! 
By (67) max |A\7 — A,” |-0 for each ((,j7= 1, 2,...,n). 
By (65) |n{*| is uniformly bounded for all 7, q 
under consideration. Thus the first term goes uni- 


formly to zero on a x <b. The second term on 

the right satisfies 

tes - 

|| (Aa —AM a ydt < | AM |-\ne — a \de 
max {° n," ni dt 


This goes uniformly to zero by (65). For the third 
term we have 


| (BS n," BO» dt < B Be “IMG” dt 


<max BS 


By the Schwartz inequality and (64) 
( | ni dt) <(b a)| ne “dt<(b—a) 


a a 


(q=1,2,.. .). 


Thus n\* 
. uy 
(q=1,2,...). Consequently, one may conc! 
from (67) that the third term goes uniformly to 
on a<az<b. 
Only the fourth term remains to be discussed 

j be temporarily held fixed. As is well known t} 
exists for each arbitrary ¢€>0 a polynomial P(z) s 
that 


dtis uniformly bounded for (j= 1, 2, 


a @ 
(Be r)—P(xr))da Ss 


16 


By (64), (66), and the Minkowski inequality it js 


true that 


da )< 2 70 


Thus by use of Schwarz’ inequality, (69), and(70 
one obtains 


nn” )at P(t)(n nn” )dt 


/@ 


i) Bin 


1 


<( [iar—rray( for aera) o 


3 
3 


<( | (ae Pydt)( | 


for allz ona<ar<h Integration by parts \ ields 


[m0 


Therefore, it follows from the uniform’ convergenct 
| of (yn —n, ) to zero that 


| P(H)(n’ nn” )dt < uniformly on a<r: 


for gq sufficiently large. 
imply that 


| B& (ng? —n°')dt-0 uniformly on a<zr<>. 





This completes the proof of the first statement of (| 








Formulas (71) and (72) 








‘)} 


fu 


an 
aun 














Let / be any element of Z* and consider 


( Wwe'lds | W'lda 
| (AG a 1% 9 da | Be 9?’ —B® _’ Vida 
max || Ayn — Asn! \da 
I, Ben’ — BY \dz' 
| IB (y n®')dx 


By a repetition of the arguments already given above 


the terms multiplied by max |) go to zero. The 
factor / may be lumped with 2 without changing the 
reasoning used in analyzing term four. Thus the 
second statement of (68) is seen to be true 


The following corollary is a direct result of lemma 


COROLLARY 1 Let the hypothe sis of lemma 15 
hold and let 
U(r) =We(r)—7—"* We(b) (73) 
h a 
, “Z r—a [{(°.,.. ea 
V (2 W 2 (t)dt WW t)dt. i4 
Ja h a 1 
The 
/ l uniformly ona<a<b (¢=1,2....,2 
} V uniformly ona <2 <b 1,2, n) 
l U weakly in L (soz 1:2....,8) 
\ Ve" uniformlvona<r<b (i=1,2,...,n). 
Let 
i(ryin tioht~~— dbgdt (=1,2, hn 
Ja . h a. a . 
ruin ( bhids ) dt 
z—a [{*" *t ° 
6 i,ds ( (7 2 i) 
of (flehan)dt Gta. 


Then u and v are elements of U’ satisfying 


67 (x,y jn) = ue 


- 


Since bg (t,y@;n@)dt and | bh (t,y@;n@)dt are 


Ja Ja 
functions of the type W(x) defined in lemma 15 for 
each admissible y@ in Y and 9’ in L?,udaz,y@;n™) | 
and o(z,y;9) are respectively of the types (73) | 
and (74). Consequently, if we take {y@} and {9@} | 
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to be the sequences which result from contradicting 
lemma 14, the hypotheses of lemma 15 will be fulfilled 


because of (62) and (65), and corollary 1 will then 
imply 
udr.y sn ulry sn 
uniformly on a<a<b 1,2 ! 76) 
OAT, ys” Ee AL LY a - 
uniformly on a<2<b (¢=1,2 " 77) 
“u(r.y n Ul(r.Yy n 
weakly in L pax 3.2. a7 78) 
; LY a] >) ryy a] 
uniformly on a<a2<b pax} 2 n 79) 
In view of (63), formulas (76) and (77) vield the 
result 
6T La? se i] (SO) 
The limit function 7 (xr) is clearly a continuous 


function on a<2 <b satisfying 


U] a n (h ad 


in 7 


To show that » is of class ¢’ that is n is 
observe that by differentiation of (SO 


Jy Iw 
%, : ] */ A bd 4 ‘ 
bh dt + Od; bh ds dt S| > 
Ja b a Ja . Ja 
Since Igy} #0, (S81) can be solved for As all 


other quantities involved in (81) are continuous, » 
must have the same property 

By the hypothesis of nonconjugate end points the 
only element of (7 satisfying (S80) is the zero element 
Therefore @and (78) and (79) may be amended 
to read 


u(r, y + ae a | ;7 ) 
0 weakly in L? ph ZF... ccetth (SR 
v(r,y@ 59) ery ;7 ) 
0 uniformly ona<r<b (41,2, n). (83) 
Consider 
lim | BT. (2, yn )6T (2, y@ 5 n@)da 
q-@,"a 
lim | “~ u,* 2u;" v -o,* ¢ da (S4) 
q-oJ a 


Although we have not previously defined all the 
notations in (S84) their meanings are obvious. By 
(82) and (83) the terms involving v, may be disre- 





garded By lemma 15 


= 
lim uc? uc?’ dr 
7 a i 


lim] fag, bq,¢ das ; m [, 0 da )( i} by.° dr )] 


lim by.* 6g.° de. 


Multiplying out, one obtains after droppmg the 
superseript (gq) from the derivatives of g 


lim | bq bgda 


lim (diy9," Ding ™ “FiyMs? Dev, 
Iu9s dun’ de 
lim dy denn? ne deslim] Gene ne da 


k 
. 


since all terms involving »® go to zero 


~— 
Let (A 


(A, Then (since by hypothesis det (9; 


matrix 
) #0) 


denote the transpose of the 


a ' 
det (G2) det (9.y ydet( 7d -) [det (9 " )}>0. 


Hence Gy 
matrix for each g=1,2, 
that 

lim | 67! (2,y@;3n" 


» 


is a real, symmetric, positive definite 
Therefore, it follows 


ST (2 y@sn@ da 
a 
lim} Gini’ ne’ dr>C>0. 


q-@J« 


However, (63) implies that this last limit is identi- 
cally zero. Consequently, our assumption leads to 
a contradiction. Therefore there must exist at least 
one neighborhood of y® in which all elements have 
nonconjugate end points 


ll. A Solution of the Fixed End Point Problem 


To conclude the discussion of Newton’s method, 
we wish to point out that with certain additional 
hypothesis we can make use of our previous results 
to obtain a solution of the fixed end point problem 
of the Calculus of Variations. 

Let f(z,y,y’) be the integrand of the integral (1). A 
strengthened Legendre condition is said to hold in the 
region FP if 

Fuge (2.Yy em >() (85) 
for all (z,y,y’) in R and all (+) #0. This implies 
that the determinant|/,,,) #0 in # and is thus 


condition that 
never vanishes in &. 


stronger than our 
determinant 4,, 


previous 


The following terminology is well known and 
be helpful. An are y is said to satisfy a strength, 
Jacobi condition if it has nonconjugate end po 
and has no subare defined on a<2<¢<b which 
conjugate end points 

Let (2) be replaced by the Euler equation 
Let the region 2? be convex in y’ and let the streng 
ened Legendre condition hold everywhere in 
The following theorem is then true 

THEOREM 6 Let y be an admissible element 
Y having nonconjugate end points Then if ti 
erists a sufficiently small constant D>O such t) 

Ti(xy o)) ‘ Dd. the seque nee of estimates de rived 
Newton’s method will cOnvErge quadratically to a 
unique are Y of class ¢ which lies within a predeter- 
mined neighborhood of y” joins the end joints f 
y", and satisfies the Euler equations. Moreover, if 7 
is such that no interior point of a<ar<b is conjugate f 
g=@, y will afford a strong re lative minimum to I(y 

By theorem 4, Newton’s method converges quad- 
ratically to an are 7 of class ¢*, which joins the end 
points of y and satisfies the Euler equations. Thy 
arc 7 also has the property of nonconjugate end 
points. This follows from lemma 14, which implies 
that the constant a, mentioned in the proof of 
theorem 4 can be chosen to have the additional 
property that no element of (y)4, has conjugate 
end points. Since by hypothesis no subare of 7 with 
initial point at s=a has conjugate end points, a 
strengthened Jacobi condition is satisfied along 7 
As it has also been assumed that a strengthened 
Legendre condition holds throughout ? and that P 
is convex in y’, we may conclude (see Bliss [2, p. 42, 
corollary 16.1]) that 7 affords a strong relative 
minimum to /(y). 

The following theorem gives further criteria that 
insure the convergence of the sequence of estimates 
derived from Newton’s method to a solution of the 
fixed end point problem. 

TueoreM 7. Let a minimum y* exist and have 
nonconjugate end points, Then there erists a constant 
a>0O such that for any initial estimate y taken from 
(y*)q the conclusions of theorem 6 hold. 

Upon recalling that a necessary condition for y* 
to be a minimum is 7(z;y*)=8, one proceeds with 
the proof just as in theorem 5. 


12. A Second Method for the Effective Solu- 
tion of the Fixed End Point Problem 


In section 2 it was remarked that in the application 
of Newton’s method, one must solve the Euler 
equations of each term of a sequence of integrals 
whose integrands are quadratic in y and y’. Thus 
it is not inappropriate to introduce at this point a 
second method for solving fixed end point problems 
which will prove to be particularly applicable to the 
case of a quadratic integrand. In an actual compu- 
tation a combination of Newton’s method, and the 
method to be presented below, might prove to 
be advantageous. 
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Henceforth the symbol Y will denote the class of 
-tuples of real valued functions on a<z<b which 
have Lebesgue square integrable derivatives. Use 
vas made of this class of functions in section 10 
The symbol U will again denote the subclass of 
lements of Y, which vanish at the end points of 
x<b. Throughout the remainder of the paper 
assume that all functions of Y are admissible and 
that the underlying region F? is convex in y and y’. 
Our previous assumptions as to the nature of the 
ntegrand f(z,y.y') are retained Moreover, we 
assume that f and its derivatives are integrable 
along ares in Y 
Assume that after starting from some initial 
estimate y® in Y one has obtained an “improved” 
estimate of a solution to the Euler equations y 
iso 1n Yy He then chooses the 7 Ist estimate to be 


no sum on 7), (S6) 


vhere 7 is an element of (° determined by the 
juation 
7 =T(x;y). (87 


Here 7 is formed for the Euler equations as illus- 
trated in section 2, and a, is a real number. 

\ procedure similar to that just deseribed has 
been announced by L. V. Kantorovitch [5] for the 
case of an integral whose integrand is quadratic in y 
and y’. However, his paper presents no proofs. 
Curry [6] has described an analogous process that 
an be used to minimize a function of n real variables. 
Other references may be found in the last paper cited. 

In order to consider the convergence of the 
sequence defined by (86) it is helpful to examine the 
expression I(y)—I(y— an) We may state the 
following lemma 

LemMaA 16. Let y be an element of Y and let n be 
oven by 87). Then. 


I(y I(y an)—a n.n.dgx ar Wdar (88) 
: a 


uf he ré 


, 


why» (2, y—aun, y’—apun’ dy 


and By, Cy are obtained from the formula for Ay by | 
replacing fy, by fy.y: and fy-,, respectively. 
saline ng | atlas | sled ted | ‘ 
An application of Taylor’s theorem results in 


I(y)—I(y= an) 


a| LA (2.yy’ met fy (asysy’)nilda—o8 | Wade. 
a ® . /@ 


By integration by parts 


PL savedsiees | 


f Syndr 





and since 


-; 


| C.n'dxz=C,[nb)—nAa)|=0, 


/@ 


Accordingly . we choose 


l - : odd 
( [ (t,y.y’) / (ssyy’ods | at 
bh—a Ja : aie: — 


Thus one may write 


the choice of C; is arbitrary 


I Wf) I(y an a Tinida a, | Wada 


a “* 
a n.n da a | Wada 
Ja Ja 


The next lemma follows directly from formula (88). 

LEMMA 17 Assume that for every number v there 
erists a number K->O and a number a*® such that if 
I(y <v and O<a<a* then 


” Wade <K [" ainide. (89) 


Ja Ja 


Then in case Tiy)<v and O<a<a* one has the 
relation 
I(y)—Ity 


an) > all ak n.n.da (90) 
i 


Formula (90) helps to prove the following theorem 

Tueorem 8. Let there exist finite numbers vy, 
K>0, a*® such that if T(y) <v and 0<a<a*, formula 
SG) holds Let y ©) be such that ] gas \< vy and let 
a and @ he chosen so that O0<a<a<1I/K Then if 
I(y) is bounded helow on the class of ares mn Y joing 
the end points of y, the sequence of estimates {y} 


defined by (86) and (87) ws such that {I(y'?)} is 
convergent. Furthermore, 
T(xsy ‘6 uniformly on a<sr<b (91) 
T’(x,y)—0 in L’. 92) 


That @ may be chosen as specified in the hypoth- 
esis is a consequence of the fact that A is finite. As 
a result of such a choice of a, one obtains the 
inequality 


I(y) Ty” aon” Tiy(@ ) 

from (90) by setting y=y, a=ao. Since [(y'")<», 
the argument can be repeated provided 9 #8. 
Hence a sequence (y“’} such that 


I(y)>T(y**) (t=0, 1,2, 
is constructed. Due to the hypothesis that /(y) is 


bounded below the sequence of definite integrals 
converges. 
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By (90) and the way | a, is chosen 


> (1 aK | n ” da 


Ty Ty" ») 


” 


-a(l—A sup a n 


where here and in similar expressions that occur 
below there is no sum on the index that denotes 
which step of the iteration process we are considering 
As the left-hand side of (93) can be made as small as 
desired by choosing 7 sufficiently large, we may 
conclude that 


lim ( [ ” dr) ()- (94) 


Conclusion (92) follows at once from (94). Since 
Y1) Is a consequence of (92), the theorem is proved. 

The problem of determining general conditions 
under which an assumption of the type made in 
lemma 17 is valid seems to be a difficult one. How- 
ever, special cases may be considered in which K is 
independent of » and a* is arbitrary. A rather 
trivial example of this occurs when the integrand is 
given by 


f(r.y.y’) y(r)yl+y”, g(r) >0 


ona <a<h, we have 


If 2A > gla 


, me | , 
, , Grn - 
Woda (1 wry . 77 
i . AL ‘ 


dyudi 


Ss K nda: 


/@ 


A much more important case is that in which the 
integrand is quadratic, that is, 
ry + Cy(2 yy 


f(z.y.y’)=—Ay(r)yy,+2B 


Dixy +k (r)yi+F f(r). (95) 


For f(z,y,y') given by (95) we obtain upon integrating 
out yp 


Wde LA, (e)niny + Bye) ney’ + Cy (a2) a nda 
7 a 
<K] ninids 
Ja 


for a suitably chosen constant A provided 7,(a) 
n,(6)=0. This last result follows from the Schwarz 
inequality and the fact that the integrand is a 
quadratic form in » and »’ having continuous 
coefficients. 

By setting y=y", »=y7-—y°, a=—1 in the 
computations employed to establish lemma 16, one 
obtains the formula 


Ty)=I(y)+ | Tihasy \nidrt+ | | Wdzr, (96) 
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which is true for all 7 in Y which join the end p 


of y®. This equation is helpful in the proofs | 
piven below. 


LemMa 18. Let there exist a constant k>O sue} 
Wdr>k n.n.dx 
Ja Ja 
for all n in U7 and all Y in Y. Then I (y) is bour 


below on the class of elements of Y joining the 
points of some arbitrary fixed y in Y. Furthern 
given a constant v there exists a constant M such 
for each y m Y jon ing the end points of ye 
I(y)<v 


for wi 


ninida<M. S 


In (98) we have g=y—y" 
By Schwarz’ inequality 


| Ti(xsy ; \n. 
Ja 


( "Tasy IT (2; y" da an n nda ) 


Thus by (96) and (97) 


Iiy)>Ty os h ( [ ainidr) +k [nintds 


h? h Be de iF 
E+ 57 (4 [ainda] 


> Ty f 


T(y) 


h 
th 


provided y joms the end points of y. This proves 
that J(y) is bounded below on the class of functions 
joining the end points of y®. The remainder of the 
lemma is also a ready consequence of this last 
inequality. 

Condition (97) does not seem to be very widely 
applicable. However, in case the integrand is 
quadratic the strengthened Jacobi condition in con- 
junction with the strengthened Legendre condition 
implies the existence of a constant k>O such that 
(97) holds. 


Lemma 19. Let (y) be a sequence in Y such that 


yy, uniformly on a<r<b, ((=1,2,...." 
99) 
y?'-»7/ almost uniformly (s=1,2. . n 
(100) 
Ti(2;y)-0 almost uniformly (t=1,2, ...,) 
101) 
Then 
T(x;97)=0 almost everywhere (i=1,2,...,n 

102 








Col 


the 





i- 
1 
it 


~ 








Let J denote the interval a<zr<b. Since 7’ isin 
it is summable and hence finite almost everywhere. 
Hence given e>0 there exists by (100) a set S, of 
measure less than «/2 such that y{?’—-7 uniformly 
un J—S, for ( . # n) Of course, S; must con- 
iin all the infinites of 7’. The summability of 7’ 
lso implies the existence of a set S, of measure less 
han e4 such that 7 is bounded on J—S, for 

1.2... .n). Therefore, if Ss is an open set of 
ieasure less than e€ containing S;+ SS), the following 
atements are valid: 


a) Z—S, is closed and bounded, 


b) 7 is bounded on I—S; (¢=1,2,. © .n), 
c) y -y, uniformly on I—S; (1,2, sn), 
d) y -y, uniformly on J—S; (¢=1,2, wn). 


By (a) through (d) it can be concluded that 


2.” 3 ~f (1,4, 9) uniformly on I S, 


and similarly for f,. Hence it readily is verified that 
Diecy@)-Ti(z,9) uniformly on I—S, for (7 
3. n), that is, almost uniformly and thus almost 
everywhere. Conclusion (102) is then obtainable 
from (101 

If one wishes, he may express (102) in the form 
\x;7)=—8. This follows directly by integration. 

LEMMA 20. 
102 Then if there exists k->0O such that (97) holds, 
7 affords a unique minimum to I(y) on the class of 
functions in y which yom the end points of 7%. 

Let y be any element of Y joining the end points 
of 7. By (96 


-) 


Ti(rsylyi—yDdr- ’ Wdr. 


Ja 


By hypothesis 77(2;7)—0 almost everywhere and 
y,—¥, isin L*. Thus the first integral on the right 
is zero by Hélder’s inequality. Therefore, if y #7, 


I(y)—1H) Wdr>k] (yi —9Y—Wdzr>0. 
Jf “a 
( ‘onsequently . 
I(y)>I(9). (103) 


\ssume that y also satisfies (102). Then by revers- 
ing the roles of y and % one obtains J(y)>I(y). 
This contradicts (103), and so implies that the as- 
sumption that y#7 satisfies (102) 1s false. 
cannot be a minimum. 

Tueorem 9. Let there exist finite numbers v, K>0, 
a®* such that if I (y)<v and 0<a<a*<1/K, formula 
89) holds 1 
and @ be chosen so that O<a<a<a*. Then if there 
erusts a constant k>O such that (97) is satisfied, the 
conclusions of theorem 8 are satisfied. Furthermore, 


there exists a function % such that for the sequence 


Let y® be such that I(y™)<v and let a@ | 


Let y he an ele ment of y satisfying 


Hence y 
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y} described in theorem 8 we have 


y°—7 ~=uniformly on a<zr<b 
y'—7' in the mean of order 2. 


The limit G affords a unique minimum to I(y) on the 
class of functions in Y joining the end points of y 

By lemma 18, J/(y) is bounded below on the class 
of y in Y joining the end points of y®. Since all of 
the other hypotheses of theorem 8 have been as- 
sumed, the conclusions of that theorem must be valid. 

As the sequence {y} of theorem 8 is such that 
each y‘ joins the end points of y® and 

Ty T(y) <v | eS 

formula (98) of lemma 1S is applicable to each term. 
This formula implies the existence of \/>0 such that 
for all n, m 


( 7° 7 ge ‘dx ) < M, 


(104) 


(n,m 


where 7 y"? — yy. For n>m, I(y =) -I(y ead 
is negative. Thus from substitution in the main 
formula of the proof of lemma 18 it follows that 


i 


0<k [n, m)’ pinem)’ dy < h mf ["s mem gen da )> 


where 


mi—?> @, 


h,, ( [rey ™)T, (z; 9") da ) 0. 


Consequently, the existence of 7 such that y'°’->7’ 
in the mean of order 2 follows from 104. This last 
result implies y“?-—>y uniformly on a<r2<h. By 
recalling (92) and taking subsequences one finally 
obtains a subsequence {y@ } such that 


y@—y uniformly on a<r <b 


yy?’ -—»7// almost uniformly 


T’ (2; y@)->@ almost uniformly. 
Consequently (by lemma 19) 7 satisfies (102). That 
7 is unique is a consequence of lemma 20. Hence the 
theorem is proved. It is not difficult to show that 
y (xr) is of class ce’. 

It is clear that these results apply to the case when 
f is a polynomial of degree two in y;. 


13. An Example of Numerical Computation 
Consider the problem of minimizing the integral 


I(y) | “yyi Ly de, y>0, 
ry 


whose Euler equation is 


» srt ‘ 
y J (105) 
Y 
To apply Newton’s method write (105) as 
T=1- y" yy” 0 
Then 
6T = 2y'n’ —y"n— yn" 
and the differential equation 7+67=0 can be put 


into the form 


n” =An+ Bn’ +C, 
where 
y" 21’ 1+ y’ 9 
Au™, B y, ¢ y 
i] Y 4 


The general solution of (105) is given by 


a 
4 


y b cosh : (106) 


Given two points (7,,y;) and (2, 42) three cases can 
occur 

(1) Two ecatenaries (106) join the 
In this case one will be minimizing and 
not 

(2) One catenary (106) joins the 
This catenary is not minimizing 

(3) No catenary (106) joins the 
The examples for computation were 
this in mind. 


given points. 
the other will 


given points. 
given points. 


selected with 


EXAmPLe I: 
(fo, Yo) are respectively 


The points (4, ¥%), 
(—2, 3.086) and (2, 3.086 The solutions are 
») J 
y=2 cosh 5 
z 
) 1.404 cosh ; 
’ 1.404 


where the first is minimizing and the second is not. 
The results of computation are presented below in 
tables 1, 2, and 3. In these tables as well in tables 
4, 5, y denotes the exact solution, y“ the initial ap- 
proximation, and y“ the approximation at the ith 


step. The symbol C; denotes the value of 
ee 
c -¥f 
Y F 
along the curve y“. Along a solution of (105) we 
have c=0. The function v, found in tables 2 to 5 is 


a solution of the equation é¢e=0, which results from 
taking the variation of the equation ¢,—0. The 
evaluation of v, can be arranged to be a by-product 
of the computation of y“*". It follows from the 
theory of conjugate points that y will be a minimizing 
arc if for ¢, sufficiently small rv, does not vanish more 
than once. If v, vanishes twice, y is not minimizing. 





It is of interest to note that convergence to a 
tion of the Euler equations that was not a mini) 


was obtained. This shows the necessity of ma 
additional assumptions (as was done in sectior 
in order to insure convergence to a minimum. 


EXxamp.e II: 


The points (z,, ¥%:), (2, y), are, respectiy 
(0, 2.5894) and (5, 5.9284). The solutions are 


. I 
y=2 cosh 


1.5806 


1.8883 cosh 5 , 
L.SSS3 


4 


The first are is minimizing and the second is not 
It will be noticed that although the two solutions 
are quite close to each other convergence to each ons 
was obtained. 


Exampce III: 


The points (2,,y,) and (x,y) are, respectively 
(0,2), (5,3). In this case there is no solution. After 


an initial estimate consisting of the line joining (0,2 


and (5,3), the second iteration resulted in points 
(x,y) lying outside of the domain in which y>0 


The material in this section was supplied to th: 
author by M. R. Hestenes, under whose auspices thi 
computations were carried out. 





TABLE | 
2 cosh 1/2 5 86 e<r<?2 
y Y] y 
0 086 2 688 2132 2015 1.064 2 000 m4 se) OS ft 
2 3.086 2447 2.140 208% 2004 2 010 iM OST ow oor “ 
| 43.085 24644 2167 20% 2.034 2 040 24 OST va ™* “ 
| + 6 3.088 2 407 2218 2106 2 O85 2 OO ” ne a2! LD 
| 
| +.8 3.0865 2.542 2.277 2.177 (2.157 2 162 24 u6 Oz on 
+1.0 3.0868 2600 2300 2 20 124 101 O25 o4 Ly 
+1.2 3.086 2672 248 2.383 $24 107 Os ona “ 
| +1.4 3.088 2 24 114 ES 3 0s ” 
| +1.6 3.088 a4 121 ana (ay 
| +1.8 3. One 324 128 Oos O08 “ 
| 
| +20 3.086 3.086 3.088 3.086 3. 086 Os, 24 136 042 OF; “ 
| 
TABLE 2 
| y=2cesh 2 y' 1.88643 2 
| z y y y y ¢ r 
| - 
| 0 1.886) 2038 200 2000 070 010 ool oD 
| +.2 1.808; 2063 2013 200 O86 m0 001 ” 
+.4 1a 2. 092 2.04083 2 040 O83 om ool 
+. 6 1.904 2. 140 2.003 2 1 034 009 ool 
+. 8 2. 078 2. 208 2. 165 2. 142 00s 009 001 #2 
+1.0 2. 186 2. 237 2 258 2 255 022 008 on « 
+1.2 2. 318 2. 407 2.373 2. 371 O55 0o7 ool 8 
+1.4 2.474 2. 40 2.512 2 510 OsSY 008 ool 74 
+1.6 2.6 2 2. 676 2. 67 124 O05 oo! ¢ 
+1.8 2. 858 2. 878 2 8A7 2 see 158 0038 ool 4s 
+20 3. OR 3. O86 +. O86 3. O86 191 02 ool 
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